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15.  Supplementary  Notes 


167/ Abstract  yj 

— pilhis  paper  reviews  ongoing  efforts  aimed  at  ^he  solution  of  basic  problems  deal- 
ing with  ozone  chemistry  and  atmospheric  mass  transports  in  representative  two- 
dimensional  (altitude  and  latitude)  photochemical  models  of  the  troposphere  and 
stratosphere.  Emphasis  is  placed  on  the  impact  of  recent  measurements  of  the  re- 
action rate  coefficients  for  H02’  + NO  -^OH  + NO2  and  HO2  + 03^  0H  + 2 0£>  on'v(a) 


the  balance  between  ozone  production  andcfestruction  on  both  ozone  concentrations 
and  columns  at  middle  northern  latitudes  during  summer  for  an  oxygen-hydrogen- 
nitrogen  -(0-H-N)> atmosphere;  (b)  the-- calculation  of.  NOx  concentrations  in  the  upper,, 
troposphere  which  are  important  for  the  effects  of  high-a.lti.tude  subsonic  flight 
on  the  ozone  column,'  and  (c)  estimates  of  the  tropospheric(NOx;dnd  OH 
concentrations  as  well  as  the  removal  rates  of  water-soluble “species  in  the 
lower  troposphere  from  considerations  of  the  tropospheric  ozone  and  carbon  monoxide 
budgets.  The  severe  limitations  of  arbitrary  extensions  of  available  2-D  parameter- 
izations  of  the  dynamics  to  an  upper  troposphere  and  lower  stratosphere  assumed  to 
be  perturbed  by  large  emissions  of  engine  effluents  are  brought  out  from  numerical 
solutions  of  the  primitive  equations  for  instantaneous  and  continuous  sources  of 
inert  material  in  the  lower  stratosphere.  The  paper  includes  “preliminary  2-D 
results  based  on  the  revised  ozone  chemistry  for  the  trend  of  subsonic  and  super- 
sonic flight  effects  on  ozone*  and  (b)  some  recommendations  for  future  work 
concerning  the  2-D  treatment  offohe  dynamics  for  the  perturbed  upper  troposphere  and 
lower  stratosphere . \ 
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This  paper  reviews  ongoing  efforts  aimed  at  the  solution 
of  basic  problems  dealing  with  ozone  chemistry  and  atmospheric 
mass  transports  in  representative  two-dimensional  (altitude  and 
latitude)  photochemical  models  of  the  troposphere  and  strato- 
sphere. Emphasis  is  placed  on  the  impact  of  recent  measurements 
of  the  reaction  rate  coefficients  for  HO^  + NO  ■*  OH  + NO^  and 
HO^  + 0^  -*•  OH  + 2 0^  on  (a)  the  balance  between  ozone  produc- 
tion and  destruction  on  both  ozone  concentrations  and  columns 
at  middle  northern  latitudes  during  summer  for  an  oxygen-hydrogen- 
nitrogen  (0-H-N)  atmosphere,  (b)  the  calculation  of  N0x  concen- 
trations in  the  upper  troposphere  which  are  important  for  the 
effects  of  high-altitude  subsonic  flight  on  the  ozone  column, 
and  (c)  estimates  of  the  tropospheric  N0X  and  CH  concentrations 
as  well  as  the  removal  rates  of  water-soluble  species  in  the 
lower  troposphere  from  considerations  of  the  tropospheric  ozone 
and  carbon  monoxide  budgets.  The  severe  limitations  of  arbi- 
trary extensions  of  available  2-D  parameterizations  of  the 
dynamics  to  an  upper  troposphere  and  lower  stratosphere  assumed 
to  be  perturbed  by  large  emissions  of  engine  effluents  are 
brought  out  from  numerical  solutions  of  the  primitive  ea.uations 
for  instantaneous  and  continuous  sources  of  inert  material  in 
the  lower  stratosphere.  The  paper  includes  (a)  preliminary  2-D 
results  based  on  the  revised  ozone  chemistry  for  the  trend  of 
subsonic  and  supersonic  flight  effects  on  ozone,  and  (b)  some 
recommendations  for  future  work  concerning  the  2-D  treatment  of 
the  dynamics  for  the  perturbed  upper  troposphere  and  lower 
stratosphere . 


iii 


ACKNOWLEDGMENTS 


The  author  is  indebted  to  several  modelers  of  atmospheric 
phenomena,  who  exhibited  their  scientific  spirit  both  in  regard 
to  open  access  to  their  numerical  data  and  self-criticism  of 
the  basic  assumptions  behind  their  data.  In  alphabetical  order, 
they  are  Dr.  Paul  J.  Crutzen  of  the  National  Center  for  Atmo- 
spheric Research,  Dr.  Jack  Fishman  of  the  Colorado  State  Univer- 
sity, Dr.  j.D.  Mahlman  of  the  Geophysical  Fluid  Dynamics  Labora- 
tory (NOAA),  and  Dr.  George  Widhopf  of  the  Aerospace  Corporation. 
Although  they  have  reviewed  a draft  of  this  paper  and  provided 
valuable  comments,  the  contents  of  this  paper  are  the  sole  re- 
sponsibility of  the  author. 

The  author  is  also  indebted  to  Dr.  G.D.  Robinson,  Prof. 

J.R.  Holton,  and  Dr.  R.C.  Oliver  for  valuable  comments  and 
constructive  criticism  of  a previous  draft  of  this  paper. 


CONTENTS 


Abstract  ill 

Acknowledgments  v 

Overview  and  Summary  S-l 

1.  INTRODUCTION  1 

1.1  Atmospheric  Characteristics  Relevant  5 

to  Aviation  Effects 

1.2  Mechanisms  for  the  Modeling  of  Aviation  13 

Effects 

1.2.1  Absorption  of  Solar  Radiation  13 

1.2.2  Rate  Coefficients  of  Chemical  14 

Kinetics  Processes 

1.2.3  Heterogeneous  Processes  for  the  15 

Tropospheric  Sinks 

1.2.4  Water  Vapor  (H2O)  Sources  and  Sinks  16 

1.2.5  Mass  Transports  in  the  Reference  and  16 

Perturbed  Tropospheres  and  Stratospheres 

1.3  Summary  of  Pending  Changes  in  the  Modeling  19 

of  Aviation  Effects  on  Ozone 

2.  TWO-DIMENSIONAL  MODELING  OF  AVIATION  EFFECTS  21 

ON  OZONE 

2.1  Results  and  Characteristics  of  the  1976  25 

Crutzen  Model 

2.2  Results  and  Characteristics  of  the  1976  35 

Widhopf  Model 

3.  IMPACT  OF  CHANGES  IN  OZONE  CHEMISTRY  ON  2-D  47 

MODELING 

3.1  Underdestruction  of  Ozone  with  1977  Chemistry  49 

3.2  Estimates  of  Tropospheric  N0X  and  H2O2  Rainout  55 

Sinks 

3.2.1  Assessment  of  Tropospheric  N0X  in  the  58 
Northern  Troposphere 

3.2.2  Ozone  Change  Due  to  Subsonics  with  62 

1977  Chemistry 

3.2.3  Heterogeneous  Removal  Rates  of  H2O2  67 

in  Tropospheric  Sinks 

vii 


3.3  Comparison  of  2-D  Model  Tropospheric  Results 

3.4  Critical  Chemical  Reactions  for  Aviation 
Effects 

ASSESSMENT  OP  THE  2-D  FORMULATION  OF  THE 
ATMOSPHERIC  DYNAMICS 

4.1  Basic  Assumptions  in  the  2-D  Parameterizations 
of  the  Dynamics 

4.1.1  Decoupling  of  Temperature  and  Dynamics 
from  the  Chemistry 

4.1.2  Hypothesis  of  Large-Scale  Diffusion 
in  the  Stratosphere 

4.2  Evaluation  of  the  Extension  of  the  Diffusion 
Hypothesis  to  the  Perturbed  Atmosphere 

4.2.1  Dynamic  Meteorology  of  Eddy  Transports 
in  the  Stratosphere 

4.2.2  Mass  Transports  of  Inert  Material  in 
the  Reference  and  Perturbed  Atmospheres 

4.3  Statistical-Dynamical  Formulation  of 
2-D  Dynamics 

REFERENCES 


I 


FIGURES 


1 . 

Approximate  estimates  of  gas  budgets  of 

the  world 

7 

• 

2. 

Chemical  mechanisms  in  the  stratosphere 
sphere 

and  tropo- 

9 

§ 

3. 

Ozone  column  in  reference  atmosphere  as  a function 
of  latitude  and  month  in  the  Northern  and  Southern 
Hemispheres 

17 

44.. 

Isolines  of  the  daily  average  of  UV  solar  radiation 
flux  at  a typical  wavelength  of  300.14  nm  during 
winter  (January)  and  summer  (July) 

26 

5 . 

Simulated  ozone  concentration  and  ozone 
growth  with  decreasing  altitude  for  the 
atmosphere  at  l)5°N  during  summer  (August 

column 

reference 

30) 

28 

6. 

Simulated  vertical  profiles  of  daylight. 

zonal 

33 

average  of  the  volume  mixing  ratio  of  indicated 
species  in  the  reference  atmosphere  at  l»5°N  during 
summer  (August  30)  by  the  Crutzen  model  with  1976 
chemistry 

7.  Relative  effects  of  individual  subsonic  or  super- 
sonic flight  on  the  atmospheric  ozone  column 
during  summer  (August  30)  after  six  years  of  fleet 
operations 

8.  Comparison  of  simulated  and  observed  ozone  mixing 
ratio  at  35°N  during  June 

9.  Simulated  vertical  profiles  of  214-hour,  zonal  Ji 2 

average  of  the  volume  mixing  ratio  of  indicated 
species  in  the  reference  atmosphere  at  llO°N  during 
summer  (August  15)  by  the  Widhopf  model  with  1976 
chemistry 

10.  Effects  of  a combined  fleet  of  subsonic  and  super-  it h 
sonic  aircraft  on  the  atmospheric  ozone  column 

during  summer  (June  15)  and  fall  (October  15)  after 
five  years  of  fleet  operations 

11.  Illustration  of  the  impact  of  the  Howard-Evenson  51 

reaction  rate  coefficient  on  the  simulated  ozone 

mixing  ratio  for  the  reference  atmosphere 


lx 


Assumed  tropospheric  N0X  profiles  for  sensitivity 
analysis  of  the  tropospheric  ozone  and  CO  budgets 
to  N0X  and  HO 

Comparison  of  N0X  calculations  in  the  troposphere  by 
2-D  models  utilizing  1976  (Crutzen  and  Widhopf)  and 
revised  (Widhopf  and  Fishman-Crutzen)  chemistries 

Illustration  of  H2O  parameterization  in  the  Crutzen 
2-D  model  at  a grid  point  ([TgpD 

Simulation  of  water  vapor  distribution  in  the  ref- 
erence atmosphere  by  the  Crutzen  model  as  a function 
of  latitude  in  the  Northern  Hemisphere 

Simulation  of  water  vapor  distribution  in  the 
reference  atmosphere  by  the  Widhopf  model  as  a 
function  of  latitude  in  the  Northern  Hemisphere 

Meridional  Circulation  as  a function  of  season 

Schematic  representation  of  the  slope  (a)  of  the 
mixing  length  path  necessary  for  countergradient 
eddy  fluxes  in  the  Prandtl  mixing  length  hypotheses 

Simulation  of  meridional  circulation  by  a three- 
dimensional  model  of  the  troposphere-stratosphere- 
mesosphere  (0  < z < 70  km)  utilizing  the  condition 
of  geostrophic  equilibrium  and  the  spectral  repre- 
sentation of  the  dependent  variables 

Simulation  of  6-lines  by  a three-dimensional  passive 
tracer  model  utilizing  finite  differences  of  the 
primitive  equations  for  an  instantaneous  injection 
of  material  at  36°N  and  65  mbar  during  January 

Stratospheric  tracer  mixing  ratio  as  a function  of 
latitude,  altitude  and  season  for  winter  (December) 
of  the  fourth  year  after  an  assumed  instantaneous 
release  of  tracer  material  at  35°N  and  65  mbar 

Interactions  at  the  injection  level  (65  mbar)  among 
the  eddy  and  meridional  motions  in  the  dispersion 
of  tracer  material  (Eq.  19)  in  the  stratosphere 
during  winter  and  summer 

Enhancement  of  tracer  mixing  ratio  as  a function  of 
latitude  and  altitude  during  winter  (December)  of 
the  fifth  year  for  a simulated  SST-like  continuous 
release  of  tracer  material  as  per  Eq.  20 

Enhancement  of  tracer  mixing  ratio  as  a function  of 
latitude  and  altitude  during  summer  (August)  of  the 
sixth  year  for  a simulated  SST-like  continuous  re- 
lease of  tracer  material  as  per  Eq.  20 


■ 


126 


25.  Eddy  transports  and  the  local  rate  of  change  of 
mixing  ratio  (Eq.  19)  at  65  mbar  during  winter 
(December)  of  the  fourth  year  after  an  instantaneous 
release  of  passive  tracer  material 

26.  Eddy  transports  and  local  rate  of  change  of  mixing  128 
ratio  (Eq.  19)  at  65  mbar  during  winter  (December) 
after  5 years  of  an  SST-like  continuous  release  of 
passive  tracer  material 


1.  Mechanisms  (M)  for  production  and  destruction  10 

of  chemical  species 


2a.  Daylight,  zonal  averages  of  ozone,  NgO,  N0X,  30 

N2O5,  NO3  and  HNO3  species  in  the  reference 
atmosphere  as  was  simulated  by  the  2-D  Crutzen 
photochemical  model  with  1976  chemistry  at  45°N 
during  late  summer  (August  30) 


2b. 

Daylight,  zonal  averages  of  H2O,  H2O2J  OH,  CO 
and  CHi}  species  in  the  reference  atmosphere  as 
was  simulated  by  the  2-D  Crutzen  photochemical 
model  with  1976  chemistry  at  45°N  during  late 
summer  (August  30) 

31 

3. 

1976  reaction  rate  coefficients  for  an  oxygen- 
hydrogen-nitrogen  atmosphere 

36 

4. 

24-hour,  zonal  average  of  species  in  the  reference 
atmosphere  as  simulated  by  Widhopf,  et  al.  (1977) 
at  40°N  during  late  summer  (August  15) 

41 

5. 

Ratio  of  values  based  on  Howard-Evenson  rate 
(daily  averages)  to  those  corresponding  in 

Table  4 (24-hour  average)  at  40°N  during  late 
summer 

53 

6. 

Sensitivity  of  ozone  budget  to  prescribed  N0X 
concentrations  for  the  Northern  Hemisphere 

61 

7. 

• 

Sensitivity  of  CO  budget  to  prescribed  N0X 
concentrations  for  the  Northern  Hemisphere 

61 

8. 

Sensitivity  of  tropospheric  ozone  production 
to  prescribed  N0X  concentrations  for  the 

Northern  Hemisphere 

64 

9. 

Competition  of  reactions  involving  the  HO2 
radical  for  important  subsonic  flight  conditions 

66 

10. 

Sensitivity  of  OH  to  removal  rates  of  H2O2  for 
the  assumed  N0X  profiles  A and  B in  the  Northern 
Hemisphere 

68 

11. 

Average  residence  time  of  water  vapor  in  the 
troposphere  as  a function  of  latitude 

74 

xiii 


80 


I 

f 


12.  Sensitivity  of  the  sign  of  the  ozone  change 
to  rate  coefficients  of  R21  and  R13 

13.  Sample  eddy  parameters  for  Crutzen  model  during  99 

summer  and  winter  in  Northern  Hemisphere 

14.  Sample  eddy  parameters  for  Widhopf  models  during  100 
summer  (August)  in  Northern  Hemisphere  and  winter 

in  Southern  Hemisphere 

15.  Sample  eddy  parameters  derived  from  poleward  103 

transport  of  ozone 

16.  Primitive  equations  for  meridional  circulation  108 

17.  Primitive  equations  for  eddy  motions  109 

18.  Slope  of  mixing  length  path  as  given  by  129 

primitive  equations,  [a]  (10-^) 

19.  Vertical  eddy  coefficients  as  given  by  primitive  131 
equations,  Kzz  (10^  cm^/sec) 

20.  Horizontal  eddy  coefficients  as  given  by  132 

primitive  equations,  (10 10  cm^/sec) 


xiv 


OVERVIEW  AND  SUMMARY 


This  paper  is  limited  to  a review  of  the  status,  as  of 
mid-1978,  of  ongoing  efforts  aimed  at  the  solution  of  basic 
problems  which  are  involved  in  the  development  of  two-dimen- 
sional (altitude  and  latitude)  photochemical  models  of  the 
stratosphere  and  troposphere.  The  interest  in  these  two- 
dimensional  (2-D)  photochemical  models  stems  from  the  following 
factors : 

• A need  to  forecast  the  probable  effect  of  future  world- 
wide, high-altitude  traffic  of  subsonic  and  supersonic 
aircraft  on  atmospheric  ozone  as  a function  of  altitude, 
latitude,  and  season.  High-altitude  flight  can  have 

an  impact  on  ozone  through  the  interactions  of  NC>x  (NO 
and  NO 2 ) and  H20  engine  effluents  with  a rather  large 
family  of  trace  chemical  species  that  affect  either 
directly  or  indirectly  the  ozone  distributions  in  the 
stratosphere  and  troposphere. 

• A need  for  proper  interpretation  of  field  observations 
of  trace  chemical  species  affecting  ozone  as  a function 
of  altitude,  latitude,  and  season.  These  results  are 
important  for  the  gaining  of  physical  insights  for 
both  photochemical  modeling  of  the  atmosphere  and  plan- 
ning subsequent  field  observations. 

The  scope  of  this  paper,  therefore,  excludes  considerations 
of  either  the  more  limited  1-D  models  (altitude  only)  or  the 
currently  less  practical  3-D  models  (Including  longitude)  for 
the  assessment  of  the  latitudinal  and  seasonal  effects  on  ozone 
from  high-altitude  flight.  Furthermore,  this  paper  describes 


the  ongoing  development  of  2-P  modeling  by  considering  in  some 
depth  two  representative  models  of  special  interest  to  the  High 
Altitude  Pollution  Program  (HAPP),  i.e.,  the  Crutzen  (Hidalgo 
and  Crutzen,  1977)  and  the  Widhopf  (Widhopf  et  al.,  1977,  1978) 
models . 

The  basic  problems  in  the  development  of  atmospheric  photo- 
chemical models  have  involved  the  numerical  simulation  of  com- 
plex couplings  among  atmospheric  radiation,  chemistry,  and 
dynamics  in  the  stratosphere  and  troposphere.  The  complexity 
of  these  phenomena , coupled  with  the  absence  of  adequate , rele- 
vant , laboratory , and  field  data , as  well  as  the  cost  limitations 
of  available  computers  have  imposed  severe  constraints  on  the 
development  of  photochemical  models.  A review  of  the  development 
of  representative  2-D  photochemical  models  has  become  important 
due  to  recent  drastic  changes  in  the  understanding  of  the  ozone 
chemistry  as  a result  of  a direct  measurement  of  the  reaction 
rate  for  HO^  + NO  -*■  OH  + NO.,  by  Howard  and  Evenson  (1977). 

A summary  of  2-D  developments  as  of  mid-1978  is  given  below 
for  (a)  a reference  atmosphere  defined  as  one  without  anthropo- 
genic effects  on  ozone  from  either  aircraft  effluents  or  surface 
emissions  of  chlorof luoromethanes . This  latter  condition  has 
implied  omission  of  the  more  recent  chlorine  chemistry  in  2-D 
models,  (b)  the  formulation  of  the  dynamics  for  an  atmosphere 
assumed  to  be  perturbed  by  large  emissions  of  aircraft  effluents, 
and  (c)  the  probable  trends  of  aviation  effects  on  ozone  due  to 
the  revised  understanding  of  the  atmospheric  chemistry.  This 
summary  of  2-D  developments  is  then  followed  by  some  recommended 
areas  for  future  work. 

| 

S.l  REFERENCE  ATMOSPHERE 

Significant  efforts  have  been  concerned  with  the  problem 
of  calculating  excessive  ozone  in  the  reference  stratosphere 
from  use  of  the  Howard-Evenson  (1977)  rate  coefficient  for  HO., 


S-2 


+ NO  -*•  OH  + NO^.  A 2-D  model  (WIdhopf  and  Glatt,  1978),  utiliz- 
ing this  1977  rate  coefficient  along  with  other  1977  revisions 
of  the  chemistry,  cannot  calculate  enough  ozone  destruction  to 
obtain  a balance  between  ozone  production  and  destruction  that 
matches  the  observed  ozone  in  the  stratosphere.  Results  from 
this  model,  based  on  a revised  1977  chemistry,  yield  (a) 
ozone  concentrations  that  are  larger  than  those  based  on  pre- 
vious (1976)  chemistry  by  as  much  as  28  percent  at  35  km  alti- 
tude at  middle  northern  latitudes  during  summer  (Fig.  11  and 
Table  5),  (b)  ozone  columns  for  the  combined  stratosphere  and 
troposphere  that  were  from  20  to  30  percent  larger  than  either 
observations  or  those  obtained  with  previous  (1976)  chemistry, 
and  (c)  tropospheric  ozone  columns,  representing  about  10  per- 
cent of  the  combined  tropospheric  and  stratospheric  columns, 
that  are  significantly  higher  than  the  corresponding  values 
based  on  1976  chemistry  (Tables  4 and  5),  but  still  lower  than 
observations.  Subsequent  1-D  model  results  incorporating  a 1978 
measurement  of  the  rate  coefficient  for  H02  + 0^  -*•  OH  + 2 0^ 
(Howard,  1978),  indicated  a decrease  (by  about  9 percent)  of  the 
reference  1-D  ozone  column  as  a consequence  of  using  the  1978 
instead  of  the  1977  rate  coefficient  for  HO,,  + 0,  -+  OH  + 2 02 
(Zahniser  and  Howard,  1978).  Similar  considerations  in  the 
Widhopf  2-D  model  indicate  that  the  reference  ozone  column  at 
middle  northern  latitudes  during  summer  also  decreases  when 
using  the  1978  instead  of  the  1977  rate  coefficient  for  H0o  + 

0^  + OH  + 2 0o  in  combination  with  the  1977  rate  coefficient 
for  HO^  + NO  -*■  OH  + NO,,;  however,  this  later  2-D  result  repre- 
sents an  ozone  column  that  is  still  about  10  percent  higher 
than  the  corresponding  observations  at  middle  latitudes  durir.r 
summer.  Furthermore,  this  2-D  ozone  column  decrease  is  caused 
mainly  by  a corresponding  ozone  decrease  in  the  lower  strato- 
sphere; i.e.,  the  calculated  ozone  imbalance  in  the  middle 
stratosphere  appears  to  be  sensitive  mainly  to  the  1977  rate 
coefficient  of  HO,,  + NO  -*■  OH  + NO.,  in  an  oxygen-hydrogen-nitrogen 
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( O-H-N ) atmosphere  (Fig.  11).  Therefore,  this  preliminary  2-D 
result  suggests  a need  for  a reassessment  of  the  remaining  re- 
action rate  coefficients  for  an  O-H-N  atmosphere  as  well  as  the 
consideration  of  chlorine  chemistry  in  the  middle  and  upper 
stratosphere . 

The  use  of  the  Howard-Evenson  (1977)  rate  coefficient  for 
HO^  + NO  -*■  OH  + NO^  in  the  Widhopf  2-D  model,  together  with  1977 
chemistry  and  a revised  rainout  sink  decreases  the  calculated. 

N0X  in  the  upper  troposphere  by  more  than  an  order  of  magnitude 
from  the  corresponding  value  based  on  1976  chemistry.  This 
result  indicates,  then,  calculated  N0x  concentrations  in  the 
upper  troposphere  that  are  in  far  better  agreement  than  before 
with  available  limited  observations  (Noxon,  1978).  A correct 
calculation  of  the  NO  background  in  the  upper  troposphere  is 
of  particular  interest  for  the  effect  of  subsonic  flight  on  the 
ozone  column. 

Estimates  of  the  N0x  concentrations  in  the  troposphere 
have  been  made  by  Fishman  and  Crutzen  (1978)  using  a 2-D 
tropospheric,  nondynamic,  chemical  model.  These  estimates 
are  based  mainly  on  (a)  the  sensitivity  of  the  tropospheric 
ozone  balance  between  photochemical  production  and  destruc- 
tion to  N0x  amount  (Table  6),  (b)  the  sensitivity  of  the 
tropospheric  CO  budget  to  NO  amount  through  the  effect  of  NO 
on  OH  (Table  7),  and  (c)  a revised  chemistry  and  radiation  as 
given  by  the  Howard-Evenson  (1977)  rate  coefficient  for  HO^  + 

NO  -*•  OH  + NO2,  a third  body  component  for  the  rate  coefficient 
of  CO  + OH  + M -*■  H + CO2  + M,  the  NASA  (1977)  recommended 
chemistry,  and  the  temperature  dependence  of  the  quantum  yield 
in  the  photolysis  of  02  + hv  -*■  0(1D)  + O2.  The  results  of  these 
sensitivity  calculations  yield  NO  concentrations  in  the  upper 
troposphere  at  middle  northern  latitudes  that  are  in  rather 
good  agreement. with  the  foregoing  Widhopf  2-D  model  results 
based  on  1977  chemistry  (Fig.  13;  Widhopf,  1978)  as  well  as  the 
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Noxon  observations . This  sensitivity  study  of  the  CO  budget 
also  provides  estimates  for  the  unknown  heterogenous  removal 
rate  of  water-soluble  species  in  the  troposphere.  It  is  thus 
estimated  that  the  inverse  heterogenous  removal  rate  constant 
of  H2O2  in  the  troposphere  is,  at  most,  of  the  order  of  1A  hr 
(Table  10).  Finally,  it  is  important  to  note  that  the  use  of 
the  1977  rate  coefficient  for  HO  2 + NO  -*•  OH  + NOp  gives  more 
reasonable  values  than  before  for  the  photochemical  sources  and 
sinks  in  the  CO  budget. 


S . 2 DYNAMICS  OF  THE  PERTURBED  ATMOSPHERE 

The  present  extension  of  arbitrary  2-D  parameterizations 
of  the  essentially  unknown  dynamics  in  the  reference  strato- 
sphere to  describe  the  mass  transports  of  effluents  in  the 
stratosphere  perturbed  by  aircraft  effects  has  no  scientific 
basis.  The  present  parameterizations  of  the  2-D  dynamics  for 
the  reference  stratosphere  are  based  on  two  arbitrary  assump- 
tions (Reed  and  German,  1965):  (a)  uncoupling  of  the  mass 

transports  by  the  meridional  (2-D)  circulation  from  the  eddy 
mass  transports  by  macroscale  turbulence,  which  is  not  justified 
by  the  equations  of  motion  (Tables  16,  17);  and  (b)  adaptation 
of  the  Prandtl  mixing  length  hypothesis  for  the  eddy  transport 
of  momentum  in  microscale  flow  to  describe  the  mass  transport 
by  macroscale  turbulence  in  the  atmosphere.  The  adaptation  of 
this  hypothesis  to  the  stratosphere  is  based  on  assumptions  con- 
cerning the  unobservable  magnitude  of  a mixing  length  path  angle 
(a)  defined  by  the  ratio  of  the  vertical  to  the  horizontal  eddy 
transports  (Fig.  18).  The  use  of  the  basic  3-D  conservation 
equations  for  the  mass  transport  from  an  instantaneous  release 
of  a passive  tracer  in  the  stratosphere  indicates  that  the  mag- 
nitudes of  the  mixing  length  path  angles  are  not  as  given  in  the 
2-D  parameterizations  of  the  stratospheric  dynamics  (Tables  13 
through  15  and  18).  The  mixing  length  path  angles  from  the 
conservation  equations,  when  coupled  with  the  mixing  length 


hypothesis,  would  then  lead  to  negative  values  of  the  eddy  dif- 
fusion coefficients  for  the  reference  stratosphere,  which  are 
physically  unacceptable  (Tables  19,  20).  Furthermore,  the  rela 
tive  behavior  of  the  a-parameter  is  not  the  same  in  the  refer- 
ence and  perturbed  upper  tropospheres  as  well  as  stratospheres, 
as  indicated  by  the  3-D  solutions  for  the  mass  transport  from 
a continuous  release  (i.e.,  SST  or  supersonic  transport-like 
sources)  of  a passive  tracer  (Table  18) . The  significance  of 
these  results  is  as  follows: 

• Any  empirical  parameterization  of  the  2-D  dynamics 
for  the  reference  stratosphere  based  on  appropriate 
observations  would  indeed  be  acceptable,  even  if  they 
were  to  utilize  concepts  other  than  the  mixing  length 
approach.  Two  relevant  points  are  as  follows:  (a) 
such  observations  have  been  difficult  to  achieve  due 
to  the  smallness  of  the  vertical  transports  in  the 
stratosphere  (Oort  and  Rasmusson,  1971),  and  (b)  the 
available  2-D  parameterizations  for  the  reference 
atmosphere,  based  on  the  mixing  length  concept  are  not 
internally  consistent  at  every  latitude  in  the  strato- 
sphere and  upper  troposphere  (Section  4.1.2).  Hence, 
just  as  photochemical  models  identify  their  use  of 
reaction  rate  coefficients,  they  must  also  document 
their  use  of  eddy  parameters  in  the  lower  atmosphere 
(0  < z < 50  km) . 

• The  extension  of  any  empirical  2-D  dynamics  for  the 
reference  atmosphere  to  the  significantly  enhanced 
mass  eddy  transports  of  engine  effluents  in  the  flight 
corridors  of  the  perturbed  upper  troposphere  and  lower 
stratosphere  is  not  justified  by  the  conservation 
equations  (Figs.  25  vs  26,  Tables  18  through  20),  which 
are  the  only  means  to  describe  the  eddy  transports  of 
mass  in  the  atmosphere  perturbed  by  aircraft  effluents. 
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It  is  important  to  note  that  the  ozone  perturbations 
due  to  the  distribution  of  NC  effluents  in  the  upper 
troposphere,  as  calculated  from  the  1976  Crutzen  2-D 
model  of  the  dynamics  and  chemistry  (Hidalgo  and  Crutzen, 
1977),  are  expected  to  increase  due  to  subsequent  revi- 
sions in  the  chemistry  alone  (Section  3.2.2).  Hence, 
it  also  becomes  important  to  revise  the  1976  extensions 
of  the  dynamics  of  the  reference  atmosphere  to  the 
perturbed  conditions. 

S . 3 AVIATION  EFFECTS 

Reliable  assessments  of  aviation  effects  as  a function  of 
projected  air  traffic,  latitude,  and  season  for  subsonic  and/or 
supersonic  high-altitude  flight  are  not  available  because  2-D 
models  have  not  yet  incorporated  validated  (a)  chemistry  for 
an  O-H-N  and  chlorine  atmosphere  and  (b)  dynamics  for  the  per- 
turbed upper  troposphere  and  lower  stratosphere.  Therefore, 
the  available  preliminary  2-D  results  may  be  suggestive  more 
of  a trend  than  absolute  assessments  of  probable  aviation  effects 
on  atmospheric  ozone. 

Previous  2-D  estimates  based  on  1976  chemistry  by  Hidalgo 
and  Crutzen  (1977)  of  the  ozone  column  increase  due  to  NO^ 
emissions  from  an  assumed  1990  traffic  of  subsonic  aircraft  in 
the  upper  troposphere  (Fig.  7)  will  become  larger  from  the  use 
of  the  Howard-Evenson  (1977)  measurement  of  the  rate  coefficient 
for  HO.,  + NO  -*■  OH  + N02  (Table  9).  Furthermore,  the  previous 
2-D  estimate  based  on  1976  chemistry  by  Widhopf,  et  al.  (1977) 
of  the  peak  ozone  column  increase  by  NO  emissions  from  a .1990 
combined  fleet  of  subsonic  and  supersonic  transports  (Fig.  10) 
becomes  larger  by  a factor  of  about  2.5  from  the  use  of  (a)  the 
1977  rate  coefficient  for  H02  + NO  -+  OH  + NO?,  (b)  the  NASA 
(1977)  recommended  rates  for  other  reactions,  (c)  a modified 
dynamics  for  the  middle  and  upper  stratosphere,  (d)  the  simul- 
taneous emissions  of  N0x  and  H,,0  effluents,  and  (e)  a larger 
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SST  fleet  by  a factor  of  3 (Widhopf  and  Glatt,  1978).  However, 
this  result  for  the  combined  fleet  is  indicative  more  of  the 
direction  than  the  magnitude  of  the  effect  of  aircraft  effluents 
on  the  ozone  column  from  the  1977  revisions  of  the  1976  ozone 
chemistry.  This  is  so  because  of  the  following  constraints  on 
this  result:  (a)  use  of  the  1977  instead  of  the  1978  rate  co- 
efficient for  H0o  + 0^  -*■  OH  + 2 0? , (b)  use  of  the  1977  (none) 
instead  of  the  1978  temperature  dependence  for  the  rate  coeffi- 
cient of  H0n  + NO  ■+  OH  + NO.-,,  and  (c)  the  calculated  imbalance 
between  ozone  production  and  destruction  in  the  reference  strato- 
sphere (Fig.  11).  Nevertheless,  the  significance  of  these  pre- 
liminary considerations  for  individual  subsonic  (Table  9)  and 
combined  subsonic  and  supersonic  fleets  is  to  make  the  modeling 
of  the  upper  troposphere  for  aviation  effects  more  important 
than  ever  before. 

Furthermore,  the  absolute  magnitude  of  current  2-D  fore- 
casts of  aviation  effects  on  the  ozone  column  are  inconclusive 
due  to  the  following  factor's:  (a)  lack  of  experimental  verifi- 
cation of  the  1978  rate  coefficients  in  several  chemical  re- 
actions for  the  effects  of  subsonic  and  supersonic  cruise  flight 
in  the  upper  troposphere  and  lower  stratosphere,  respectively; 

(b)  lack  of  adequate  data  to  determine  the  tropospheric  budgets 
of  species  such  as  NO  , OH,  and  CO.  It  is  important  to  note 
that  measurements  of  species  with  a short  lifetime  (e.g.,  OH) 
require  simultaneous  measurements  of  related  species  (e.g.,  H^O, 
0^ , N0x,  Ho0^,  CO,  and  CH^)  to  permit  proper  Interpretation  of 
their  variability;  and  (c)  lack  of  a validated  dynamics  for  the 
perturbed  upper  troposphere  and  lower  stratosphere  as  a result 
of  extensions  of  arbitrary  parameterizations  of  the  dynamics  of 
the  reference  atmosphere  to  an  atmosphere  perturbed  by  assumed 
large  emissions  of  engine  effluents.  The  foregoing  factors 
become  important  for  the  forecast  of  small  aviation  effects  on 
atmospheric  ozone,  which  are  the  result  of  the  difference 
between  two  large  numbers;  i.e.,  of  the  ozone  columns  of  the 
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perturbed  and  reference  (without  aircraft  effects)  atmospheres. 
The  absolute  macnitude  of  forecasts  of  larger  aviation  (or  any 
other  anthropogenic)  effects  on  the  ozone  column  are  also  open 
to  question,  because  current  2-D  (or  1-D)  photochemical  models 
cannot  account  for  feedback  effects  on  ozone  changes  from 
changes  in  the  temperature  and  dynamics  of  the  perturbed  strato- 


sphere . 

S . 4 RECOMMENDED  AREAS  FOR  FUTURE  WORK 

Although  significant  attention  continues  to  be  given  to 
the  experimental  measurement  of  critical  reaction  rates  for  the 
assessment  of  anthropogenic  effects  on  atmospheric  ozone,  there 
are  other  problem  areas  that  merit  consideration,  as  described 
below . 

A sensitivity  analysis,  based  on  1-D  models,  is  needed  to 
establish  the  relative  importance  on  the  ozone  effects  from 
aviation  of  the  uncertainties  in  the  knowledge  of  reaction  rates 
and  reaction  mechanisms  in  the  ti’opospheric  ozone  chemistry. 

This  sensitivity  analysis  is  becoming  of  interest  because  2-D 
results,  based  on  1977  revisions  of  the  ozone  chemistry,  indi- 
cate an  increased  importance  of  the  effect  of  subsonic  flight 
in  the  upper  troposphere  on  the  ozone  column. 

The  extension  of  the  2-D  parameterizations  of  the  dynamics 
for  the  reference  atmosphere  to  the  dynamics  perturbed  by  assumed 
large  emissions  of  engine  effluents  in  the  upper  troposphere  and 
lower  stratosphere  should  be  evaluated  by  comparing  2-D  distri- 
butions of  passive  tracers  with  corresponding  zonal  averages  of 
3-D  numerical  solutions  for  aircraft-like  sources  of  such  tracers 
in  the  upper  troposphere  and  lower  stratosphere,  some  of  which 
are  beginning  to  become  available  from  the  fiFDL  model  (Mahlman, 
1978)  . 


Statistical-dynamical  2-D  formulations  of  the  dynamics  are 
needed  for  the  transports  of  mass  in  the  perturbed  upper  tropo- 
sphere and  stratosphere  for  assumed  large  emissions  of  engine 
effluents  from  high-altitude  subsonic  and  supersonic  flight. 

Such  statistical-dynamical  approaches,  when  applied  to  the 
troposphere,  have  been  able  to  simulate  the  three-cell  structure 
of  the  troposphere  (Fig.  19).  The  development  of  a statistical- 
dynamical  2-D  model  for  the  stratosphere  would  avoid  the  use  of 
the  Prandtl  mixing  length  hypothesis,  but  will  stress  the  current 
meteorological  understanding  of  the  stratosphere.  Because  of 
the  relative  difficulty  of  this  endeavor,  an  early  advent  of  a 
valid  2-D  statistical-dynamical  formulation  of  the  stratospheric 
dynamics  will  also  depend  on  the  priority  set  for  such  endeavor 
by  dynamic  meteorologists. 


1.  INTRODUCTION 


During  the  last  few  years,  a great  deal  of  effort  has  been 
devoted  to  understanding  the  potential  impact  of  high-altitude 
aircraft  flights  on  the  atmospheric  ozone  column.  The  interest 
in  ozone  stems  from  its  strong  absorption  of  UV-B  solar  radia- 
tion at  wavelengths  (X)  in  the  range  290  < X < 320  nm.  A de- 
crease of  the  ozone  column  by  aircraft  engine  effluents,  for 
example,  would  increase  the  UV-B  radiation  flux  at  the  ground; 
an  effect  that  could  have  complex  ramifications  on  the  biosphere 
(Grobecker  et  al.,  1974  ; National  Academy  of  Sciences,  1975). 

Plight  regimes  of  current  interest  include  future  subsonic 
and  supersonic  transports  (SSTs)  operating  at  cruise  altitudes 
(z)  between  6 < z < 14  km  and  15  < z < 19  km,  respectively  (e.g., 
Oliver  et  al.,  1977).  The  mechanisms  by  which  high-altitude 
flight  would  have  an  impact  on  atmospheric  ozone  are  (1)  the 
deposition  of  N0x  and  water  vapor  (H^O)  engine  effluents  in  the 
flight  corridors,  where  N0X  stands  for  both  nitric  oxide  (MO) 
and  nitrogen  dioxide  (I^);  (2)  the  eventual  global  dispersion 
by  air  motions  of  such  engine  effluents  throughout  the  lower 
atmosphere  (0  < z < 50  km),  and  (3)  the  photochemical  and  chemi- 
cal kinetics  reactions  of  the  aircraft  NO  and  H„0  effluents 

X c 

with  species  that  control  directly  or  indirectly  the  formation 
and  destruction  of  atmospheric  ozone.  While  any  detrimental 
effect  of  high-altitude  flights  on  atmospheric  ozone  could,  of 
course,  be  reduced  through  worldwide  regulations  on  the  NO  and 
H^O  engine  emissions,  the  assessment  of  the  effect  of  engine 
effluents  on  ozone  is  essential  in  order  to  (a)  prevent  any 
costly,  unnecessary  regulation;  and  (b)  persuade  governments,  on 
an  international  basis,  of  the  need  for  any  necessary  regulations. 
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The  complexity  of  atmospheric  phenomena  involving  UV  solar 
radiation,  photochemistry  as  well  as  chemical  kinetics,  and  air 
motions  represents  a very  difficult  obstacle  to  the  goal  of 
achieving  an  accurate  assessment  of  the  effects  of  aircraft  ef- 
fluents on  atmospheric  ozone.  While  the  overall  scientific 
methodology  for  the  ozone  problem  (relative  to  that  for  climate 
change)  appears  to  be  in  hand  (e.g.,  Hidalgo,  ±978),  the  sources 
of  basic  difficulties  are,  in  general,  as  follows: 

• Lack  of  adequate  laboratory  data  to  help  develop 
complete  and  valid  models  of  the  photochemistry  and 
gas-phase  chemical  kinetics  in  important  regions  of 
the  lower  atmosphere  (0  < z < 50  km). 

• Lack  of  adequate  laboratory  data  to  establish  the 
heterogeneous  chemistry  of  water-soluble  trace  species 
or  sinks  for  such  species  in  the  lower  troposphere. 
These  sinks  are  important  because  they  can  have  sig- 
nificant effects,  for  example,  on  the  long-term  con- 
centrations of  NOx  in  the  lower  atmosphere. 

• Lack  of  adequate  atmospheric  statistics  to  define  the 
global  distributions  of  relevant  trace  chemical  species 
in  the  lower  atmosphere  as  a function  of  latitude, 
altitude,  and  season.  These  statistical  data  are  impor 
tant  to  help  model  (a)  the  chemistry  by  providing  con- 
straints on  ratios  and  budgets  of  atmospheric  chemical 
species,  (b)  the  transports  of  such  species  by  air 
motions,  and  (c)  the  variability  of  their  distributions 
due  to  the  variability  of  both  the  UV-B  radiation  and 
air  motions,  for  any  given  month  or  season  during  suc- 
cessive annual  cycles. 

• The  limitations  in  current  computer  capabilities  (i.e., 
speed  and  cost)  when  attempting  to  perform  global  cal- 
culations of  the  lower  atmosphere  with  an  without  air- 
craft effluents. 


2 


W'h- 


A forceful  appreciation  of  the  foregoing  constraints  can 
perhaps  be  obtained  by  considering  the  following  additional 
factors:  (a)  the  smallness  of  the  probable  threshold  value 

for  permissible  (i.e.,  unregulated)  decreases  in  the  ozone 
column  due  to  NO^  and  H^O  aircraft  effluents  (perhaps  0.5  per- 
cent, Grobecker  et  al.,  197^),  (b)  the  calculated  small  changes 
in  the  ozone  column  by  aircraft  effluents,  changes  that  are  the 
result  of  the  difference  between  two  large  numbers;  i.e.,  the 
ozone  column  perturbed  by  the  aircraft  effluents  and  a refer- 
ence ozone  column  that  is  unperturbed  by  such  effluents.  It 
should  be  noted  that  this  factor  is  of  particular  interest 
whenever  the  temperature  field  and  mass  transports  in  the  per- 
turbed atmosphere  are  actually  different  from  those  in  the 
reference  atmosphere,  and  (c)  the  long-term  transient  nature 
of  the  reference  ozone  column  without  aircraft  effects,  since 
the  reference  ozone  column  is  a variable  that  depends  on  other 
time-dependent  natural  (i.e.,  variability)  and  anthropogenic 
effects.  The  latter  stem  from  emissions  at  the  ground  level 
of  chlorof luoromethanes  (National  Academy  of  Sciences,  1976); 
nitrous  oxide  resulting  from  use  of  industrial  fixed  nitrogen 
fertilizers  (Crutzen,  1976;  McElroy  et  al.,  1976;  Liu,  1977); 
and  carbon  monoxide  (CO),  which  interacts,  ultimately,  with 
methane  (CH^),  and  water  vapor  in  the  stratosphere  (Crutzen 
and  Fishman,  1977). 

It  is  therefore  to  be  expected  that  the  foregoing  four 
types  of  basic  constraints  have  brought  about  severe  limitation 
on  the  development  of  adequate  photochemical  models  of  the  low- 
er atmosphere  (0  < z < 50  km).  For  example,  under  the  time 
constraints  of  the  recent  Climatic  Impact  Assessment  Program 
(CIAP),  it  was  only  possible  to  develop  a 3-D  model  based  on 
(a)  low  geographical  resolution  (~  1500  km),  (b)  an  auxiliary 
2-D  model  for  the  distribution  of  NO  effluents,  and  (c)  very 
restricted  (now  obsolete)  197^  chemistry  (Cunnold  et  al.,  1975» 
1977;  Hidalgo,  1976).  In  the  same  time  scale  of  CIAP,  it  was 


also  possible  to  develop  a 3-D  model  with  a high  geographical 
resolution  (~  250  km)  for  the  study  of  the  dispersion  of  a 
chemically  inert  tracer  in  the  lower  atmosphere  (z  < 31  km), 
a tracer  that  was  assumed  to  dissipate  ultimately  only  in  the 
rainout  sinks  of  the  lower  atmosphere  (Mahlman,  1973,  1975). 

To  circumvent  the  severe  restrictions  brought  about  by 
computer  limitations  on  the  considerations  of  the  complex 
chemistry,  it  also  became  necessary  to  consider,  during  CIAP, 
the  development  of  2-D  and  1-D  photochemical  models  (Grobecker 
et  al.,  197^ ; Hidalgo,  1976).  The  development  of  2-D  models 
was  based  on  considerations  of  the  time  scales  for  the  chemistry 
and  air  motions,  which  suggested  that  the  longitudinal  effects 
may  be  disregarded,  as  a first  approximation,  for  aviation 
effects  on  ozone  in  time  scales  significantly  longer  than  the 
diurnal  cycle  introduced  by  the  earth's  rotation.  The  2-D 
photochemical  models,  then,  had  to  describe  the  radiation, 
chemistry,  and  air  motions  in  the  lower  atmosphere  in  terms 
of  longitudinal  or  zonal  averages  of  appropriate  variables 
at  each  latitude  and  altitude.  The  further  development  of  the 
more  practical  2-D  models  in  a time  scale  subsequent  to  CIAP, 
has  continued  to  put  emphasis  on  the  couplings  among  radiation, 
chemistry,  and  dynamics  in  the  lower  atmosphere. 

A main  objective  of  this  paper  is  therefore  to  describe 
the  status,  as  of  mid- 1 97 8 , of  the  development  of  representa- 
tive two-dimensional  photochemical  models  of  interest  to  the 
current  High  Altitude  Pollution  Program  (HAPP)  of  the  Federal 
Aviation  Administration  (FAA).  As  indicated  previously,  atmos- 
pheric photochemical  models  are  indispensable  tools  to  simulate 
the  lower  atmosphere  (0  < z < 50  km)  with  and  without  aircraft 
effluents.  An  output  of  these  2-D  models  is  then  the  assessment 
of  the  potential  effects,  subject  to  the  four  basic  constraints 
noted  before  of  subsonic  and  supersonic  flights  on  the  ozone 
column  as  a function  of  latitude,  season,  and  both  the  location 
and  intensity  of  present  and  near-future  world  airline  traffic. 
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A review  of  the  development  of  2-D  photochemical  model  has 
become  important  at  this  time  because  of  recent  drastic  changes 
in  the  description  of  the  relevant  ozone  chemistry . Important 
review  considerations  are  also  the  several  arbitrary  parameteri- 
zations  for  the  essentially  unknown  dynamics  in  the  reference- 
stratosphere  and  the  arbitrary  extensions  of  such  paraineteri- 
zations  to  the  stratosphere  and  upper  troposphere  perturbed  by 
assumed  large  emissions  of  engine  effluents. 

The  scope  of  this  paper  excludes  direct  consideration  not 
only  of  the  3-D  but  also  of  the  1-D  (altitude  only)  photochemi- 
cal models,  since  the  latter  are  useful  only  to  provide  esti- 
mates of  aviation  effects  on  the  global  and  annual  averages  of 
the  ozone  column  for  a globally  and  seasonally  uniform  air 
traffic  at  given  flight  altitudes.* 

1.1  ATMOSPHERIC  CHARACTERISTICS  RELEVANT  TO  AVIATION  EFFECTS 

A clearer  understanding  of  the  basic  problems  in  HAPP  re- 
quires a brief  review  of  the  relevant  atmospheric  characteris- 
tics of  the  lower  atmosphere  between  the  ground  and  50  km 
altitude.  Basic  characteristics  of  interest  are  the  variation 
of  temperature  with  altitude  as  well  as  the  trace  chemical 
gaseous  species  that  can  interact  directly  or  indirectly  with 
atmospheric  ozone. 

As  described  in  the  rather  extensive  body  of  literature 
(e.g.,  CIAP  Monograph  1),  the  lower  atmosphere  between  0 and 
50  km  consists  of  the  troposphere  and  stratosphere . The  tro- 
posphere is  the  lowest  atmospheric  layer  that  extends  from  the 


1-D  photochemical  models  have,  however,  been  a powerful 
tool  for  examining  the  sensitivity  of  aviation  effects 
to  the  uncertain  knowledge  of  the  chemistry. 
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from  use  of  the  Howard-Evenson  (1977)  rate  coefficient  for  H02 
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ground  to  the  tropopause.  The  location  of  the  tropopause  Is 
at  about  16  km  at  low  latitudes,  0,  (i.e.,  30°N  < 0 < 30°S) 
which  decreases  to  about  half  as  much  at  the  colder  high  lati- 
tudes. The  troposphere  is  characterized  by  a negative  vertical 
temperature  gradient,  which  tends  to  vanish  at  the  tropopause. 

The  layer  immediately  above  the  tropopause  is  the  stratosphere, 
which  extends  upwards  from  the  tropopause  to  the  stratopause 
at  about  50  km  altitude.  The  stratosphere  is  characterized  by 
a negligible  vertical  temperature  gradient  in  the  lower  stra- 
tosphere (z^  < z 20  km,  where  z denotes  the  tropopause 
altitude)  and  a positive  vertical  temperature  gradient  at  the 
middle  and  upper  stratosphere.  The  positive  vertical  gradient 
tends  to  vanish  at  the  stratopause.  It  is  important  to  indi- 
cate that  the  vertical  temperature  gradients  in  the  stratosphere 
are  produced  by  the  absorption  of  UV  solar  radiation  by  ozone. 
Therefore,  the  altitude  regime  for  subsonic  cruise  (6  < z < 1^  km) 
extends  to  the  upper  troposphere  and  even  lower  stratosphere  at 
middle  latitudes,  whereas  that  for  SSTs  (15  < z < 19  km)  is 
located  in  the  lower  stratosphere. 

Figure  1,  after  Crutzen  (1978),  summarizes  schematically 

the  atmospheric  chemical  species  of  main  interest  to  HAPP. 

This  figure  shows  the  origin  as  well  as  the  estimated  world's 

budgets  for  the  relevant  gases  in  the  ozone  (0^)  chemistry; 

i.e.,  nitrous  oxide  (N20) , oxides  of  nitrogen  (N0x),  water 

vapor  (H20),  methane  (CH^),  carbon  monoxide  (CO),  hydrogen  (H2), 

methyl  chloride  (CH^Cl),  and  fluorocarbons  (CFCl^,  CF2C10). 

The  arrows  in  the  figure  indicate  whether  a given  chemical  is 

being  removed  from  or  released  to  the  atmosphere.  The  numbers 

at  the  ground  level  are  their  respective  fluxes  in  molecules 
2 

per  cm  sec.  The  numbers  in  the  troposphere  are  estimated 
volume-mixing  ratios  (i.e.,  concentrations  relative  to  air) 
whereas  the  percent  values  given  at  the  tropopause  indicate 
the  fraction  that  reaches  the  stratosphere.  The  percent  values 
at  the  bottom  of  the  figure  indicate  the  relative  anthropogenic 
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source  for  each  of  the  indicated  species.  The  bent  arrows 
touching  the  ground  indicate  the  effect  of  the  tropospheric 
sinks  for  such  species. 

Figure  2,  also  after  Crutzen  (1978),  indicates  the  mecha- 
nisms or  chemical  reactions  that  either  create  (inward  arrows) 
or  destroy  (outward  arrows)  particular  chemical  species  in  the 
stratosphere  and  troposphere.  These  mechanisms  are  given  in 
Table  1,  which  shows  the  photochemical  and  chemical  kinetics 
processes  from  the  absorption  of  solar  radiation  (hv)  at 
several  wavelengths  (X)  and  collisional  processes  among  chemical 
species,  respectively.  The  mechanisms  in  Table  1 reflect  the 
understanding  of  the  ozone-related  chemistry  as  of  mid-1978.# 

In  the  case  of  ozone,  for  example,  arrow  1 indicates  the  crea- 
tion of  ozone  (mainly  in  the  tropical  middle  stratosphere)  by 
photolysis  or  absorption  of  UV  radiation  (Ml,  Table  1);  whereas 
arrow  Ma,  b)  indicates  the  creation  of  ozone  in  the  troposphere 
from  the  oxidation  of  CO  and  CH^  (MiJa,  b).  Similarly,  arrow  2 
indicates  the  catalytic  destruction  of  ozone  in  the  stratosphere 
by  nitric  oxide  (NO),  atomic  chlorine  (Cl),  hydrogen  (H), 
hydroxyl  radical  (OH),  and  atomic  bromine  (Br).  Arrow  3 (a,  b, 
c)  indicates  the  destruction  of  ozone  in  the  troposphere  by 
metastable  atomic  oxygen  [0(1D)],  H,0,  the  catalytic  OH-HO^ 
cycle,  and  that  of  NO.  and  NO,. 

Figures  1 and  2 indicate  that  a significant  source  for  N0x 
in  the  reference  stratosphere  is  given  by  the  surface  emissions 
of  NpO  (which  are  not  well  known),  its  inert  propagation  from 
the  ground  to  the  stratosphere,  and  mechanism  6 (Table  1)  in 
the  stratosphere.  Other  sources  of  N0x  in  the  reference  atmo- 
sphere are  anthropogenic  sources  at  the  ground  as  well  as 


Because  of  rather  rapid  changes  in  the  understanding  of 
the  chemistry,  a time  reference  (given  in  years)  has  become 
necessary  in  the  study  of  anthropogenic  effects  on  atmo- 
spheric ozone. 
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TABLE  1 

MECHANISMS  (M)  FOR  PRODUCTION  AND 

DESTRUCTION  OF  CHEMICAL 

SPECIES  (Source:  Crutzen, 

1978) 

Ml 

02  + hv  -►  2 0 

o + o2  + m-o3  + m 

X < 240  nm 

3 02  - 2 03 

M2 

03  + hv  -*•  0 + Og 

X = NO,  Cl,  H,  OH,  Br 

- 

o + xo  - x + o2 

X + 03  - X0  + 02 

2 °3  - 3 °2 

M3a 

03  + hv  - 0( 7D)  + 02 

O^D)  + H20  -►  2 OH 

X < 310  nm 

03  + H20  - 2 OH  + 02 

M3b 

OH  + 03  - H02  + 02 

H02  + 03  - OH  + 2 02 

2 °3  - 3 °2 

M3c 

NO  + 03  - N02  + 02 

no2  + o3  - no3  + o2 

N03  + hv  - NO  + 02 

2 °3  - 3 °2 


M4a  CO  + OH  + M - H + C02  + M 

H + 02  + M - H02  + M 
H02  + NO  - OH  + N02 

N02  + hv  -►  NO  + 0 X < 400  nm 

Q + 02  + M-»-03  + M 

CO  + 2 02  -*•  C02  + 03 


10 


_____ 


taneous  emissions  of  N0x  and  H„0  effluents,  and  (e)  a larger 
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TABLE  1.  (cont'd) 


M4b*  CH4  + OH  - CH3  + H20 

CH3  + 02  + M - CH302  + M 
CH302  + NO  - CHjO  + N02 

ch3o  + o2  -»  ch2o  + ho2 

CH20  + hv  - H + HCO  As  350  nm 

H + 02  + M - H02  + M 
HCO  + 02  - H02  + CO 
H02  + NO  ♦ OH  + N02 

N02  + hv  -*■  NO  + 0 A < 400  nm 

o + o2  + m-o3  + m 

CH4  + 8 02  - 4 03  + CO  + H20  + 2 OH 

M5  N20  + 0^0)  -*•  N2  + 02 

N20  + hv  - N2  + 0 A < 337  nm 

M6  N20  + 0(]D)  - 2N0 


M7  X + e,  hv  -*•  X*  + 2e,  e 

X = n2,  o2,  0 

(Particles,  precipitation,  protons,  electrons, 
aurorae,  extreme  UV,  X-rays)  solar,  galactic, 
cosmic  rays. 

M8  NO  + hv  - N + 0 A < 191  nm 

N + NO  - N2  + 0 

2N0  - N2  + 2 0 

M9  N02  + aerosol,  hydrometeors  -*■  products 

N02  + 0H(+  M)  - HN03(+  M) 

HN03  + aerosol , hydrometeors  + products 

M10  Lightning,  F » 1 (10),  cm2/s 

Mila  CH4  +0H  - -*■  CO,  H2  - C02,  H?0 

*The  net  result  in  M4b  is  obtained  by  considering  multiples  of  the  last 
three  reactions;  i.e. , 4 times  the  last  two,  and  3 times  the  remaining  one. 
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TABLE  1.  (cont'd) 


Ml  lb 

CH4  + 0 ( 1 D ) -....-  CO,  H2  - C02, 

H20 

M12 

CH,0  + hv 

- CO  + H2 

X 

2 

M13 

C02  + hv 

- CO  + 0 

M14 

CO  + 0H(+  M) 

-*  C02  + H(+  M) 

Ml  5a 

H2  + 0 ( 1 D ) 

- H + OH 

M15b 

H + H02 

- H2  + °2 

M16 

H,  + OH 

- H + H,0 

z 

2 

M17* 

CH3C1  + OH 

-<■  CH2C1  + H20. . 

C1X 

X s 350  nm 


CFC1 3 + hv  - CFC12  + Cl....  - 3C1X 
CF2C1 2 + hv  - CF2C1  + Cl....  - 2C1 X 
CH3CC13  + OH  - CH2CC13  + H20..-+  3C1X 
C1X  = Cl  + CIO  + HC1  + Cl N03 

(In  troposphere:  rain,  washout,  aerosol  attachment.) 


Crutzen,  Isaksen  and  McAfee  (1978) 
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natural  sources  (M7,  10),  and  sinks  (M8,  9),  In  the  troposphere 
and  stratosphere.  Figures  1 and  2 also  indicate  that  the  strat- 
osphere is  a relatively  dry  medium  with  natural  sources  pro- 

vided by  tropical  upwelling  (as  a function  of  longitude)  of  tro- 
pospheric moisture  as  well  as  reactions  of  CH^  with  OH  and  0(1D) 
(i.e.,  Mila,  b).  The  fluorocarbons  (CFCl^  and  CF2C12)  are  of 
particular  importance  in  a reference  atmosphere  without  aviation 
effects  alone,  because  of  their  chemical  destruction  of  strato- 
spheric ozone  (M17,  2)  and  their  interaction  with  the  NO^  air- 
craft effluents. 


1.2  MECHANISMS  FOR  THE  MODELING  OF  AVIATION  EFFECTS 

Because  of  the  four  types  of  basic  constraints  in  the  assess- 
ment of  aviation  effects  on  the  ozone  column,  as  indicated  pre- 
viously, a difficult  challenge  has  been  that  of  achieving  a rath- 
er accurate  modeling  of  the  global  distributions  of  the  chemical 
species  shown  in  Figs.  1 and  2 for  a reference  troposphere  and 
stratosphere  without  aircraft  effluents  and  emissions  of  fluoro- 
carbons. The  main  mechanisms  that  must  be  considered  in  such 
modeling,  as  well  as  some  illustration  of  the  uncertainties  in- 
volved in  each  mechanism,  are  as  follows: 

1.2.1  Absorption  of  Solar  Radiation 

The  absorption  of  solar  radiation  (hv)  by  oxygen,  ozone, 
nitrogen  dioxide,  nitric  acid,  etc.  takes  place  at  different 
wavelengths  (X),  as  indicated  in  Table  1.  Parameters  of  inter- 
est here  are  the  photolysis  rates  (as  defined  subsequently  in 
Section  3.2)  for  mechanisms  (M)  such  as  those  in  Ml,  M2,  and  M3a. 

An  example  of  the  uncertain  knowledge  of  some  of  these  photolysis 
rates  is  given  by  the  reaction  0^  + hv  -*•  0('*'D)  + 0^  in  M3a;  re- 
cent laboratory  measurements  (Arnold  et  al.,  1977;  Moortgat  et 
al.,  1977)  have  shown  that  the  quantum  yield  (defined  also  in 
Section  3.2)  of  metastable  atomic  oxygen,  0(1D),  for  ozone  photo- 
lysis between  300  and  320  nm  decreases  with  decreasing  temperature. 
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Previous  calculations  had  used  quantum  yield  data  found  at  room 
temperature  (e.g.,  Moortgat  and  Warneck,  1975). 

1.2.2  Rate  Coefficients  of  Chemical  Kinetics  Processes 

The  rate  coefficients  of  the  chemical  kinetics  processes 
(denoted  by  k's)  establish  the  speed  or  rate  of  change  of  the 
concentrations  (molecules/cm  ) of  given  chemical  species  in 
either  the  two  or  three  (involving  the  M species)  body  reactions 
in  Table  1.  The  uncertain  knowledge  of  these  rate  coefficients 
can  be  illustrated  by  recent  measurements  of  these  coefficients 
for  (a)  H02  + NO  -*  OH  + N02  (e.g.,  M4a,  b),  (b)  H02  + 0^  OH  + 

2 02  (e.g.,  M3b),  and  (c ) CO  + OH  + M -+■  H + C02  + M (e.g.,  M4a). 
Recent  developments  for  these  reactions  are  as  follows: 

• A direct  laboratory  measurement  of  the  rate  coefficient 
for  H0?  + NO  OH  + N0~  yielded  a value  of  k = (8.1  + 

1.5)  10  cm'Vmolecule-sec#  at  a temperature  of  296°K 
(Howard  and  Evenson,  1977).  This  result  turned  out  to 
be  higher  by  a factor  of  about  40  than  the  previous 
accepted  value  (2.2  x 10  J),  which  was  used  in  the 

1976  assessments  of  the  effect  of  high-altitude  flight 
on  the  ozone  column  (Hidalgo  and  Crutzen,  1977;  Widhopf 
et  al.,  1977).  This  1977  result  for  the  rate  coefficient 
of  the  H02  + NO  -*•  OH  + N02  reaction  has  had  a profound 
impact  on  the  modeling  of  the  reference  ozone  column  as 
a function  of  latitude  and  season.  It  is  important  to 
indicate  that  the  Howard-Evenson  result  for  the  rate  co- 
efficient of  H02  + NO  -*■  OH  + N02  has  been  verified  at 
other  laboratories  (Burrows  et  al.,  1978;  Kaufman  and 
Reimann,  1978;  Margitan  and  Anderson,  1978).  More  recent 
developments  concerning  this  rate  coefficient  have  yielded 
its  temperature  dependence  as  k = (3.3  ± 0 . 7) 10"12exp[ ( 254 
±50)/T]  cm^/molecule-sec  (Howard,  1978).  A preliminary 

3 —12  3 —1  —1 

i.e.,  k = (8.1  ± 1.5)10  " cnr  molecule-  sec-  , which  has  been 

broken  up  as  shown  for  editorial  reasons. 

14 


result  of  the  temperature  dependence  for  this  rate  co- 
efficient in  the  range  270  s T £ 425°K  has  also  yielded 
k = (4.6  ± 1 . 3 ) 10-12  exp  [(200  ± 90)/T]  cm3/molecule-sec , 
where  the  uncertainties  represent  double  standard  devi- 
ations (Ming-Taun  Leu,  1978). 

• A more  recent  measurement  of  the  rate  coefficient  for 

H02  + 03  - OH  + 2 02  yielded  k = (1.4  ± 0.4)10-l4exp 

[-(580  ± 100)/T]  cm3/molecule-sec  (Howard,  1978).  This 

result  is  faster  than  the  NASA  (1977)  recommended  rate  1 

-]  4 

coefficients  (7-3  x 10  exp  -1275/T)  by  factors  of 
3.3  to  4.5  at  the  (standard)  temperature  range  of  the 
lower  (~220°K)  and  middle  (~240°K)  stratosphere. 

• For  the  rate  coefficient  of  the  CO  + OH  -*•  C02  + H re- 
action, evidence  has  been  found  recently  of  a three- 
body  component  (in  addition  to  the  two-body  component) 
that  becomes  important  at  the  high  pressure  of  the  tro- 
posphere (Sie  et  al.,  1976;  Cox  et  al.,  1976;  and  Chan 
et  al. , 1977)  • 

1.2.3  Heterogeneous  Processes  for  the  Tropospheric  Sinks 

Heterogeneous  chemical  processes  determine  the  removal  rate 
of  water-soluble  species  such  as  NO  , nitric  acid  (HNO.,),  hydro- 

* J 

gen  peroxide  (H202),  nitrogen  pentoxlde  (N20^),  etc.  from  the 
troposphere,  and  become  important  in  the  modeling  of  their  global 
distributions.  The  N0x  distributions  in  the  reference  tropo- 
sphere and  stratosphere  are  important  in  establishing  the  mag- 
nitudes of  (a)  the  perturbations  of  atmospheric  NOx  by  the  cor- 
responding aircraft  effluents,  and  (b)  the  consequent  changes 
in  the  reference  ozone  column  by  such  effluents.  The  uncertainty 
in  the  magnitude  of  the  observed  NOx  concentrations  in  the  tro- 
posphere is  rather  large,  because  of  the  variability  of  polluted 
and  unpolluted  conditions  at  the  earth's  surface  for  different 
geographical  locations.  Nitric  acid  is  an  important  sink  for 
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N0X  through  mechanism  9 (Table  1),  whereas  hydrogen  peroxide  is 
important  for  the  global  distributions  in  the  troposphere  of 
species  such  as  OH,  CH^ , and  CO  (Fishman  and  Crutzen,  1978). 
Therefore,  the  modeling  of  the  heterogeneous  processes  in  the 
tropospheric  sinks  becomes  important  for  the  assessment  of  air- 
craft effects  on  the  ozone  column. 

1.2.4  Water  Vapor  (H?0)  Sources  and  Sinks 

The  water  vapor  sources  and  sinks  in  the  troposphere  and 
stratosphere  are  important  in  determining  the  H20  distributions 
in  the  reference  stratosphere,  which  are  in  turn  important  to 
establish  (a)  the  magnitude  of  the  perturbations  of  H^O  by  the 
corresponding  engine  effluents,  and  (b)  the  consequent  changes 
in  the  reference  ozone  column  by  such  effluents.  The  uncertainty 
in  the  magnitude  of  the  observed  H^O  concentrations  in  the  strato- 
sphere is  significant,  because  of  the  uncertain  knowledge  of  the 
tropospheric  sources  and  sinks  or  the  transport  of  H20  across 
the  tropopause  due  to  (a)  the  temperature  variations  as  a func- 
tion of  longitude  in  the  upwelllng  branch  of  the  Hadley  cell  at 
low  latitudes,  and  (b)  the  subsidence  or  descending  transports 
as  a function  of  latitude. 

1.2.5  Mass  Transports  in  the  Reference  and  Perturbed  Tropo- 
spheres and  Stratospheres 

The  importance  of  the  mass  transports  by  the  air  motions 
in  the  stratosphere,  for  example,  comes  about  from  the  facts 
that  (a)  they  export  ozone  from  the  production  regions  at  low 
latitudes  to  the  northern  high  latitudes.  This  mechanism  be- 
comes most  evident  during  the  polar  night  of  the  winter  hemi- 
sphere, where  the  ozone  column  is  about  twice  as  large  as  that 
at  the  equatorial  latitudes  (Fig  3,  Diitsch,  1971);  and  (b)  they 
likewise  establish  the  stratospheric  water  vapor  distributions, 
which  are  a result  of  the  H20  transports  across  the  tropopause 
or  tropospheric  sources  and  sinks  in  addition  to  Mll(a,  b)  in 
Fig.  2.  The  prevailing  practice  in  the  2-D  parameterization  of 
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FIGURE  3.  Ozone  column  in  reference  atmosphere  as  a function  of  latitude 
and  month  in  the  Northern  and  Southern  Hemispheres.  Numbers 
are  in  Dobson  units  or  m-atm-cm;  i.e.,  if  the  ozone  column 
would  be  compressed  to  standard  pressure  and  temperature,  its 
height  at  90°N  during  spring,  for  example,  would  be  0.440  cm. 
Note  the  asymmetry  in  the  ozone  columns  in  the  Northern  and 
Southern  Hemispheres.  (Source:  Dutsch,  1971) 


the  stratospheric  dynamics  has  been  strictly  empirical  and  based 
on  extensions  of  the  mixing  length  hypothesis  for  microscale 
flow  to  the  macroscale  turbulence  in  the  atmosphere.  The  use 
of  this  hypothesis  brings  about  a set  of  eddy  coefficients, 
which  must  be  determined  from  either  the  matching  of  the  distri- 
butions of  water  vapor  and  ozone  in  the  lower  stratosphere  or 
extrapolations  of  tropospheric  transports  of  heat  as  supported 
by  the  matching  of  fragmentary  observations  of  the  dispersion 
of  nuclear  debris  in  the  stratosphere.  The  characteristics  of 
the  current  2-D  parameterizations  of  the  reference  stratospheric 
dynamics  may  be  judged  from  the  following  factors:  (a)  the 
theoretically  unsound  decoupling  of  the  advection  and  eddy  mo- 
tions. The  former  represent  longitudinal  averages  of  the  north- 
ward and  vertical  winds,  whereas  the  latter  represent  the  depar- 
ture from  such  averages;  (b)  the  consequent  combination  for  the 
same  atmosphere  of  a finite  set  of  different  statistics  for  the 
eddy  motions  (as  given  by  different  2-D  models)  with  the  same 
advection  motion;  and  (c)  the  fact  that  the  1977  change  in  the 
reaction  rate  coefficient  for  the  HO^  + NO  -*■  OH  + N02  reaction 
as  described  earlier,  might  also  require  readjustments  of  the 
previous  parameterizations  of  the  eddy  motions  based  on  the 
ozone  and  water  vapor  distributions.  It  must  be  noted  that 
while  a given  empirical  fit  of  a set  of  observations  may  be 
adequate  to  describe  such  data,  it  is  also  true  that  arbitrary 
empirical  fits  of  observations  may  not  necessarily  be  adequate 
for  extensions  of  such  data  to  the  atmosphere  perturbed,  for  air- 
craft effluents.  This  latter  argument  is  consistent  with  one  of 
the  constraints  indicated  earlier;  i.e.,  the  calculated  small 
changes  of  the  ozone  column  by  aviation  are  given  by  a difference 
between  two  large  ozone  columns  for  the  reference  and  perturbed 
atmospheric  conditions. 


1.3  SUMMARY  OF  PENDING  CHANGES  IN  THE  MODELING  OF  AVIATION 

EFFECTS  ON  OZONE 

The  foregoing  considerations  imply  a need  for  significant 
pending  revisions  in  the  2-D  modeling  of  aviation  effects  on 
ozone  as  a consequence  of  results  from  ongoing  research  in  (a) 
the  chemical  mechanisms  in  Table  1 and  (b)  the  numerical  solu- 
tions of  the  3-D  equations  of  motion  for  instantaneous  or  con- 
tinuous releases  of  passive  tracers  in  the  lower  stratosphere, 
to  be  described  later.  The  basic  pending  revisions  as  of  mid- 
1978  of  representative  2-D  photochemical  modeling  of  the  tropo- 
sphere and  stratosphere  are  as  follows: 

• Incorporation  of  chlorine  chemistry  in  the  stratosphere, 
i.e.,  mechanisms  M17,  2 (with  x * Cl)  in  Table  1;  which 
becomes  important  in  defining  a reference  atmosphere  for 
aviation  effects  that  takes  into  account  the  effects 
from  emissions  of  fluorocarbons  at  the  earth's  surface 
(Fig.  2).  The  use  of  rate  coefficients  for  the  chlorine 
chemistry  must  consider  the  NASA  (1977)  recommended 
values  together  with  results  of  ongoing  research  as 
they  become  available.  However,  any  significant  ozone 
decrease  from  fluorocarbons  would  also  require  consid- 
eration of  feedback  effects  from  the  ozone  decrease  on 
the  temperature  and  vertical  transports  of  a 2-D  refer- 
ence stratosphere  with  fluorocarbons. 

■ • Use  of  the  Howard-Evenson  (1978)  rate  coefficients  for 
(a)  the  reaction  H0o  + NO  -*  OH  + N0o,  including  the 

C. 

temperature  dependence  which  has  been  measured  recently, 
and  (b)  the  reaction  H0o  + 0,  -*■  OH  + 2 0_. 

• Use  of  the  third  body  component  for  the  rate  coefficient 
of  the  reaction  CO  + OH  + M ■+■  CO,,  + H + M as  well  as  the 
temperature  dependence  for  the  quantum  yield  in  the  photo- 
lysis rate  of  0^  + hv  -*•  0(1D)  + 0^- 
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• Validation  of  the  modeling  of  the  water  vapor  distri- 
butions in  the  reference  stratosphere,  which  is  impor- 
tant for  the  assessment  of  the  simultaneous  effects  of 
N0x  and  Ho0  engine  effluents  on  the  ozone  column.  While 
parameterizations  of  the  Ho0  distributions  in  the  ref- 
erence atmosphere  are  available,  the  understanding  of 
the  actual  Ho0  budget  in  the  stratosphere  is  not  yet 
conclusive . 

• Other  basic  considerations  involve  the  2-D  modeling  of 
the  dynamics  in  the  upper  troposphere  and  stratosphere, 
which  ought  to  consider  the  development  of  more  basic 
approaches  that  do  not  rely  on  utilizations  of  the  mix- 
ing length  hypothesis  for  microscale  flow  to  describe 
the  macroscale  turbulence  of  the  atmosphere.  A physi- 
cally correct  formulation  of  the  dynamics  in  the  upper 
troposphere  and  stratosphere  is  of  special  interest  for 
applications  to  the  perturbed  atmosphere. 

Subsequent  sections  attempt  to  elaborate  on  the  foregoing 
pending  revisions  of  the  2-D  modeling  of  the  troposphere  and 
stratosphere.  These  sections  include  (a)  a brief  review  of  the 
most  recent  2-D  model  results  for  HAPP,  which  are  based  on  1976 
chemistry,  (b)  the  impact  of  recent  revisions  of  the  ozone 
chemistry  on  2-D  modeling,  and  (c)  assessment  of  the  2-D  for- 
mulation of  the  atmospheric  dynamics  for  the  atmosphere  perturbed 
by  aircraft  effluents. 


2.  TWO-DIMENSIONAL  MODELING  OF  AVIATION  EFFECTS  ON  OZONE 


The  first  task  in  the  2-D  photochemical  modeling  of  aviation 
effects  on  the  reference  ozone  column  i_  to  attempt  to  simulate 
numerically,  with  a fair  degree  of  precision,  the  meridional 
(i.e.,  latitude-altitude)  distributions  of  the  important  chemical 
species  for  the  ozone  chemistry  in  the  reference  troposphere  and 
stratosphere  (e.g..  Figs.  1 and  2).  This  task  involves  the  use 
of  the  radiative,  chemical,  and  dynamical  mechanisms  for  the 
production  and  destruction  as  well  as  transport  by  air  motions 
not  only  of  ozone,  but  also  of  every  chemical  constituent  that 
interacts  either  directly  or  indirectly  with  ozone.  The  produc- 
tion and  destruction  mechanisms  consist  of  the  gas-phase  chemical 
reactions  given  in  Table  1 and  the  heterogenous  chemical  pro- 
cesses in  the  tropospheric  sinks  (Fig.  1),  while  the  transports 
by  air  motions  must  take  into  account  those  by  the  advection  and 
eddy  motions.  The  accuracy  of  the  numerical  simulations  of  the 
reference  troposphere  and  stratosphere  may  be  assessed  from  com- 
parisons between  the  calculated  and  observed  vertical  profiles 
as  a function  of  latitude  and  season  for  as  many  species  as  data 
are  available. 

The  second  task  in  the  2-D  modeling  of  aviation  effects  on 
the  reference  ozone  column  is  to  prescribe  the  aircraft  NO^  and 
Ho0  emissions  in  the  flight  corridors  of  the  upper  troposphere 
and  lower  stratosphere.  The  results  of  these  considerations  are 
the  determination  of  the  aircraft  NO  and  H-,0  effluents  as  a 

A <_ 

function  of  altitude  and  latitude  for  each  type  of  aircraft. 

These  emissions  can,  of  course,  also  be  a function  of  time  (years) 
as  given  by  projections  of  future  changes  in  the  structure  of  the 
world's  airline  fleets  and  traffic. 

There  are  two  2-D  photochemical  models  of  special  Interest 
to  HAPP,  because  they  have  utilized  aircraft  emissions  of  interest 
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to  the  FAA . They  are  the  Crutzen  model  at  the  National  Center 
for  Atmospheric  Research  (NCAR)  and  the  Widhopf  model  at  the 
Aerospace  Corporation.  The  Crutzen  model  has  been  used  to  study 
the  effect  of  individual  subsonic  or  supersonic  fleets,  assumed 
to  operate  by  themselves  at  cruise  altitudes  in  the  range  of 
10.8  £ z s 18  km,  so  as  to  isolate  the  effect  of  flight  altitude 
and  aircraft  type  on  the  total  ozone  column.  The  Widhopf  model 
has  been  used  to  study  the  effect  of  a combined  fleet  of  subsonic 
and  supersonic  aircraft  operating  simultaneously  at  altitudes 
in  the  range  6 £ z £ 19  km.  The  need  for  at  least  two  different 
2-D  photochemical  models  stems  from  their  complexity,  since  such 
results  help  to  identify  model-dependent  effects  which  are  un- 
related to  the  aircraft  effluents. 

A unique  characteristic  of  the  foregoing  2-D  models  was 
their  incorporation  of  the  ozone  production  mechanisms  in  the 
troposphere  due  to  CO  and  CH^  reactions  with  OH  (i.e.,  M^a,  b. 
Table  1).  As  a consequence  of  these  mechanisms,  results  from  the 
Crutzen  model  indicated  that  subsonic  flight  in  the  upper  tro- 
posphere may  produce  a small  increase  of  the  reference  ozone 
column  at  the  northern  middle  latitudes  of  the  assumed  heaviest 
airline  traffic  (Hidalgo  and  Crutzen,  1977).  These  results  were 
considered  tentative  in  view  of  the  uncertainty  in  some  of  the 
1976  reaction  rate  coefficients  as  well  as  in  the  heterogenous 
removal  rates  of  chemical  species  in  the  rainout  sinks  of  the 
lower  troposphere.  The  magnitude  of  the  increase  in  the  refer- 
ence ozone  column  by  subsonic  flight  was  apparently  large  enough 
to  provide  a net  small  increase  of  the  ozone  column  for  a com- 
bination of  a large  subsonic  and  a small  supersonic  aircraft 
fleet  (Widhopf  et  al.,  1977).  The  results  from  these  2-D  models 
contradicted  previous  assessments  of  the  net  effect  of  high- 
altitude  aircraft  which  had  been  based  on  the  use  of  1-D  global 
models  that  neglected  the  ozone  production  mechanisms  in  the 
troposphere  (Grobecker  et  al.,  197^;  National  Academy  of  Sciences, 
1975). 
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As  indicated  earlier,  the  foregoing  results  from  the  Crutzen 
and  Widhopf  models  did  not  incorporate  chlorine  chemistry  (i.e., 
M17,  2,  Table  1)  in  the  reference  stratosphere.  A question  of 
interest  for  aircraft  effects  was  then  the  relative  roles  of  the 
tropospheric  ozone  production  (M4a,  b)  and  stratospheric  ozone 
destruction  (M2,  with  x = NO  and  Cl  = Br  = 0)  mechanisms  as  in- 
dicated in  Fig.  2.  Thus,  even  for  subsonic  flight  alone  (e.g., 
at  10.8  km),  these  two  mechanisms  would  come  into  play  as  the 
aircraft  N0x  effluents  propagated  eventually  from  the  upper 
troposphere  to  the  middle  stratosphere.  The  Crutzen  model  in- 
dicated that  the  relative  roles  of  these  two  mechanisms  depended 
on  the  flight  altitude  of  aircraft  fleets  assumed  to  operate 
individually  at  a given  cruise  altitude.  It  was  found  that  the 
CO  and  CH^  mechanisms  (M^a,  b)  in  the  troposphere  would  prevail 
over  that  of  the  N0x  catalytic  cycle  (M2)  in  the  stratosphere 
for  subsonic  flight  at  10.8  km,  but  that  the  opposite  would  be 
true  in  the  operation  of  SSTs  alone  at  18  km  altitude  (Hidalgo 
and  Crutzen,  1977).  For  combined  instead  of  individual  aircraft 
fleets,  the  CO  and  CH^  mechanisms  in  the  troposphere  would  still 
dominate  over  that  of  the  N0x  catalytic  cycle  for  assumed  condi- 
tions of  a rather  overwhelming  subsonic  traffic  relative  to  that 
of  SSTs  (Widhopf  et  al.,  1977). 

It  must  be  emphasized  that  the  foregoing  Crutzen  model 
results  for  the  individual  aircraft  fleets  were  based  on  the 
1976  estimate  of  the  rate  coefficient  for  the  reaction  HO^  + NO  -*• 
OH  + NO^,  an  estimate  that  turned  out  to  be  smaller  (i.e.,  slower) 
by  a factor  of  about  40  than  the  value  determined  subsequently  by 
direct  measurements  of  HO^  (Howard  and  Evenson,  1977).  The  use 
of  the  slower  rate  coefficient  for  the  above  reaction  also 
yielded  calculated  N0X  concentrations  in  the  reference  troposphere 
that  were  significantly  larger  than  some  subsequent  NO^  observa- 
tions (Noxon,  1978)-*  As  indicated  in  a subsequent  section,  the 

The  NO2  observations  may  be  interpreted  as  N0X  (i.e.,  NO  + NO2) 
because  the  ratio  NO/NO2  is  zero  for  night  conditions  and 
approximately  unity  for  daylight  conditions. 
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effect  from  the  much  faster  rate  coefficient  for  HO,,  + NO 
OH  + NOp  on  the  perturbation  of  the  ozone  column  by  subsonic 
aircraft  is  to  increase  the  effects  of  the  CO  and  CHl4  mechanisms 
(M4a,  b)  as  well  as  to  decrease  the  effect  of  the  stratospheric 
N0x  catalytic  cycle.  Furthermore,  the  Howard-Evenson  (1977) 
rate  coefficient  yields,  now,  lower  background  NO  concentrations 
in  the  reference  troposphere  which  produce  larger  NO^  perturba- 
tions than  before  for  the  same  level  of  NO  emissions  by  sub- 

x 

sonic  aircraft.  The  probable  net  effect  of  the  faster  Howard- 
Evenson  rate  coefficient  for  subsonic  flight  alone  appears  then 
to  be  an  augmentation  of  the  previous  increase  in  the  ozone  column 
at  northern  middle  latitudes.  For  supersonic  flight  alone,  the 
weakening  of  the  stratospheric  N0x  catalytic  cycle  would  tend  to 
diminish  the  previous  decrease  in  the  ozone  column  at  the  north- 
ern middle  and  high  latitudes. 

The  results  for  the  combined  fleet  of  subsonic  and  super- 
sonic aircraft  in  the  Widhopf  model  were  based  on  a reaction 
rate  coefficient  for  the  HO^  + NO  OH  + NO,,  reaction  that  was 
nearly  identical  to  that  in  the  Crutzen  model.  Hence,  the  same 
qualitative  arguments  given  above  for  individual  fleets  indicate 
that  for  a combined  fleet  of  subsonic  and  supersonic  aircraft, 
the  faster  Howard-Evenson  reaction  rate  coefficient  may  also 
augment  the  previous  net  increase  in  the  ozone  column  at  the 
northern  middle  latitudes  for  the  same  combination  of  N0x 
emissions . 

The  foregoing  qualitative  considerations  suggest  a probable 
combination  of  a stronger  effect  from  the  tropospheric  mechanism 
of  ozone  production  with  a weaker  effect  from  the  stratospheric 
mechanism  of  ozone  destruction  than  before.  The  consequence  of 
this  result  is  to  make  the  2-D  modeling  of  the  reference  and 
perturbed  troposphere  nearly  as  important  as  that  of  the  stra- 
tosphere, in  spite  of  the  fact  that  the  contribution  of  unper- 
turbed tropospheric  ozone  to  the  total  ozone  column  is  only  of 
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the  order  of  10  percent  or  so.  For  this  reason,  it  becomes  of 
interest  to  review  the  available  2-D  results  for  the  effects 
of  aviation  on  ozone,  as  determined  from  the  Crutzen  and  Widhopf 
models,  which  are  based  on  1976  chemistry. 

2.1  RESULTS  AND  CHARACTERISTICS  OF  THE  1976  CRUTZEN  MODEL 

Because  of  the  rather  massive  amount  of  data  generated  by 
2-D  models  for  both  the  reference  and  perturbed  (by  aircraft 
effluents)  tropospheres  and  stratospheres,  it  becomes  necessary 
to  examine  such  data  by  emphasizing  situations  of  major  interest. 
These  latter  conditions  may  be  taken  to  be  the  northern  middle 
latitudes  during  late  summer  on  the  following  grounds:  (a)  most 
of  the  world's  airline  traffic  is  usually  assumed  to  be  concen- 
trated at  latitudes  in  the  range  15°N  & 0 £ 55°N,  with  peak  NC>x 
and  H^O  emissions  occurring  at  about  A5°N  (Hidalgo  and  Crutzen, 
1977;  Widhopf  et  al.,  1977),  (b)  the  consequent  changes  in  the 
reference  ozone  column  from  subsonic  flight  are  largest  in  the 
latitude  range  25°N  s 0 s 45°N,  changes  that  tend  to  be  of  the 
same  order  of  magnitude  during  any  season.  However,  the  changes 
in  the  ozone  column  from  a combined  fleet  of  subsonics  and  super- 
sonics  show  that  they  tend  to  be  larger  during  late  summer  than 
in  early  summer;  (c)  the  intensity  of  the  UV-B  solar  radiation 
is  largest  during  summer,  as  illustrated  by  Figs,  ^a  and  *Jb 
(Venkateswaran  et  al.,  1975)  for  winter  and  summer.  These 
figures  indicate  that  the  UV-B  radiation  at  a typical  wavelength 
at  the  northern  middle  latitudes  decreases  by  more  than  two  orders 
of  magnitude  during  winter  as  compared  with  summer.  These  re- 
sults are  also  typical  of  those  in  Western  and  Northern  Europe 
(Venkateswaran  et  al.,  1975),  and  (d)  outdoor  human  activity 
(with  the  consequent  exposure  to  UV-B  radiation)  increases  during 
the  longest  daylight  conditions  of  the  summer  season  at  middle 
latitudes . 
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FIGURE  4.  Isolines  of  the  daily  average  of  UV  solar  radiation  flux  at  a typical  wavelength  of 
300.4  nm  during  winter  (January)  and  summer  (July).  Numbers  on  left-hand  side  are 
in  microwatts/m2-nm,  whereas  those  on  the  right  are  in  photons/m^-sec-nm.  (Source: 
Venkateswaran  et  al , 1975) 


Figure  2 indicated  schematically  that  ozone  is  generated 
mainly  in  the  stratosphere  by  mechanism  1,  and  the  downward 
direction  of  the  arrow  (in  that  figure)  showed  that  zone 
descended  towards  the  earth's  surface.  The  ozone  concentration 
as  a function  of  altitude  is  then  the  net  result  of  the  chemical 
mechanisms  of  ozone  production  and  destruction,  as  indicated  in 
Fig.  2,  and  the  ozone  transports  by  the  air  motions.  Figure  5 
shows  quantitative  results  for  ozone  in  the  reference  atmosphere 
at  45°N  latitude  during  late  summer  (August  30)  as  was  simulated 
by  the  Crutzen  model.  This  figure  shows  both  the  ozone  concen- 
tration and  the  growth  of  the  ozone  column  with  decreasing  alti- 
tude. The  ozone  concentration  reaches  a peak  value  near  23  km 
(34  mbar  ),  when  the  ozone  column  has  reached  about  52  percent 
of  its  total  value.  The  negative  vertical  gradient  of  the 
stratospheric  ozone  concentration  above  about  25  km  is  controlled 
by  the  interplay  between  mechanism  Ml  and  M2  (Fig.  2)  as  well  as 
the  air  motions.  The  ozone  production  mechanism  (Ml)  becomes 
weaker  with  decreasing  altitude  due  to  the  decay  of  atomic  oxygen 
with  decreasing  altitude.  For  the  conditions  of  Fig.  5»  the 
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calculated  concentration  of  atomic  oxygen  drops  from  5*7  x 10 
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per  cm0  at  the  stratopause  to  4 x 10  per  cnv3  at  20  km.  Com- 
parisons of  the  characteristic  times  for  the  chemistry  and  dyna- 
mics indicate  that  the  former  are  short  relative  to  that  of  the 
dynamics  in  the  upper  and  middle  stratosphere,  while  the  opposite 
is  true  in  the  lower  stratosphere  (e.g.,  Johnston,  1975).  There- 
fore, the  chemistry  is  fast  and  in  near  equilibrium  in  the  upper 
and  middle  stratosphere  but  slow  in  the  lower  stratosphere,  where 
ozone  has  been  considered  in  the  past  to  be  nearly  chemically 
inert.  These  different  relative  roles  of  the  chemistry  and 
dynamics  as  a function  of  altitude  in  the  stratosphere  will  be 
used  in  subsequent  sections  for  considerations  of  either  the 
parameterizations  of  the  dynamics  in  the  lower  stratosphere  or 
the  consistency  of  the  chemistry  itself.  Hence,  the  positive 
vertical  gradient  of  the  stratospheric  ozone  concentration  below 
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FIGURE  5.  Simulated  ozone  concentration  (bottom  scale)  and  ozone  column  growth  with  decreasin 
altitude  (top  scale)  for  the  reference  atmosphere  at  45°N  during  summer  (August  30) 
(Source:  Hidalgo  and  Crutzen,  1977) 


about  23  km  in  Fig.  5 is  the  result  of  the  dominant  role  of  the 
ozone  transports  by  the  air  motions.  The  results  for  the  buildup 
of  the  ozone  column  with  decreasing  altitude  (Fig.  5)  shows  that 
(a)  about  3 percent  of  the  ozone  column  is  created  above  4o  km 
(2.7  mbar)  and  (b)  the  contribution  of  the  tropospheric  ozone 
(below  10  km)  is  about  9 percent  of  the  total  ozone  column, 
which  has  a value  of  328.6  m-atm-cm  (or  Dobson  units)  at  45°N 
latitude  during  late  summer.  The  former  result  for  the  upper 
stratosphere  becomes  of  interest  in  attempts  to  validate  the 
understanding  of  the  N0x  chemistry  by  taking  advantage  of  the 
locally  large  perturbations  of  ozone  during  the  infrequent  solar 
proton  events  (i.e.,  M7,  Table  1;  Heath  et  al.,  1977);  whereas 
the  latter  result,  for  the  contribution  of  tropospheric  ozone 
to  the  ozone  column,  is  of  interest  in  considerations  of  the 
effect  of  high-altitude  subsonic  flight  on  the  ozone  column. 

Table  2a  shows  the  vertical  distributions  of  the  daylight, 

zonal  averages  of  0,  and  NO  related  species  as  simulated  by 

the  Crutzen  model  with  1976  chemistry  at  45°N  latitude  during 

late  summer.  The  first  ozone  column  represents  its  concentra- 

■5 

tion  in  molecules  per  cm  . The  subsequent  columns  provide  the 
volume  mixing  ratio  (i.e.,  the  concentration  of  a given  chemical 
constituent  normalized  with  that  of  air)  of  ozone,  nitrous  oxide, 
N0x  (defined  here  as  NO  + N02  + 2N20<-  + N0^)>  nitrogen  pentoxide 
(N20j-),  NO^,  and  nitric  acid  (HNO^).  The  numbers  in  parentheses 
denote  exponents  for  the  given  value.  For  example,  the  peak 
ozone  concentrations  at  23.41  km  given  as  4.38(12)  denotes 
4.38  x 10  molecules/crrr . Table  2b  shows  similar  data  for 
water  vapor,  H0x,  CO,  and  CH^.  The  first  H20  column  is  the 
simulated  relative  humidity,  whereas  the  subsequent  columns  rep- 
resent the  volume  mixing  ratio  for  water  vapor,  hydrogen  peroxide 
(H202),  OH,  CO,  and  CH^. 
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TABLE  2a.  DAYLIGHT,  ZONAL  AVERAGES  OF  OZONE,  NoO,  NOx,  N2O5,  NO3  and  HNO3 
SPECIES  IN  THE  REFERENCE  ATMOSPHERE  AS  WAS  SIMULATED  BY  THE  2-D 
CRUTZEN  PHOTOCHEMICAL  MODEL  WITH  1976  CHEMISTRY  AT  45°N  DURING 
LATE  SUMMER  (AUGUST  30)* 


Z(km) 

[03] 

°3 

n2o 

N°x 

N2°5 

NO, 

1 

hno3 

54.39 

4.23(10) 

3.70(-6) 

1.51 (-8) 

1.00(-8) 

0 

0 

3 . 84  (-12) 

52.13 

2.88(10) 

1 . 94  ( - 6 ) 

1.51(-8) 

1 . 51 ( -8 ) 

0 

0 

3 . 84 (-12) 

49.85 

4.69(10) 

2 . 42 ( -6 ) 

1.81 ( -8) 

1.71(-8) 

0 

0 

8 . 00 ( - 12 ) 

47.61 

7.52(10) 

2 . 96( -6) 

2 . 32 ( -8) 

1 - 77 ( -8) 

0 

0 

1 . 73( -11 ) 

45.38 

1.34(11) 

3.92(-6) 

2. 99( -8) 

1.77 (-8) 

0 

0 

4 . 07 (-11) 

43.18 

2.44(11) 

5 . 29 ( -6 ) 

3. 84 (-8) 

1.71(-8) 

0 

0 

9 . 35 (-11) 

41.02 

4.39(11) 

7 . 00 ( -6 ) 

4 . 89 ( -8 ) 

1 . 62 ( -8 ) 

0 

0 

2 . 03( -10) 

38.89 

7.26(11) 

8. 50(-6) 

6. 12 ( -8 ) 

1 . 49 ( -8 ) 

0 

0 

3 . 76 ( - 10 ) 

36.81 

1.13(12) 

9 . 67 ( -6 ) 

7 . 44 { —8 ) 

1.36 (-8) 

1 . 28 ( - 10 ) 

3 . 35  ( - 

•12) 

4 . 90 (-10) 

34.78 

1.63(12) 

1 . 02 ( -5) 

8. 74( -8) 

1 . 24 ( -8) 

3 . 35 ( - 10 ) 

4 . 6 1 ( - 

12) 

7.08(-10) 

32.79 

2.17(12) 

9 . 98 ( -6 ) 

1 . 00( -7 ) 

1. 12 ( -8 ) 

5 . 65 ( - 10) 

5 . 18  ( - 

12) 

9 . 65 ( - 10 ) 

30.82 

2.73(12) 

9. 23 ( -6 ) 

1. 13 ( -7 ) 

9. 89 (-9) 

7 . 95 ( - 10 ) 

4.99(- 

•12) 

1 . 32 { —9) 

28.92 

3.29(12) 

8. 18 ( - 6 ) 

1 • 26 ( -7 ) 

8 . 45 ( -9 ) 

9.40(-10) 

4 . 23  ( - 

■12) 

1 . 79( -9) 

27.08 

3.81(12) 

6.96(-6) 

1 • 42 ( —7 ) 

6 . 79 ( -9 ) 

9 . 27 ( - 10) 

3. 17(- 

•12) 

2 . 34 ( -9 ) 

25.24 

4.19(12) 

5.65(-6) 

1 - 59 ( -7 ) 

5 . 03( -9 ) 

7 . 73( - 10 ) 

2.25(, 

■12) 

2 . 79 ( -9 ) 

23.41 

4.38(12) 

4 . 37  ( - 6 ) 

1- 77 (-7) 

3 .42 ( -9 ) 

5 . 54 ( - 10 ) 

1 . 45  ( - 

•12) 

2 . 93 ( -9 ) 

21.61 

4.33(12) 

3. 19 ( -6 ) 

1 - 96 ( -7 ) 

2. 1 3 ( -9 ) 

3 . 34 ( - 10 ) 

8. 13(  - 

•13) 

2.71 ( -9 ) 

19.81 

3.91(12) 

2. 15 ( -6 ) 

2. 1 3 ( -7 ) 

1 . 20 ( -9 ) 

1 . 73 ( - 10 ) 

4 . 25 ( - 

•13) 

2 . 14 ( -9 ) 

18.03 

3.15(12) 

1 . 28 ( -6 ) 

2 . 27 ( -7 ) 

6 . 08 ( - 10 ) 

6.761*11) 

1 . 68  ( - 

•13) 

1 . 42 ( -9 ) 

16.25 

2.13(12) 

6.41(-7) 

2 . 37 ( - 7 ) 

2.86(-10) 

1 . 99 (-11) 

5 . 02  ( - 

•14) 

7 . 49 ( -10) 

14.47 

1.20(12) 

2 . 74 ( -7 ) 

2 . 42 ( -7 ) 

1 . 43  ( - 10 ) 

6 . 40 ( - 12 ) 

1 . 67  ( - 

•14) 

3.43( -10) 

12.66 

7.17(11) 

1 • 25 ( -7 ) 

2 . 46 ( - 7 ) 

8 . 43( - 11) 

0 

0 

2 . 04 ( - 10 ) 

10.82 

5.18(11) 

6.99-8) 

2.48-7 

6.54-11 

0 

0 

1 . 65 ( - 10 ) 

8.90 

5.47(11) 

5 . 66 ( -8 ) 

2.49(-7) 

6 . 36 (-11) 

0 

0 

1 . 59 ( - 10 ) 

6.88 

6.02(11) 

4.97(-8) 

2 . 49 ( -7 ) 

7 . 12( - 11 ) 

0 

0 

1 .47 ( -10 ) 

4.75 

7.09(11) 

4.63(-8) 

2.50(-7) 

1 . 49 ( - 10 ) 

0 

0 

3. 15  ( -10 ) 

2.52 

8.60(11) 

4 . 50 ( -8 ) 

2. 50(-7 ) 

3.90(-10) 

0 

0 

9.67(-10) 

1.94 

9.07(11) 

4.47(-8) 

2.50(-7) 

5.281-10) 

0 

0 

1. 13( -9) 

1.36 

9.34(11) 

4. 33(-8) 

2 . 50 ( -7 ) 

8 . 78 ( - 10 ) 

0 

0 

1 . 20( -9 ) 

0.78 

9.26(11) 

4.02(-8) 

2.50(-7) 

1.53(-9) 

0 

0 

1 . 09 ( -9 ) 

0.21 

8.67(11) 

3. 53(-8) 

2.50(-7) 

2 . 82 ( -9 ) 

0 

0 

5.11(-10) 

The  first  ozone  column  is  the  concentration  (molecules/cnr) ; the 
subsequent  columns  denote  the  volume  mixing  ratio.  Numbers  in 
parentheses  denote  powers  of  10;  e.g.,  the  concentration  of  ozone 
at  54.39  km  is  4.23  x 10l0  molecules/cm^. 


30 


rt 


TABLE  2b.  DAYLIGHT,  ZONAL  AVERAGES  OF  H2O,  H2O2,  OH,  CO  AND  CH4  SPECIES 
IN  THE  REFERENCE  ATMOSPHERE  AS  WAS  SIMULATED  BY  THE  2-D  CRUTZEN 
PHOTOCHEMICAL  MODEL  WITH  1976  CHEMISTRY  AT  45°N  DURING  LATE 
SUMMER  (AUGUST  30)* 


7.2C-7) 
5.0(-7) 
4.4(-7) 
4.3(-7) 
7.0(-7) 
1.11-6) 
2 . 0 ( -6 ) 
5.0(-6) 

Hill 

6.0(-5) 
1.3-4) 
3 . 0 ( -4 ) 
7.5(-4) 
1 . 6 ( -3) 
3 . 4 ( -3 ) 

8.7- 3) 
1 - 9 ( - 2 ) 
4. 7(-2) 
i.l(-l) 

Hi:!! 

3.5(-l) 

5.4(-l) 

5.1- 1) 
5.8(-l) 

4.8- 1 
4 . 9 1 - 1 ) 

5.2- 1) 
6.0(-l) 
7.3(-l) 


h2o 

H2°2 

01 

CO 

4 . 44  ( 

-6 

8. 15( 

-11) 

1.00( 

4.44( 

-6 

5.37 

-10) 

5.52 

-10) 

7 36  ( 

4.46 

-6 

6.05 

-10) 

4.99 

-10 

6.30( 

4.46 

-6 

7.02 

-10) 

4.47 

-10 

k 5.81 ( 

4.45 

-6 

8.73 

-10) 

4.05 

-10 

' 5 . 48 ( 

4.39 

-6 

1.03 

-9) 

3.37 

-10 

5. 1 9 ( 

4.32 

-6 

1.20 

-9) 

2.52 

-10 

4.85( 

4.23-6) 


4. 12(-6) 
4. 10(-6) 
4 . 12 ( -6) 
4. 16(— 6) 


4 . 45 ( -6 ) 


1 . 28( -5 ) 
1 . 88 ( -5 ) 
3.54(-5) 
9 . 48 (- 5 ) 
2 . 05( -4 ) 
7.32-4) 
2.231-3) 
6.06-3) 
7. 111-3) 
8.59(-3) 
1.04 (-2) 
1.34(-2) 


1 . 45( -9) 

1 . 86( -9 ) 

2 . 50 ( -9 ) 

3 . 06 ( -9 ) 
3.34-9) 
3.19(-9) 
2.70(-9) 
2.09(-9) 

1 . 55 ( -9) 
1.14-9) 

8 . 69 ( - 10 ) 

6.34(-10) 

4.37 (-10) 

3.77-10) 

4 . 24 ( - 10) 

6.59-10 

7 . 34 ( - 10) 

1.44-9) 

6.27-10) 

2.45(-9 

2.32(-9) 

1.41-9) 

5 . 12 ( - 10 ) 
5. 12(-10) 


1.65(-1Q) 


2 . 97 (-11) 
1.61  -ll) 
9.03(-12 
5. 23(-12) 
3.14-12) 
1.97 (-12) 
1.30(-12 
9 . 04 ( - 13) 
6. 11(-13) 
3. 67 ( - 13) 
2 . 14 (-13) 
1 . 35 ( - 13 ) 
1.15(-13) 


8.60(-14) 
1.17-13) 
7 . 60 ( - 14 ) 
1.48-13 
1 . 45 (-13) 
1.31-13) 
1 . 07 ( - 13 ) 
1 . 0( — 13) 


4 . 46 ( -8 ) 
4. 06 (-8) 


3.63(-8) 
3.24(-8) 
2.89(-8) 
2 . 60 (-8 ) 


2.38(-8) 
2 . 22 ( -8 ) 
2 . 17 ( -8 ) 
2.28(-8) 
2. 64 (-8) 
3 . 35  ( -8 ) 
4.34(-8) 
5 . 52  ( -8 ) 
7.24(-8) 
9. 07 (-8) 
9.  S8(-3) 
1 . 09 ( -7 ) 
1 • 24 ( -7 ) 
1.46-7) 
1.56(-7) 
1.69-7) 
1.85 (-7) 
2. 02 (-7) 


2 . 31 ( -7 ) 
2 . 31 (-7 ) 
2 . 53( -7 ) 
2. 88 (-7) 
3. 34 (-7) 
3 . 90 { -7 ) 
4.55-7) 
5 . 26 ( -7 ) 
5 . 96 ( -7 ) 
6 . 58 ( -7 ) 
7 . 15 ( —7 ) 
7 . 70( -7 ) 
8.27(-7) 
8 . 92 ( -7 ) 
9.68-7) 
1 .05(— 6) 
1 . 13 ( —6 ) 
1 . 22 ( -6 ) 
1.29 (-6) 
1 . 34 (-6) 
1 . 38 ( -6 ) 
1 .40(-6) 
1 . 43 ( -6 ) 
1.44(-6) 
1.44-6) 
1 . 45 ( -6 ) 
1 .47 ( -6 ) 
1 .47 ( —6 ) 
1.48(-6) 
1.49 (-6) 
1.50(-6) 


The  first  H?0  column  denotes  the  relative  humidity,  and  the  subsequent 
ones  the  volume  mixing  ratio. 
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Figure  6 is  a plot  of  the  data  for  the  volume  mixing  ratios 
in  Tables  2a,  b.  The  abscissa  of  the  figure  extends  over  13 
orders  of  magnitude  as  split  at  the  bottom  for  the  solid  lines 
(for  N0x,  HNO^,  OH,  and  H2O2)  and  at  the  top  for  the  dashed 
lines  (N20,  CO,  CH^,  0^,  H20).  This  figure  shows  that  the  ozone 
mixing  ratio  has  a broad  peak  near  35  km  (5.5  mbar),  which  is  of 
interest  in  considerations  of  the  heating  of  stratospheric  air 
by  the  absorption  of  UV  radiation  by  ozone.  The  NgO  mixing  ratio 
is  nearly  constant  in  the  troposphere  and  decreases  with  increas-  ! 

ing  altitude  in  the  stratosphere  as  a consequence  of  the  mech- 
anisms M5  and  M6  (Table  1).  The  N0x  mixing  ratio  in  the  tropo- 
sphere decreases  with  increasing  altitude  from  the  high  values  of 
the  NC>x  emissions  at  the  earth's  surface  to  a broad  minimum  in 
the  upper  troposphere.  The  N0x  mixing  ratio  increases  in  the 
stratosphere  as  a consequence  of  mechanisms  that  start  with  M6 
(Table  1).  The  OH  profile  is  of  interest  because,  through  the 
reaction  OH  + N02(+M)  -*•  HN0^(+M),  it  is  (a)  a temporary  reservoir 
for  NO  that  weakens  the  stratospheric  NO  cycles  (M2,  M3c), 
and  (b)  a permanent  tropospheric  sink  due  to  the  solubility  of 
HNO^  in  water  (M9).  The  tropospheric  OH  profile  is  also  impor- 
tant because  OH  is  a sink  for  CO  and  CH^  (M4a,  b).  Furthermore, 

OH  scavenges  not  only  CO  and  CH^,  but  also  N0n  and  S0o  (M9)  as 
well  as  CH  Cl  F (M17).  Finally,  Fig.  6 shows  that  the  simu- 

X y Z 

lated  volume  mixing  ratio  of  Ho0  drops  rapidly  in  the  troposphere 

d -2 

and  lower  stratosphere  from  the  surface  value  of  1.3  x 10  to 
about  4 x 10  ^ in  the  middle  stratosphere. 

Figure  7 shows  results  of  high-altitude  flight  on  the  refer- 
ence ozone  column  as  obtained  from  the  Crutzen  model  with  1976 
chemistry  (Hidalgo  and  Crutzen,  1977).  This  figure  shows  the 
computed  change  in  the  reference  ozone  column  as  a function  of 
latitude  for  assumed  individual  fleets  of  subsonic  and  SST  air- 
craft after  6 years  of  fleet  operations.  The  former  is  assumed 
to  be  a large  fleet  operating  in  the  year  1990  at  10.8  km,  and 
the  latter  a small  fleet  operating  at  18  km.  The  total  NO 
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FIGURE  6.  Simulated  vertical  profiles  of  daylight,  zonal  average  of  the  volume  mixing  ratio  of 
indicated  species  in  the  reference  atmosphere  at  45°N  during  summer  (August  30)  by 
the  Crutzen  model  with  1976  chemistry.  (Source:  Hidalgo  and  Crutzen,  1977) 


4-?7--8-9 


FIGURE  7.  Relative  effects  of  individual  subsonic  (top)  or  supersonic 
(bottom)  flight  on  the  atmospheric  ozone  column  during  summer 
(August  30)  after  six  years  of  fleet  operations.  NO-  emissions 
were  for  hypothetical  fleets  assumed  to  operate  in  the  1990 
time  frame.  (Source:  Hidalgo  and  Crutzen,  1977) 
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emissions  (as  NO^)  were  2.06  x 10^  kg  per  year  for  the  subsonlcs 
and  about  10  percent  of  this  value  for  the  SSTs . The  dashed 
lines  in  the  figure  indicate  the  assumed  latitudinal  distribu- 

Q 

tions  of  the  N0x  emissions  in  units  of  10  gm  of  NO^  per  year. 

It  is  thus  seen  that  the  bulk  of  the  NOx  emissions  were  assumed 
to  take  place  at  the  northern  middle  latitudes.  This  figure 
shows  that  subsonic  aircraft  produce  near  145°N  latitude  a 
maximum  ozone  increase  of  0.76  percent,  as  the  result  of  the 
increase  in  the  reference  ozone  column  (Fig.  5)  from  328.6  to 
the  perturbed  value  of  331.1  Dobson  units.  The  corresponding 
NOx  perturbation  at  10.8  km  was  about  300  percent  (Hidalgo  and 
Crutzen,  1977 ) • 

The  numerical  data  from  the  Crutzen  model  for  the  reference 
atmosphere  was  obtained  after  15  years  of  Integrations.  The 
basic  structure  of  the  model  was  based  on  the  conservation  of 
chemical  species,  dynamics  that  were  empirically  constructed  so 
as  to  match  the  H„0  and  0-.  observations  in  the  lower  stratosphere, 
and  rate  constants  for  the  chemical  mechanisms  in  Table  1 as 
given  in  Table  3.  The  model  was  completed  by  specifying  boundary 
conditions  at  the  ground  and  near  the  tropopause  (approximately 
km)  as  well  as  the  heterogenous  removal  rates  for  the  tro- 
pospheric rainout  sinks  of  NC>X,  HNO^,  H^O^,  and  CH^0oH.  Figure 
8 shows  a comparison  of  the  simulated  and  observed  0^  mass  mix- 
ing ratio,  a result  that  is  constrained  (as  that  of  H^O)  in  the 
lower  stratosphere  by  the  empirical  adjustments  of  the  2-D  eddy 
coefficients.  Further  considerations  of  the  2-D  dynamics  will 
be  given  in  a subsequent  section. 

2.2  RESULTS  AND  CHARACTERISTICS  OF  THE  1976  WIDHOPF  MODEL 

The  basic  formulation  of  the  Widhopf  2-D  model  was  the  same 
as  in  the  Crutzen  2-D  model,  except  for  the  empirical  representa- 
tion of  the  eddy  motions.  Instead  of  using  the  Ho0  and  0^  dis- 
tributions, as  in  the  Crutzen  model,  the  Widhopf  model  was  based 
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TABLE  3.  1976  REACTION  RATE  COEFFICIENTS  FOR  AN  0XYGEN-HY0R06EN-NITR06EN  ATMOSPHERE* 

SOURCE:  HIDALGO  AND  CRUTZEN,  1977 


SI 
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AIR  PRESSURE,  mbai 


CRUTZEN  MODEL,  BASED  ON  1976 
CHEMISTRY,  JUNE  30,  35°  N 
KRUEGER  ET  AL  . ROCKET  DATA, 
JUNE  18,  1970,  34°  N 


MIXING  RATIO,  Mgm/gm 


FIGURE  8.  Comparison  of  simulated  and  observed  ozone  mixing  ratio  at 
35°N  during  June.  The  simulated  values  in  the  lower  and 
middle  stratosphere  were  constrained  by  their  use  in  the 
parameterization  of  the  eddy  transports.  The  observations 
are  given  by  the  data  of  Krueger  et  al.  (1973).  (Source: 
Hidalgo  and  Crutzen,  1977) 


ALTITUDE,  km 


on  (a)  extrapolations  to  the  stratosphere,  up  to  latitudes  of 
50°N,  by  Luther  (1973)  of  the  tropospheric  transports  of  heat 
(Oort  and  Rasmusson,  1971),  and  (b)  matching  of  parts  of  the 
time  history  in  the  dispersion  by  the  stratosphere  of  excess 
C-14  that  resulted  from  the  1961-1962  nuclear  weapon  detona- 
tions in  the  atmosphere  at  approximately  70°N.  This  2-D 
parameterization  of  the  dynamics  was  then  tested  by  reproduction 
of  segments  of  the  time  history  of  the  dispersion  of  other 
atmospheric  radioactive  debris;  i.e.,  W-185  and  Zr-95.  Further- 
more, the  ozone  distributions  in  the  reference  atmosphere  were 
calculated  utilizing  1976  chemistry;  these  results  indicated  good 
agreement  with  observations  as  a function  of  both  altitude, 
latitude,  and  season  (Widhopf  et  al.,  1977).  It  must  then  be 
noted  that  this  parameterization  of  the  eddy  transports  in  the 
reference  stratosphere  is  Independent  of  the  ozone  chemistry. 

The  chemistry  in  the  two  models  was,  in  general,  equivalent; 
especially  since  the  adoption  of  the  tropospheric  ozone  produc- 
tion mechanisms  in  the  Widhopf  model  was  based  on  earlier  experi- 
ence with  the  Crutzen  model.  However,  it  is  of  interest  to 
indicate  some  differences  in  the  important  reaction  rates  in 
the  stratosphere  as  used  in  these  two  models. 


Five  reaction  rates  in  the  stratospheric  HOx~NOx  cycles 
have  been  identified  as  being  the  most  important  ones  of  those 
in  Table  3 for  changes  of  the  ozone  column  by  SSTs . These  re- 
sults are  based  on  the  use  of  a 1-D  model  and  on  considerations 
of  the  uncertainty  in  the  magnitude  of  the  calculated  changes 
in  the  ozone  column  by  SSTs  due  to  the  uncertainty  in  the  mag- 
nitude of  the  reaction  rate  coefficients  (Duewer  et  al.,  1977). 
With  the  notation  in  Table  3,  these  reactions  are: 


R(21) 
R(l4) 
R(13) 
R(  7 ) 
R(15) 


OH  + H02 
H02  + NO 


h2o  + o2 

OH  + N0n 


03  + H02  -*■  OH  + 2 02 
O^  + OH  -*■  H02  + 02 


OH  + N02(+M) 


hno3(+m) 


A comparison  of  the  reaction  rate  coefficients  for  the  fore- 
going reactions,  as  used  in  the  Crutzen  and  Widhopf  2-D  models, 
is  as  follows: 


Reaction 

R(21) 
R(l4) 
R(  13) 
R(  7) 

R(15) 


Crutzen  Model 

5.0  x 10-11 
2.2  x 10-13 

1.0  x 10-13exp ( -1250/T) 
1.6  x 10_12exp(-1000/T) 
4.2  x 10~22exp(-170/T) (M) 
4(10)-11(M)  + 1 . 6 ( 10 ) -2T 


Widhopf  Model 


2.0  x 10 
2.3  x 10' 


2.76  x 10~X 3exp ( 880/T ) 
1.166  x 10 1 xp ( 220/T ) + [ M ’ 


Besides  the  different  formulation  of  the  eddy  motions  and 
difference  in  reaction  rate  coefficients  in  the  Widhopf  and 
Crutzen  models,  other  significant  sources  of  differences  between 
these  models  are:  (a)  the  formulation  of  the  Ho0  distributions 
in  the  troposphere  and  stratosphere,  which  were  either  specified 
(Widhopf)  or  parameterized  (Crutzen),  (b)  the  magnitude  of  the 
removal  rates  for  the  tropospheric  rainout  sinks,  (c)  the  use 
of  the  height  (Widhopf)  instead  of  the  pressure  (Crutzen)  as 
the  vertical  coordinate,  the  latter  being  in  better  accord  with 
airline  flight  practice,  and  (d)  the  numerical  integrations  of 
the  continuity  equations  for  the  conservation  of  the  chemical 
species.  While  the  Widhopf  model  solves  for  each  continuity 
equation,  the  Crutzen  model  utilizes  a family  concept  that  must 
invoke  chemical  equilibrium  among  short-lived  species.  Some  of 
the  possible  limitations  introduced  by  the  latter  approximations 
are  described  by  Chang  (1977). 

Table  4 and  Fig.  9 show  results  for  the  reference  atmosphere, 
similar  to  those  described  above  for  the  Crutzen  model,  that  were 
obtained  from  the  Widhopf  model.  However,  in  contrast  with  the 
daylight  averages  of  the  Crutzen  results  (Table  2 and  Fig.  6), 
the  corresponding  Widhopf  values  (Table  4 and  Fig.  9)  are  for  a 
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TABLE  4.  24-HOUR,  ZONAL  AVERAGE  OF  SPECIES  IN  THE  REFERENCE  ATMOSPHERE  AS  SIMULATED 
BY  WIDHOPF,  et  al . (1977)  AT  40°N  DURING  LATE  SUMMER  (AUGUST  15)* 


VOLUME  mixing  ratio 


Simulated  vertical  profiles  of  24-hour,  zonal  average  of  the  volume  mixing  ratio  of 
indicated  species  in  the  reference  atmosphere  at  40°N  during  summer  (August  15)  by 
the  Widhopf  model  with  1976  chemistry.  (Source:  Widhopf,  1978) 


24-hour  average;  a difference  that  becomes  important  for  species 
with  a short  life  relative  to  the  diurnal  cycle  (i.e.,  NO,  N02, 

, NO^,  H202,  and  OH).  Therefore,  a comparison  of  the  results 
in  Pigs.  6 and  9 must  be  limited  to  the  species  with  a relatively 
longer  life.  However,  even  then,  the  interpretation  of  any 
difference  in  such  results  must  take  into  account  the  difference 
in  the  chemistry  as  illustrated  previously  by  the  reactions  R15 
and  R21.  These  results  indicate  the  following: 

• The  1976  Widhopf  HNO^  profiles  are  somewhat  higher  than 
the  Crutzen  values  above  about  15  km  altitude.  For 
example,  at  24  km  altitude,  the  Widhopf  HNO^  is  higher 
than  the  Crutzen  value  by  about  11  percent. 

• In  spite  of  the  higher  Widhopf  HNO^  concentrations, 
the  chemistry  used  in  this  model  tends  to  give  slightly 
higher  ozone  concentrations  than  the  Crutzen  values 
above  35  km  altitude. 

• The  Widhopf  ozone  concentrations  in  the  troposphere  are 
significantly  lower  than  the  Crutzen  values  by  a factor 
as  low  as  approximately  0.15  at  10  km  altitude.  As  indi- 
cated subsequently,  the  2-D  Crutzen  ozone  concentration 
for  these  conditions  is  lower  than  observations  by  about 
40  percent  (Fishman  and  Crutzen,  1978),  which  means  that 
the  Widhopf  tropospheric  ozone  values  are  also  signifi- 
cantly lower  than  observations. 

• The  Widhopf  CO  concentrations  in  the  stratosphere  are 
lower  than  the  Crutzen  values  by  a factor  of  approximately 
0.3  at  30  km  altitude. 

A further  comparison  of  the  Widhopf  and  Crutzen  results  is 
given  in  a subsequent  section,  which  includes  the  Fishman-Crutzen 
2-D  results  for  the  troposphere  based  on  the  Howard-Evenson  rate 
coefficient  and  the  three-body  component  for  the  rate  coefficient 
of  the  reaction  CO  + OH  (+M)  -*•  C02  + H + (M). 


Figure  10  shows  results  for  a combined  fleet  of  subsonics 

and  SSTs  as  were  obtained  from  the  Widhopf  model  utilizing  1976 

chemistry  (Widhopf  et  al.,  1977).  The  'total  NO  emissions 

q x 

used  in  these  results  were  2.81  x 10^  kg  of  N02  per  year, 
which  was  distributed  in  altitude  and  latitude.  The  bulk 
of  the  N0x  emissions  (97  percent)  was  assumed  to  take  place  at 
altitudes  below  15  km,  so  that  the  SST  component  of  the  N0X 
emissions  in  the  altitude  range  15  s z s 19  km  was  a relatively 
small  fraction  of  the  total  N0x  emissions.  The  latitudinal  dis 
tribution  of  the  N0x  emissions  (Fig.  7)  was  similar  to  that  of 
the  subsonics  for  the  Crutzen  model. 
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FIGURE  10.  Effects  of  a combined  fleet  of  subsonic  and  supersonic  air- 
craft on  the  atmospheric  ozone  column  during  summer  (June  15) 
and  fall  (October  15)  after  five  years  of  fleet  operations. 
N0X  emissions  were  for  projected  fleets  assumed  to  operate  in 
the  1990-time  frame.  (Source:  Widhopf,  et  al . , 1977) 


It  should  be  noted  that,  based  on  1-D  experience  (e.g., 
Oliver  et  al. , 1977;  Broderick,  1977),  the  above  slower  rate  for 
R(21)  in  the  Widhopf  (as  compared  with  that  of  the  Crutzep)  model 
would  tend  to  decrease  the  corresponding  perturbations  of  the 
reference  ozone  column  by  SSTs. 

Even  without  consideration  of  the  subsequent  1977  revision 
of  the  1976  chemistry  by  the  Howard-Evenson  measurement  of  the 
reaction  rate  for  R14,  important  caveats  in  regard  to  the  results 
in  Figs.  7 and  10  are  as  follows:  (1)  they  omitted  the  simula- 
tion of  chlorine  chemistry  in  the  atmosphere;  i.e.,  of  the  inter- 
actions among  the  HO  - NO  - CIO  cycles  in  an  oxygen-hydrogen- 

XX  X 

nitrogen-chlorine  atmosphere;  (2)  they  could  not,  at  the  time, 
incorporate  the  simultaneous  effect  of  the  aircraft  H^O  efflu- 
ents, which  could  again  play  a significant  role  in  the  coupling 
of  the  stratospheric  HO  - NO  - CIO  cycles;  and  (3)  they  had 

XX  X 

to  omit  consideration  of  the  peroxynitric  acid  (HO^NO^)  molecule 
or  reaction  HO2  + NC^  ■+  HO2NO2.  Although  it  was  considered 
likely  that  this  reaction  takes  place  in  the  atmosphere,  vir- 
tually nothing  was  known  of  its  kinetics  (Hidalgo  and  Crutzen, 
1977).  Furthermore,  it  should  be  noted  that  no  sensitivity 
analysis  is  yet  available  concerning  the  effect  of  subsonics  on 
the  calculated  change  of  the  ozone  column  from  the  uncertain 
knowledge  of  the  reaction  rate  coefficients,  and  even  some  of 
the  chemical  mechanisms,  in  the  troposphere. 
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3.  IMPACT  OF  CHANGES  IN  OZONE  CHEMISTRY  ON  2-D  MODELING 


Table  3 identified  the  1976  photochemistry  and  chemical 
kinetics  of  the  troposphere  and  stratosphere  that  has  been  used 
in  2-D  models  for  an  oxygen-hydrogen-nitrogen  (O-H-N)  atmosphere. 
This  formulation  must  not  only  be  extended  to  include  chlorine 
chemistry  (i.e.,  M2,  M17,  Table  1)  and  molecules  such  as  H02  NC>2, 
but  the  complexity  of  some  of  the  chemical  cycles  in  this  table 
is  such  as  to  produce  significant  uncertainties  on  their  mech- 
anisms and/or  rates.  These  complexities  have  prevented  direct 
measurements  of  some  of  these  reaction  rate  coefficients,  which 
had  then  to  be  estimated  theoretically  from  appropriate  ratios 
involving  other  reactions  with  better  known  rate  coefficients. 

As  indicated  in  the  previous  section,  it  has  been  possible  to 
identify  the  most  important  reactions  in  Table  3 (i.e.,  R7,  13, 
14,  15,  21)  for  the  effects  of  SSTs  on  ozone  for  an  O-H-N 
atmosphere.  Ongoing  research  has  recently  achieved  direct  mea- 
surements of  most  of  these  important  rate  coefficients.  Never- 
theless, because  of  the  remaining  uncertainty  in  the  knowledge 
of  the  tropospheric  and  stratospheric  chemistry,  as  well  as  the 
empirical  2-D  formulation  of  the  essentially  unknown  strato- 
spheric dynamics,  it  is  important  to  indicate  that  experience 
with  photochemical  models  has  brought  out  their  rather  severe 
inherent  limitations  due  to  the  availability  of  a significant 
number  of  degrees  of  freedom  for  the  tuning  of  a limited  set 
of  observations.*  Using  a reversed  chronological  order,  a brief 
summary  of  this  experience  is  as  follows: 


Because  of  the  four  types  of  basic  limitations  indicated 
previously  (Section  1.0),  it  is  important  to  emphasize  that 
most  of  these  inherent  limitations  would  also  be  applicable 
to  3-D  and  1-D  models. 


PHSCSOLNO  PAG* 


Tho  tvi’i'tit  de'vo 1 opine  nt  of’  a laser-magnetic  resonance 
detection  te'chnlque  ha:;  allowed  direct  measurements 
of  HO,  in  the  react  Ions  HO,  + NO  ► OH  + NO-,  ( H 1 h ) and 
0,  » HO.,  » Oil  + 0,  (H13).  This  technique  lias  yielded 

the  new  (1977,  1978)  rate  coefficients  for  these  reac- 
t Ions  as  described  bot’ore.  The  1977  result  for  H 1 h is 
a factor  of  about  *10  faster  than  the  197b  values  used 
In  the  C rut  ■/.on  (■’•■'  x 10  * ')  and  Widhopf  (0.3  x 10 
C-P  models.  The  1977  direct  measurement  of  the  rate 
coefficient  for  H 1 9 significantly  upset  previous  re- 
sult:;, which  had  Implied  an  accurate,  model-independent 
matching  of,  for  example,  either  vertical  ozone  pro- 
file's at  various  latitudes  and  seasons  (Widhopf,  Glatt, 
and  Kramer,  1977)  or  those  of  methane  and  nitric  acid 
(Hidalgo  and  drut/.en,  1°77).  It  then  becomes  evident 
that  artificial  comb  I nut  ions  of  react  ion  rate  coeffi- 
cient:; have'  yielded  an  adequate'  matching  be' tween  the' 
calculated  ami  e'bsorveei  di  st  r 1 but  l ems  of  chemical 
species  In  tlu'  re? fere'ne’.e  atmosphere. 

Initial  e'xpe'r  I me  til  a I results  (burrows  e't  al.,  1977) 
and  upper  limit  cents  l derat  ions  (Chang  ami  Kaufman, 

1 977 ; Howard,  1977)  for  the'  rate'  coef  f 1 clent  of 
OH  + lie?,  » 11,0  + 0 , (R.?l),  have  led  to  a recommended 
value  of  3.0  x 10  ^ cm  V molecule-sec  for  this  coef- 

ficient (NAHA,  1977).  This  recommended  value  is  nearly 
an  orde'r  of  magnitude  slower  than  the  197J|  value 
(.'.0  x 10“ !°)  used  In  the'  Climatic  Impact  Assessment 
I'rejgram  (CIAP).  The  twlslon  e'f  this  rate'  cetef f 1 c lent 
haei  a significant  Impact  e'n  the  prevleuts  CIAP  re'sults 
for  the  probable  e'ffrcts  e'f  hlgh-altltude'  flight  on 
the'  ozone'  column  (dreibe'ckor  e't  al.,  197*i;  National 
Acade'my  of  He  1 e'ttce's , 1979).  Mase'd  eui  1-D  model  results, 
the’  revision  of  this  rate  coefficient  produced  a sig- 
nificant eh'wngradl ng  of  the'  CIAP  assessment  of  the 


potential  effects  of  SSTs  on  the  reference  ozone  col- 
umn (Broderick,  1977;  Oliver  et  al.,  1977). 

• Inclusion  of  the  reaction  OH  + NOn(+M)  -*•  HNO^C+M)  (R15), 
in  the  CIAP  assessment  of  the  effects  of  SST  on  the 
ozone  column  brought  about  a downgrading,  by  a factor 

of  about  2,  of  a previous  1-D  assessment  of  such  ef- 
fects; which  had  been  made  without  the  use  of  this 
reaction  (Chang,  Hidalgo,  and  Johnston,  1975). 

• A pioneering  effort  on  the  2-D  modeling  of  atmospheric 
ozone  found  it  even  feasible  to  match  the  ozone  column 
as  a function  of  latitude  and  season  through  the  use 
of  only  the  0X  Chapman  cycle  (R1  to  R5 , Table  3)  and 
suitable  parameterization  of  the  eddy  motions 
(Frabhakara,  1963).  The  emphasis  of  this  effort  was 
to  bring  out  the  role  of  the  dynamics  in  the  strato- 
spheric distribution  of  ozone. 

The  foregoing  experience  indicates  that,  although  signi- 
ficant progress  has  been  made  during  recent  years  in  measuring 
important  reaction  rates  for  an  O-H-N  atmosphere,  the  quanti- 
tative 2-D  assessment  of  high-altitude  flight  on  the  ozone  col- 
umn has  yet  to  meet  satisfactorily  the  rigorous  standards  of 
science,  a result  that  is  a consequence  of  the  basic  constraints 
indicated  in  Section  1.  This  experience  also  indicates  that  con- 
sensus of  expert  human  Judgment  unsupported  by  critical  experi- 
ments has  not  been  able  to  resolve  the  very  complex  geophysical 
puzzles  Involved  in  the  consideration  of  the  ozone  (and  climate) 
change  problems.  Had  such  an  approach  been  used,  for  example 
in  the  time  of  Copernicus,  a most  probable  consensus  would  have 
been  that  the  earth  did  not  move. 

3.1  UNDERDESTRUCTION  OF  OZONE  WITH  1977  CHEMISTRY 

The  original  0x  cycle  (R1-R5,  Chapman  1930)  was  proposed 
to  explain  the  characteristics  of  the  ozone  column.  However, 

H9 


the  subsequent  gathering  of  more  accurate  ozone  data  developed 
a need  for  ozone  destruction  mechanisms  (i.e.,  M2,  Table  1)  so 
as  to  reduce  the  calculated  Chapman  ozone  column  to  the  observed 
values  (Johnston,  1975).  The  evolution  of  these  ozone  sinks 
started  with  the  introduction  of  the  H0x  catalytic  cycle  (i.e., 

X = OH  in  M2)  and  has  proceeded  with  the  subsequent  advents  of 
the  N0x  and  C10x  catalytic  cycles  in  the  stratosphere. 

Prior  to  the  direct  measurement  by  Howard  and  Evenson  of 
the  rate  coefficient  for  the  reaction  H0o  + NO  ■*  OH  + NO^,  the 
ozone  sinks  had  been  adjusted  in  the  Crutzen  2-D  model  with  the 
H0x  and  N0x  cycles  in  Table  3-  Likewise,  the  Widhopf  2-D  model 
had  used  these  same  cycles,  but  with  some  modification  in  the 
rate  coefficients.  However,  preliminary  results  using  the 
Howard  and  Evenson  rate  coefficient  in  the  Widhopf  model  for 
an  O-H-N  atmosphere  yield  ozone  columns  that  are  signi f leant ly 
larger  than  the  observed  values  at  all  latitudes  in  the  Northern 
Hemisphere.  Figure  11  shows  some  results  from  preliminary  in- 
vestigations using  the  Widhopf  model  that  gave  the  results  in 
Fig.  10,  except  for  the  following  modifications  (Widhopf  and 
Glatt,  1978):  (a)  change  to  the  faster  Howard-Evenson  rate 

coefficient,  (b)  changes  of  other  reaction  rates  to  the  1977 
NASA  recommended  values,  (c)  incorporation  of  the  N^O^  and  NO^ 
species,  (d)  modifications  of  the  radiation  calculations  to 
include  multiple  scattering  and  the  long  wavelength  absorption 
by  Ho0o,  (e)  change  from  the  24-hour  diurnal  to  the  daylight- 
night  average  of  the  concentration  of  chemical  species, 

(f)  extrapolations  of  the  2-D  vertical  eddy  coefficients  in  the 
middle  and  upper  stratosphere  based  on  the  Hunten  (1975)  instead 
of  the  previous  McElroy  et  al.  (1974)  1-D  vertical  eddy  coeffi- 
cients, the  former  being  about  a factor  of  two  smaller  than  the 
latter,  and  (g)  calculation,  instead  of  specification,  of  the 
water  vapor  distribution. 
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FIGURE  11.  Illustration  of  the  impact  of  the  Howard-Evenson  reaction 
rate  coefficient  on  the  simulated  ozone  mixing  ratio  for 
the  reference  atmosphere.  Data  in  Table  5 shows  an  in- 
crease of  about  28%  on  the  ozone  concentration  at  35  km 
altitude  from  use  of  the  Howard-Evenson  reaction  rate  co- 
efficient for  HO2  + NO  -*■  OH  + N02  instead  of  that  used  for 
the  results  shown  in  the  previous  figure.  The  use  of  the 
1978  Howard  rate  for  HO2  + O3  -*■  OH  + 202  does  not  correct 
the  ozone  imbalance  above  about  30  km  altitude.  Note  the 
use  of  a linear,  instead  of  a logarithmic  scale  so  as  to 
display  the  mismatch  in  the  ozone  profile  above  30  km. 
(Source:  Widhopf,  1978) 
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Figure  11  shows  that  the  Howard-Evenson  rate  coefficient 
yields  significantly  higher  ozone  concentrations  than,  for  example, 
the  Krueger  et  al.  (1973)  rocket  data  in  the  middle  and  upper 
stratosphere;  where  the  role  of  the  chemistry  relative  to  that 
of  the  dynamics  is  dominant  (Section  2.1  and  Fig.  5).  Table  5 
shows  the  modification  of  the  1976  profiles  at  40°N  during  the 
late  summer  (August)  in  Table  4 by  the  preliminary  calculations 
using  the  Widhopf  model  with  1977  chemistry.  The  results  in 
Table  5 are  given  by  the  ratio  of  the  values  based  on  1977  chem- 
istry to  those  with  1976  chemistry  in  Table  4.  Such  ratio  is 
applicable  to  species  with  a long  life  relative  to  the  diurnal 
cycle,  since  it  would  need  a correction  for  the  species  (e.g., 

OH,  etc.)  with  a short  life  relative  to  this  cycle  because 

of  the  difference  in  their  averaging.  However,  the  results  in 
Tables  2a,  4,  and  5 for  the  latter  species  can  provide  relative 
comparisons  for  the  daylight  averages  of  the  Widhopf  values  based 
on  1977  chemistry  and  the  Crutzen  results  based  on  1976  chemistry. 
These  results  indicate  the  following: 

• At  the  35-km  altitude  for  the  maximum  increase  in  the 
ozone  mixing  ratio  in  Fig.  11,  Table  5 shows  that  the 
increase  in  ozone  concentration  by  the  Howard-Evenson 
rate  coefficient  in  the  Widhopf  model  is  about  28  per- 
cent. The  tropospheric  ozone  has  been  increased  by 


factors  of  better  than  2. 


The  corresponding  increase 

18 


in  the  total  ozone  columns  is  16.5  percent;  i.e., 

from  8.1  x lO1^  molecules/cm2  (Table  4)  to  9.5  x 10 

2 

molecule/cm  (Tables  4 and  5).* 


At  the  other  latitudes  and  time  of  the  year,  this  increase 
in  the  ozone  column  can  be  20  to  30  percent  higher  than 
either  the  observed  or  calculated  value  with  the  1976 
chemistry  (Widhopf  and  Glatt,  1978). 
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TABLE  5.  RATIO  OF  VALUES  BASED  ON  HOWARD-EVENSON  RATE  (DAILY  AVERAGES)  TO  THOSE  CORRESPONDING 
IN  TABLE  4 (24-HOUR  AVERAGE)  AT  40°N  DURING  LATE  SUMMER. 

(Source:  Widhopf,  1978) 
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• Table  5 for  the  Widhopf  model  shows  that  the  HNO^ 
profiles  based  on  1977  chemistry  have  decreased  sig- 
nificantly in  the  stratosphere  and  upper  troposphere, 
whereas  those  of  H,,0  in  the  middle  and  upper  strato- 
sphere have  increased  appreciably.  The  reason  for  the 
former  effect  must  take  into  account  the  corresponding 
changes  in  the  chemistry,  including  the  consideration 
of  NO^  species  as  described  previously,  whereas  the 
latter  effect  is  due  to  the  calculation,  instead  of  the 

1976  specification,  of  the  water  vapor  distribution. 

• The  results  for  N0x  in  Table  5>  together  with  those 
in  Tables  2a  and  4 indicate  that  the  (daylight)  NC>x 
concentrations  in  the  upper  troposphere  have  decreased 
by  more  than  an  order  of  magnitude  from  the  1976 
Crutzen  value  (Table  2a)  to  the  Widhopf  value  based  on 

1977  chemistry  (Table  5).  This  Widhopf  result,  which 
is  of  interest  for  subsonic  flights,  is  a trend  toward 
the  Noxon  observed  values. 

The  underdestruction  of  ozone  from  the  use  of  the  Howard- 
Evenson  rate  coefficient  with  1977  chemistry,  as  indicated  in 
Fig.  11,  becomes  important  through  the  following  effects 
(Crutzen  and  Howard,  1978):  (a)  a faster  rate  for  H0„  + NO  + 

OH  + N02  (R14)  produces  a larger  production  of  OH  radicals 
below  about  35  km  than  before.  This  implies  a more  efficient 
conversion  of  N0x  into  the  HNO^  sink  through  the  reaction 
OH  + N02(+M)  -*■  HNO^+M)  (R15),  and  (b)  the  weakening  of  the 
ozone  destruction  by  the  reaction  OH  + 0^  -*•  H02  + 0_,  (R7)  in 
the  HO  catalytic  cycle  M3b  (Table  1).  This  latter  result  is 
a consequence  of  the  last  three  reactions  (Rl4,  R35,  R2)  in  M4b, 
which  produce  as  a net  the  reaction  H02  + 02  ■+  OH  + 0^  (R13); 
a result  that,  in  turn,  is  opposite  to  that  of  R7  in  M3b. 

Relative  results  from  a 1-D  model  by  S.C.  Liu,  utilizing 
the  1977  Howard-Evenson  rate  coefficient  for  H0_,  + NO  ■+  OH  + NO,, 
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indicated  a decrease  by  about  9 percent  in  the  ozone  column 
when  using  the  1978  instead  of  the  1977  NASA  rate  coefficient 
for  HO^  + 0^  -*■  OH  + 2 0^  (Zahniser  and  Howard,  1978).  Similar 
preliminary  considerations  using  the  Widhopf  2-D  model  also 
indicate  the  same  effects,  i.e.,  the  1978  rate  coefficient  for 
HC>2  + 0^  -*■  OH  + 2 C>2  reduces  the  imbalance  in  the  ozone  column 
to  about  10  percent  above  the  observed  value  for  the  conditions 
shown  in  Fig.  11.  However,  this  effect  on  the  ozone  column  is 
caused  by  an  ozone  decrease  in  the  lower  stratosphere,  which  is 
a region  where  the  mass  transports  become  important  (Section  2.1, 
Fig.  5)  and  the  ozone  concentrations  are  relatively  high  (Table 
4).  Hence,  the  calculated  imbalance  in  the  ozone  mixing  ratio 
in  the  middle  and  upper  stratosphere  caused  by  the  1977  rate 
coefficient  for  H0^  + NO  -*•  OH  + NO^  is  almost  unaffected  by 
the  more  recent  measurement  of  the  rate  coefficient  for 
HC>2  + 0^  OH  + 2 C>2  (Fig.  11).  These  preliminary  results 
suggest  a need  for  a reassessment  of  the  remaining  reaction 
rate  coefficients  for  an  0-H-N  atmosphere  as  well  as  considera- 
tions of  the  chlorine  chemistry  for  the  middle  and  upper  strato- 
sphere . 

3.2  ESTIMATES  OF  TROPOSPHERIC  N0x  AND  H202  RAINOUT  SINKS 

Section  2 presented  qualitative  considerations  which 
indicated  that  the  results  in  Fig.  7,  for  the  effect  of  sub- 
sonic flight  at  10.8  km  on  the  ozone  column,  will  be  increased 
as  a consequence  of  the  faster  Howard-Evenson  rate  coefficient; 
which  yields  also  a lower  calculated  tropospheric  NC>x  back- 
ground. To  provide  quantitative  considerations  of  this  effect, 
it  is  of  interest  to  describe  now  some  relevant  results  based 
on  1977  chemistry  as  obtained  by  a global,  tropospheric,  non- 
dynamical,  chemical  2-D  model  (Fishman  and  Crutzen,  1978). 

From  the  point  of  view  of  aircraft  effects  on  the  ozone  column, 
emphasis  here  is  on  the  Northern  Hemisphere  as  well  as  on 
(a)  the  indirect  assessments  of  the  amount  of  NC>x  background  in 
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the  troposphere,  (b)  a qualitative  assessment  of  the  direction 
of  the  change  of  the  ozone  column  by  subsonic  flight  in  the  up- 
per troposphere  from  the  revision  of  the  1976  chemistry  by  both 
the  Howard-Evenson  and  Noxon  N0x  data,  and  (c)  assessments  of 
the  magnitude  of  the  heterogenous  removal  rates  of  chemical 
species  in  the  tropospheric  sinks. 

The  basic  approach  used  in  the  Pishman-Crutzen  (1978)  2-D 

investigations  consists  in  examining  the  compatibility  of  assumed 

NO  backgrounds  in  the  troposphere  with  the  corresponding  0, 

and  CO  budgets  in  the  Northern  and  Southern  Hemispheres.  Figure 

2 indicated  that  these  budgets  are  essentially  controlled  by 

mechanisms  M3(a,  b,  c)  and  M^(a,  b),  while  Table  1 showed  that 

these  mechanisms  involve  the  HO  (i.e.,  OH,  H0~),  NO  (NO  and 

NO-,),  and  CH^  species.  From  Table  3>  the  considerations  of 

the  0^  and  CO  budgets  put  emphasis  on  the  odd  hydrogen  HO 
j x 

species  (R6  to  R23),  N0x  (R33>  35),  and  the  methane  oxidation 
reactions  (R24  to  R32).  Two  important  considerations  in  these 
2-D  calculations  in  the  troposphere  are  the  short  lifetimes  of 
the  H0x  species  and  their  consequent  equilibrium  with  a pre- 
scribed background  of  0^,  CO,  H~,0,  CH^  and  the  unknown  N0x. 

These  considerations  allow  the  calculation  of  the  H0x  species  by 
using  a quasi-steady  state  or  equilibrium  2-D  model  with  pre- 
scribed tropospheric  distributions  of  0^,  CO,  H_,0,  CH^,  and  NO,, 
as  a function  of  season  in  both  hemispheres.  Other  important 
consequences  of  these  considerations  are  the  lack  of  require- 
; ments  for  both  time-dependent  calculations  for  a given  season 

and  parameterizations  of  the  tropospheric  2-D  dynamics. 

Because  of  the  current  lack  of  adequate  statistical  data 
for  NO  and  HO  species,  the  basic  aim  in  these  calculations 

X X 

is  then  to  examine  the  sensitivity  of  the  tropospheric  0^  and 
CO  budgets  as  a function  of  prescribed  levels  of  N0x<  The  most 
probable  values  of  N0x  (and  HO,.)  in  the  troposphere  would  then 
be  given  by  the  assumed  N0x  amount  that  is  most  consistent  with 
the  balancing  of  the  0^  and  CO  budgets.  The  emphasis  on  con- 
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sistency  instead  of  exact  matching  of  these  budgets  stems  from 
the  rather  significant  uncertainties  in  the  current  knowledge 
of  these  budgets.  It  is  important  to  emphasize  that  these  2-D 
calculations  utilized  1977  chemistry  as  given  by: 

• The  1977  Howard-Evenson  rate  coefficient  for  R14 

• The  three-body  component  for  the  CO  and  OH  reaction, 
which  becomes  important  in  the  troposphere.  Hence,  in 
addition  to  the  two-body  reaction  R8  in  Table  3»  these 
tropospheric  calculations  used  the  reaction  CO  + OH(+M) 
C02  + H(+M)  with  a three-body  rate  coefficient  given 

by  7-3  x 10~^[M]  cm  Vmolecule-sec  (Chan  et  al.,  1977). 

• The  photolysis  rate  (J)  for  0^  + hv  ■+  0^  + O^D)  (R3a ) 
that  includes  a temperature  dependence  in  its  quantum 
yield  <j>(A).  This  latter  parameter  is  defined  by 


Qa<J>d\  , 


where  Q is  the  solar  flux,  cr  the  ozone  cross  section, 
and  A the  wavelength  in  the  range  2*90  £ A s 320  nm. 

The  temperature  dependence  of  the  quantum  yield  can 
alter  the  calculated  production  rate  of  O^D)  by  as 
much  as  20  percent  over  the  range  of  tropospheric 
temperatures . 

• The  NASA  (1977)  revisions  of  the  other  rate  coeffi- 
cients in  Table  3. 

• The  heterogenous  removal  of  H202  and  CH^O^H. 

In  addition,  these  2-D  calculations  included  the  reaction 
of  methyl  chloroform  (CH^CCl^)  with  OH;  i.e.,  CH^CCl^  + OH  -*• 

H20  + CH2  CCl^.  Other  modifications  of  the  reactions  in  Table  3 
for  this  tropospheric  model  were  as  follows:  the  dropping  of 
Rl,  R3b  and  R4  in  the  oxygen  reactions;  the  dropping  of  R6 , RIO, 
R12,  and  R20  for  the  hydrogen  species;  and  the  replacement  of 


R26  with  CH302H  + OH 


ch3o2  + h2o. 
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Figure  12  defines  the  assumed  input  levels  of  NC>x  used  in 
these  2-D  calculations.  The  profile  A is  representative  of  un- 
polluted conditions  and  it  is  consistent  with  both  the  Noxon 
(1977)  and  the  surface  observations  of  Drummond  (1977)  for  clean 

air  in  Wyoming.  This  profile  gives  a tropospheric  column 

14  2 

density  of  2.1  x 10  molecules/cm  . The  profile  B is  assumed 

to  represent  polluted  conditions  with  a corresponding  column 

1^  2 

density  of  4.6  x 10  molecules/cm  . Profile  A is  used  at  lati- 
tudes 0 £ 15°N  and  in  the  entire  Southern  Hemisphere,  whereas 
profile  B is  used  at  0 2 45°N.  Intermediate  profiles  are  used 
in  the  latitude  range  15°N  & 0 s 45°N.  These  2-D  calculations 
utilized  four  assumed  levels  of  N0x  which  were  given  by: 

• A low,  constant  vertical  profile  of  the  volume  mixing 

-12 

ratio  at  every  latitude  given  by  a value  of  10  or 

14 

a tropospheric  column  density  of  0.2  x 10  molecules/ 
cm2 

• The  profiles  A and  B or  corresponding  column  densities 
of  2.1  x 10^  and  4.6  x 10^  molecules/cm2  with  lati- 
tudinal distributions  as  described  before.  The  average 
NO  column  density  in  the  Northern  Hemisphere  is  then 

X T jl  p 

3.4  x 10  molecules/cm 

• Twice  those  of  profiles  A and  B 

• A high,  constant  vertical  profile  of  the  volume  mixing 

ratio  at  every  latitude  given  by  a value  of  10~10  or 

14 

a tropospheric  column  density  of  22.9  x 10  molecules/ 
cm2. 

Hence,  the  assumed  NOx  distributions  in  the  troposphere  of  the 
two  hemispheres  varied  by  about  two  orders  of  magnitude. 

3.2.1  Assessment  of  Tropospheric  N0x  in  the  Northern  Troposphere 

Table  6 shows  the  sensitivity  of  the  tropospheric  ozone 
production  (M4a,  b.  Table  1)  and  destruction  (M3a,  b,  c)  mech- 
anisms in  the  Northern  Hemisphere  to  the  foregoing  assumed 
amounts  of  N0X>  The  ozone  production  by  mechanisms  4a  and  4b 
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FIGURE  12.  Assumed  tropospheric  N0X  profiles  for  sensitivity  analysis  of  the  tropospheric  ozone 
and  CO  budgets  to  NQX  and  HO.  (Source:  Fishman  and  Crutzen,  1978) 


increases  with  increasing  NOx,  whereas  the  ozone  destruction  is 

somewhat  insensitive  to  the  amount  of  N0x  due  to  the  dominant 

role  of  mechanism  3a  in  the  destruction  of  ©zone.  Hence,  the 

ozone  production  mechanisms  dominate  the  net  production  of  ozone 

at  the  high  values  of  the  assumed  NO^  amounts.  However,  the 

ozone  flux  coming  through  the  tropopause  in  the  Northern  Hemi- 

11  2 

sphere  is  about  0.7  x 10  molecules/cm  -sec  (Danielsen  and 

Mohnen,  1977),  whereas  the  ozone  flux  destroyed  at  the  ground 

in  the  Northern  Hemisphere  is  also  about  0.7  x 1011  molecules/ 

2 

cm  -sec  (Fabian  and  Pruchniewicz , 1977).  Since  the  net  ozone 
production  in  the  troposphere  is  apparently  negligibly  small, 
the  results  in  Table  6 indicate  that  the  two  high  values  of  the 
assumed  NOx  amounts  are  not  probable  for  the  troposphere.  Fur- 
thermore, since  the  N0x  amount  given  by  profiles  A and  B are 
consistent  with  the  column  observations  of  Noxon  (1978)  and  the 
surface  NO  observations  of  Drummond  (1977)  for  clean  air  in 
Wyoming,  the  results  in  Table  6 indicate  that  the  NOx  given  by 
these  profiles  are  the  most  likely  ones  in  the  troposphere. 

The  bottom  row  in  Table  6 indicates  the  annual,  diurnal 
average  OH  concentrations  corresponding  to  the  assumed  NOx 
amounts.  The  OH  concentrations  play  a role  in  the  CO  budget 
through  mechanisms  M4a  and  M4b . Table  7 shows  the  sources 
and  sinks  for  the  CO  budget  for  the  Northern  Hemisphere  as 
given  by  Seiler  (1974)  for  other  than  those  from  mechanisms  4b 
and  4a,  respectively.  The  sources  account  for  6.6  x 10^  kg/yr 
for  the  contributions  to  the  CO  budget  by  anthropogenic,  oceans, 
forest  burning,  and  oxidation  of  the  hydrocarbons;  whereas 
the  sinks  account  for  3.9  x 1011  kg/yr  for  the  soil  uptake  and 
the  transport  of  CO  to  the  stratosphere.  The  table  includes 
the  calculated  2-D  values  for  the  source  and  sink  from  mechan- 
isms 4b  and  4a,  respectively,  as  a function  of  the  assumed  N0x 
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TABLE  6.  SENSITIVITY  OF  OZONE  BUDGET  TO  PRESCRIBED  NOx  CONCENTRATIONS 
FOR  THE  NORTHERN  HEMISPHERE  (Source:  Fishman  and  Crutzen,  1978) 

14  2 

Assumed  NOx  (TO  molecules/cm  ) 


Ozone  Production  (M4a,  b) 

0.2 

3,4 

6.8 

22.9 

(10  molecules/cm  sec) 

Ozone  Destruction  (M3a,  b,  c) 

0.1 

1.8 

3.3 

10.0 

(10  molecules/cm  sec) 

Net  Ozone  Production 

1.2 

1.2 

1.2 

1.4 

11  2 
(10  molecules/cm  sec) 

Annual  Diurnal  Average 

OH  Concentration  in  Northern 

-1.1 

0.6 

2.1 

8.6 

Hemisphere  (10  molecules/cm  ) 

2.9 

3.3 

3.9 

9.5 

TABLE  7.  SENSITIVITY  OF  CO 

BUDGET  TO  PRESCRIBED  N0y 

CONCENTRATIONS 

FOR  THE  NORTHERN  HEMISPHERE*  (Sources:  Seil 

er,  1974 

9 

Fishman 

and  Crutzen, 

1978) 

Assumed  NO*  (10^ 

2 

molecules/cm 

Sources  (10^  kg/year) 

0.2 

3.4 

6.8 

22.9 

Anthropogenic 

5.4 

5.4 

5.4 

5.4 

Oceans 

0.4 

0.4 

0.4 

0.4 

Forest  Burning 

0.4 

0.4 

0.4 

0.4 

Oxidation  of  Hydrocarbons 

0.4 

0.4 

0.4 

0.4 

CH4  oxidation  (M4b) 

1.7 

1.7 

2.2 

3.7 

Total  Sources 

8.3 

8.3 

8.8 

10.3 

Sinks  (lO11  kg/year) 

Soil  uptake 

3.0 

3.0 

3.0 

3.0 

Transport  to  Stratosphere 

0.9 

0.9 

0.9 

0.9 

CO  Oxidation  (M4a) 

6.7 

8.9 

11.2 

31.8 

Total  Sinks 

10.6 

12.8 

15.1 

35.7 

Imbalance  (10"  kg/year) 

-2.3 

-4.5 

-6.3 

-25.4 

For  the  Southern  Hemisphere  and  N0X  amount  of  3.4  x 10^  molecules/cm^, 
for  example,  the  corresponding  source  and  sink  total  values  are  3.6  and 
7.7  (in  units  of  10''  kg/yr),  respectively  (Fishman  and  Crutzen,  1978). 
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amounts.  The  results  in  this  table  indicate  that  the  imbalance 
between  the  total  sources  and  sinks  is  aggravated  significantly 
as  the  assumed  N0x  amount  increases,  an  effect  that  is  due  to 
the  resulting  increase  in  the  average  OH  concentration  (bottom 
row.  Table  6)  and  overdestruction  of  CO.  Therefore,  as  in  the 
case  of  the  ozone  budget,  the  two  higher  values  of  the  assumed 
N0x  amounts  are  again  not  probable  for  the  troposphere.  In 
regard  to  the  new,  significantly  lower  estimate  of  the  annual, 
diurnal  average  of  the  OH  concentrations  in  the  northern 
troposphere  (~  3*3  x 10  molecules/cm  , Table  6),  it  should 
be  noted  that  it  implies  more  reasonable  magnitudes  for  mech- 
anisms 4b  and  4a  than  before.  As  indicated  by  Seiler  (1974), 
previous  estimates  for  the  source  M4b  were  in  the  excessive  range 
from  15  x 10^  to  40  x 10"*''1'  kg/yr,  and  for  the  sink  M4a  from 
19-4  x 10^  to  50  x 10^  kg/yr.  The  significantly  lower  values 
for  the  CHjj  and  CO  oxidations  in  Table  7 stem  from  the  new 
revisions  in  the  chemistry.  Nevertheless,  the  still  significant 
imbalance  of  the  CO  budget  for  the  lower  NOx  amounts  in  Table  7 
suggests  that  (a)  each  of  the  Seiler  sources  and  sinks  may  have 
to  be  reevaluated,  and  (b)  the  tropospheric  OH  must  be  much 
better  known  as  a function  of  latitude  and  season.  It  is  impor- 
tant to  note  that  because  of  the  high  variability  of  OH,  such 
data  should  be  accompanied  by  simultaneous  measurements  of 
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3.2.2  Ozone  Change  Due  to  Subsonics  with  1977  Chemistry 

The  results  in  Fig.  7 for  the  effects  of  subsonic  flight 

at  10.8  km  at  45°N  latitude  during  summer  were  based  on  an  NO 

-11  / x 

background  given  by  a volume  mixing  ratio  of  6.5  x 10  (Table 
i),  a result  that  is  a factor  of  about  16  larger  than  that  of 
• • polluted  profile  B in  Fig.  12  (~  4 x 10-lt~).  Therefore, 

: J.-r  to  examine  qualitatively  the  direction  of  the  ozone 
* wn  In  Fig.  7 as  a consequence  of  the  Howard-Evenson 
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and  Noxon  data,  it  becomes  of  interest  to  return  to  the  results 
for  the  ozone  production  mechanism  in  Table  6.  As  indicated 
by  mechanisms  4a  and  4b  in  Table  1,  the  oxidations  of  CO  and 
CH^  produce  one  0^  molecule  for  each  CO  and  four  0^  for  each 
CHjj  molecule.  However,  a limiting  step  in  the  CH^  oxidation 
is  the  photolysis  of  CH^O.  Consideration  of  the  photolysis 
rates  in  the  photochemical  reactions  R30a,  b,  reduces  the  effi- 
ciency of  the  methane  oxidation  to  the  production  of  only  2.5 
molecules  of  0^  for  each  molecule  of  CH^.  Thus,  an  upper  limit 
for  the  production  of  0^  is  given  by  the  sum  of  the  CO  oxida- 
tion and  2.5  times  that  of  CH^.  However,  from  the  actual  model 
results,  it  becomes  possible  to  determine  a corresponding  true 
"O^  efficiency  production  factor"  for  each  of  the  four  assumed 
N0x  distributions.  Table  8 gives  a summary  of  these  results, 
which  include  the  magnitude  of  the  CO  and  CH^  oxidations,  the 
upper  limit  0^  production,  the  actual  ozone  production  that 
was  given  in  Table  6,  and  the  0^  efficiency  production  factor. 

Thus,  for  the  NO  profiles  A and  B (second  column),  the  upper 

x 11  2 
limit  0^  production  of  3-5  x 10  molecules/cm  -sec  is  given 

by  the  CO  oxidation  and  2.5  times  that  of  methane.  This  upper 

limit  0.,  production,  together  with  the  actual  ozone  production 
^11  2 

of  1.8  x 10  molecules/cm  -sec  yield,  then,  an  efficiency  factor 
of  0.5.  Important  results  in  Table  8 are  (a)  the  dominant 
role  of  the  CO  relative  to  the  CH^  oxidation  and  (b)  the  increase 
in  both  the  upper-limit  0^  production  and  0^  efficiency  produc- 
tion factor  with  increasing  N0X  amounts. 
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TABLE  8.  SENSITIVITY  OF  TROPOSPHERIC  OZONE  PRODUCTION  TO  PRESCRIBED 
NOx  CONCENTRATIONS  FOR  THE  NORTHERN  HEMISPHERE 

(Source:  Fishman  and  Crutzen,  1978) 

Assumed  NOx  (1014  molecules/cm2) 


0.2 

3.4 

6.8 

22.9 

CO  oxidation  (M4a) 

11  ? 

(10  molecules/cm  sec) 

1.8 

2.4 

3.0 

8.5 

CH.  oxidation  (M4b) 

^ 11  0 
(10  molecules/cm  sec) 

0.4 

~ 0.4 

~ 0.6 

~ 1.0 

Upper  Limit  0,  Production  (M4a,  b) 

11  ^ 2 
(10  molecules/cm  sec) 

2.8 

~ 3.5 

~ 4.4 

~ 10.9 

Ozone  Production  (M4a,  b) 

(1011  molecules/cm2sec) 

0.1 

1.8 

3.3 

10.0 

0q  Efficiency  Production  Factor 

J (M4a,  b) 

0.03 

0.5 

0.75 

0.92 

The  relative  importance  of  the  Howard-Evenson  and  Noxon 
data  in  the  ozone  production  by  NOx  effluents  in  the  upper 
troposphere  may  then  be  assessed  by  examining  the  reactions  in 
mechanisms  M4a;  i.e.. 


( R 8 ) 

CO  + OH  (+M) 

-¥ 

CO  + H (+M) 

(Rll ) 

H + 02  + M 

-+ 

HOp  + M 

(R14 ) 

H0o  + NO 

-> 

OH  + NOp 

(R35) 

N02  + hv 

-> 

NO  + O 

(R2) 

O + 0o  + M 

c 

-* 

O^  + M 

CO  + 2 Op 

c°2  + o3 

where  (R8)  includes  the  two-body  reaction  with  (+M)  = 0.  The 
production  of  ozone  takes  place  from  the  reaction  of  molecular 
oxygen  with  atomic  oxygen  (R2),  the  latter  being  produced  by 
photolysis  of  NOp  (R35).  The  role  of  N0x  and  the  Howard-Evenson 
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reaction  rate  takes  place  through  the  formation  of  NO^  from 
HC>2  (Rl4);  which  is  produced  in  the  above  chain  from  every  H 
in  (R8),  through  (Rll).  The  relative  fraction  of  the  HO^  mole- 
cules that  become  Involved  in  Rl4  may  be  determined  by  examining 
the  competing  reactions  for  H0?  in  Table  3;  i.e.,  R13,  14,  17, 

21 , and  27 : 

(R13)  C>3  + H02  -*•  OH  + 2 02 

(R14)  H02  + NO  -*■  OH  + N02 

(R17)  H02  + H02  - H202  + 02 

(R21 ) OH  + H02  H20  + 02 

(R27)  CH302  + H02  - CH302H  + 02 

This  relative  fraction  is  then  given  from  the  normalized  values 

of  k^Cx^j/Ek^tx^] > where  k^  is  the  rate  constant,  [x^l  the  con- 
centration of  the  reactants  for  H02,  and  i the  reactions  R13, 

14,  17,  21,  and  27. 

Table  9 shows  results  for  these  relative  fractions  for  the 
2-D  model  conditions  of  the  subsonic  flight  in  Fig.  7>  and  those 
given  by  the  Howard-Evenson  and  representative  N0x  profile  B 
(Fig.  12).  The  former  Includes  conditions  for  both  the  refer- 
ence and  perturbed  tropospheres  as  determined  from  the  Crutzen 
tropospheric-stratospheric  2-D  model  with  1976  chemistry.  The 
latter  are  those  from  the  Fishman-Crutzen  2-D  model  calculations 
for  the  reference  troposphere  including  the  Howard-Evenson  rate 
coefficient,  the  modified  R8  reaction  to  include  the  three-body 
component,  the  revisions  in  the  quantum  yield  for  the  photolysis 
of  0(*D),  and  the  modified  reaction  rates  as  given  by  NASA  (1977). 
The  [NO]  concentration  in  Table  9 for  the  Fishman-Crutzen  2-D 
calculations  corresponds  to  that  of  profile  B (Fig.  12)  at 
10.8  km  and  45°N  latitude.  The  results  in  Table  9 for  the 
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TABLE  9.  COMPETITION  OF  REACTIONS  INVOLVING  THE  HO?  RADICAL  FOR  IMPORTANT  SUBSONIC  FLIGHT 

CONDITIONS  (~  10.8  km,  45&N,  August  30) 
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[CH^C^]  concentrations  are  based  on  the  approximation  [CH^O^]  - 
0.5  [HO^],  results  that  are  further  constrained  by  the  signi- 
ficant uncertainty  (unmeasured)  in  the  reaction  rate  for  (R27) 
involving  H02  and  CH^O,,.  The  results  for  the  reference  and 
perturbed  atmospheres  with  1976  chemistry  indicate  the  shift 
in  the  competition  from  Rl4  for  HC>2  due  to  the  increased  N0X 
by  subsonic  aircraft.  These  shifts  indicate  (a)  an  enhancement 
in  the  competition  of  (Rl4)  for  HO^  by  a factor  of  3*3  as  a 
result  of  the  N0x  effluents,  which  produced  the  local  ozone  in- 
crease of  12  percent,  and  (b)  a suppression  of  the  R17  competi- 
tion of  H02  for  itself  by  a factor  of  1.8  as  a consequence  of 

the  NO  effluents.  The  results  in  Table  9 for  the  reference 
x 

tropospheres  indicate  an  enhancement  in  the  competition  from 
Rl4  for  H02  by  a factor  of  1.8  as  a result  of  the  revision  of 
the  1976  to  the  1977  chemistry.  Hence,  the  increased  role  of 
Rl4  in  the  competition  for  H02  in  the  reference  troposphere 
within  the  dominant  M4a  mechanism,  coupled  with  the  weakening 
of  the  stratospheric  NO  catalytic  cycle  by  Rl4,  indicates  that 
the  change  of  the  ozone  column  given  in  Fig.  7 for  the  1976 
chemistry  should  be  enhanced  by  the  1977  chemistry.  The  direc- 
tion of  this  effect  is  somewhat  corroborated  by  preliminary  re- 
sults from  the  Widhopf  model  using  the  Howard-Evenson  (1977) 
rate  coefficient,  which  indicate  that  the  previous  peak  enhance- 
ment of  the  ozone  column  in  Fig.  10  increases  from  a value  of 
about  1.4  percent  to  about  3*5  percent,  even  though  the  SST 
component  of  the  combined  fleet  was  increased  by  a factor  of 
about  3 in  the  later  calculations;  which,  however,  included  the 
effect  of  simultaneous  injections  of  N0x  and  H^O  engine  effluents 
(Widhopf  and  Glatt,  1978). 

3.2.3  Heterogeneous  Removal  Rates  of  in  Tropospheric  Sinks 

A fundamental  unknown  in  tropospheric  chemistry  is  the 
removal  rate  of  water-soluble  atmospheric  gases  by  processes 
other  than  gas-phase  reactions  (e.g.,  M9,  Table  1,  and  Fig.  2). 


Estimates  of  these  removal  rates  have  been  obtained  by  Fishman 
and  Crutzen  (1978)  from  the  sensitivity  of  the  imbalance  of  the 
CO  budget  (Table  7)  to  OH  through  the  heterogeneous  removal  rates 
of  H2O2.  These  2-D  calculations  used  the  reactions 

H ^0 2 heterogeneous  removal 
CH^C^H  -*■  heterogeneous  removal, 

both  of  which  were  assumed  to  have  the  same  removal  rate  given 
by 

k ( z )^et  = k°het  exp  C — 0 • ^6 ( z— 5 ) H , z > 5 km 

where  k°het(sec-1)  is  a constant  for  the  heterogeneous  removal 

rate  below  5 km  altitude.  The  exponential  decay  at  altitudes 

above  5 km  reproduces  the  observations  of  Davidson  et  al.  (1966) 

for  removal  of  radioactive  debris.  Table  10  gives  results  of 

these  sensitivity  analyses  for  a specified  NO  amount  with  an 

14  x 2 

average  column  density  of  3-4  x 10  molecules/cm  (Table  7); 
which  is  given  by  profiles  A or  B (Fig.  12),  depending  on  the 
latitude,  as  indicated  previously.  This  table  shows  three 
values  for  k°het,  the  largest  one  being  the  rate  used  in  all 
the  previous  results  (Tables  6 to  8).  This  table  includes  cor- 
responding values  of  ^©2  in  the  lower  (1  km)  and  upper  (10  km) 
troposphere  as  well  as  the  corresponding  annual,  diurnal  average 
of  OH  in  the  Northern  Hemisphere. 

TABLE  10.  SENSITIVITY  OF  OH  TO  REMOVAL  RATES  OF  H2O2  FOR  THE 
ASSUMED  N0X  PROFILES  A AND  B IN  THE  NORTHERN  HEMISPHERE 


(Source: 

Fishman  and 

Crutzen,  1978) 

k° 

K het 

(k°  r1 

het' 

H202  (ppb)  , 

H202  (ppb)  , 

(OH)  ave 

sec”1 

1 km 

10  km 

5 3 

10  molecules/cm 

2.0  x 10”5 

13.9  hr 

0.34 

0.23 

3.3 

5.0  x 10'6 

2.3  days 

1.40 

0.62 

4.8 

2.4  x 10”6 

4.8  days 

2.46 

0.78 

5.7 
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Prom  the  results  In  Tables  7 and  10,  it  can  be  seen  that 
the  two  slower  heterogeneous  rates  In  Table  10  would  Increase 
the  average  OH  concentration  and  aggravate  significantly  the 
imbalance  of  the  CO  budget.  Hence,  an  estimate  for  the  hetero- 
geneous removal  rate  of  H202  below  5 km  is  given  by  k°het  ^ 

2 x 10”^  per  second  or  (k°  s 13*9  hour. 

3.3  COMPARISON  OF  2-D  MODEL  TROPOSPHERIC  RESULTS 

Section  2.0  presented  some  comparisons  between  the  Crutzen 
and  Widhopf  2-D  models.  These  considerations  are  now  extended 
to  include  the  following:  (a)  N0x  concentrations  in  the  upper 
troposphere,  (b)  heterogenous  removal  rates  of  water  soluble 
species,  and  (c)  the  parameterization  of  the  water  vapor  in  the 
reference  stratosphere. 

Figure  13  shows  a comparison  between  (a)  the  tropospheric 
N0X  profiles  at  middle  northern  latitudes  during  summer  cal- 
culated by  the  Crutzen  and  Widhopf  2-D  models  with  1976  chem- 
istry, as  well  as  the  corresponding  preliminary  Widhopf  2-D 
result  based  on  1977  chemistry;  and  (b)  profiles  A and  B which 
were  found  to  be  more  consistent  with  the  balancing  of  the  ozone 
and  CO  budgets.  It  is  noted  that  the  Crutzen  and  Widhopf  pro- 
files based  on  1976  chemistry  are  closer  to  lO-'1'0,  which  was 
the  extreme  high  value  used  in  the  Fishman-Crutzen  calculations 
(Fig.  12).  Hence,  the  calculated  NO  concentrations  in  the  upper 

X 

troposphere  at  middle  northern  latitudes  based  on  1976  chemistry 
were  significantly  larger  than  the  estimated  profile  B at  this 
latitude.  Figure  13  includes  the  1976-Crutzen  N0x  profile  at 
15°N  to  show  the  effect  of  latitude  on  it,  and  which  is  to  be 
compared  with  profile  A of  the  Fishman-Crutzen  2-D  calculations. 
Again,  the  profile  at  15°N  is  significantly  larger  than  profile 
A.  However,  the  result  at  40°N  from  the  preliminary  Widhopf 
calculations  based  on  1977  chemistry  shows  rather  good  agreement 
with  profiles  A and  B,  which  are  applicable  to  latitudes 
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larison  of  N0X  calculations  in  the  troposphere  by  2-D  models  utilizing  1976  (Crutzen 
Widhopf)  and  revised  (Widhopf  and  Fishman-Crutzen)  chemistries. 


0 s 15°N  and  0 s 45°N,  respectively.  It  should  be  noted  that 
the  N0x  volume  mixing  ratio  at  the  surface  for  these  preliminary 
calculations  has  dropped  by  a factor  of  about  five  from  the  cor- 
responding 1976  Widhopf  value.  The  reason  for  this  effect  is 
a change  in  the  eddy  coefficient  at  the  surface,  which  had  to 
be  modified  by  a revision  of  the  modeling  of  the  rainout  sink. 
Therefore,  the  use  of  the  1977  Howard-Evenson  rate  coefficient 
and  the  revised  rainout  in  the  Widhopf  2-D  model  of  the  strato- 
sphere and  troposphere  yield  N0x  values  that  are  in  relatively 
good  agreement  with  the  average  NO  amount  based  on  the  0^  and 
CO  budgets,  which  is,  in  turn,  consistent  with  the  available 
limited  Noxon  data. 

Other  basic  parameters  of  interest  for  2-D  photochemical 
models  are  the  heterogeneous  removal  rates,  not  only  of 
(Table  10),  but  also  those  for  N0x,  HN03,  N^,  and  CH^H. 

These  mechanisms  have  been  parameterized  differently  in  the 
Crutzen  and  Widhopf  2-D  models.  The  former  utilized  crude 
estimates  as  given  in  Hidalgo  and  Crutzen  (1977),  while  the 
latter  did  so  using  as  a reference  the  removal  of  water  vapor 
in  the  lower  troposphere. 

The  parameterization  of  the  water  vapor  distribution  in 
the  reference  stratosphere  is  of  primary  concern  for  the  effect 
of  H^O  effluents  on  ozone,  because  of  the  relatively  dry  strat- 
ospheric medium  in  the  regions  of  significant  ozone  production 
(Figs.  5,  6).  Because  of  the  water  vapor  transports  between 
the  stratosphere  and  troposphere,  the  parameterization  of  the 
water  vapor  in  the  stratosphere  has  to  consider  the  lower 
atmosphere  (0  < z < 50  km)  as  a whole.  The  parameterization 
of  the  tropospheric  water  vapor  in  the  Crutzen  2-D  model  is 
based  on  statistical  data  for  the  2-D  dynamics  as  well  as  the 
mean  temperature  [T]  and  its  standard  deviation  [/T*2]  as  de- 
rived from  observations  by  Oort  and  Rasmusson  (1971)  and  Louis 
(197*0-  The  distribution  of  water  vapor  in  the  troposphere  is 


dependent  on  condensation  processes  and  water  removal  by  rainfall. 
The  basic  procedure  consists  in  the  calculation  of  a local  prob- 
ability of  cloud  formation,  [Pc],  at  each  grid  point  so  as  to 
allow  rainfall  anywhere  within  the  lower  troposphere.  The  pro- 
cedure is  as  shown  schematically  in  Pig.  14  and  it  consists  of 
the  following  steps  (Crutzen,  1975):  (a)  calculation  of  the 

humidity  [q]  from  the  continuity  equation  for  H^O  subject  to 
appropriate  boundary  conditions  at  the  earth's  surface  and  the 

stratopause,  (b)  consideration  of  the  saturation  humidity  [q  1, 

s 

as  a function  of  pressure  and  the  saturation  temperature  [T  ] 

s 

as  given,  for  example,  by  the  Smithsonian  meteorological  tables, 
and  (c)  a Gaussian  distribution  for  the  temperature  at  each 
grid  point  based  on  the  mean  temperature  and  its  standard  de- 
viation. From  the  intersection  of  the  calculated  humidity,  [q], 
with  the  saturation  humidity,  [q  ],  it  is  possible  to  determine 

the  saturation  temperature  [T  ].  Hence,  for  all  temperatures 

s 

[T]«[T  ],  the  air  will  be  supersaturated.  The  probability  of 
s 

cloud  formation  [P  ] can  then  be  determined  from  the  Gaussian 
temperature-distribution  function  at  each  grid-point,  as  indi- 
cated in  Fig.  14.  By  this  method,  it  is  also  possible  to  deter- 
mine the  rate  of  precipitation  at  each  grid  point.  With  know- 
ledge of  the  probability  of  cloud  formation,  it  becomes  possible 
to  estimate  the  loss  rate  [L  ] of  water-soluble  gaseous  compounds 
in  the  tropospheric  sinks.  These  estimates  were  based  on  the 
following: 


CLr(Xl)]  = - k°he  Pc[Xll  = ~ 


Vx,] 


(k°het>i 


-1 


where  the  (k°^et )^_1  were  assumed  to  be  given  by  1.5  days  for 
N0x  and  1/3  day  for  HNO^,  H000,  N2°5*  and  (Hidalgo  and 

Crutzen,  1977).  It  should  be  noted  that  this  value  for  Ho0,, 
is  consistent  with  that  estimated  subsequently  from  considera- 
tions of  the  CO  budget  (Table  10). 
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HUMIDITY,  [q] 


/SATURATED  HUMIDITY 

[qs]=  f (p,[ts]) 


PROBABILITY, 


1 


The  parameterization  of  the  tropospheric  water  vapor  in 
the  Widhopf  model  is  also  based  on  the  continuity  equation  for 
water  vapor  with  an  appropriate  2-D  dynamics,  tropospheric 
sinks,  and  boundary  conditions  at  both  the  ground  and  strato- 
pause.  The  boundary  conditions  at  the  ground  specify  the  rela- 
tive humidity.  Instead  of  using  the  local  probability  of 
cloud  formation  and  precipitation  rates  at  each  grid  point,  as 
in  the  Crutsen  model,  the  Widhopf  precipitation  rates  are 
assumed  to  be  average  values  in  the  troposphere,  and  this  cal- 
culation utilizes  the  concept  of  residence  time,  which  is  de- 
fined as  a function  of  latitude  by  the  ratio  of  the  amount  of 

2 

precipitable  water  in  a vertical  column  of  air  (g/cm  ) to  the 

2 

average  rate  of  precipitation  (g/cm  - year).  These  residence 
times  are  based  on  the  Bannon  and  Steele  (1957)  maps  of  preci- 
pitable water,  which  yield  the  results  given  in  Table  11  (Junge, 
1963): 

TABLE  11.  AVERAGE  RESIDENCE  TIME  OF  WATER  VAPOR  IN  THE  TROPOSPHERE 
AS  A FUNCTION  OF  LATITUDE 

(Source:  Junge,  1963) 

Latitude  (degrees) 

80-90  70-80  60-70  50-60  40-50  30-40  20-30  10-20  0-10 

Residence 

Time  (days)  (15.0)  (13.4)  8.7  6.2  6.4  8.7  12.0  11.2  8.1 


The  results  in  this  table  are  given  as  a function  of  lati- 
tude, and  they  indicate  that  the  average  residence  time  of  water 
vapor  in  the  troposphere  is  about  10  days.  The  calculation  of 
the  precipitation  rates  in  the  Widhopf  model  is  then  given  by 
the  ratio  of  the  local  amount  of  precipitable  water  to  the  resi- 
dence time  for  the  water  vapor  column  as  a function  of  latitude. 
The  parameterization  of  the  heterogeneous  removal  rates  of  water- 
soluble  species  was  assumed  to  be  given  in  terms  of  the  residence 
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times  for  water  vapor  (Table  11).  The  residence  times  for  NO 
and  H202  were  specified  as  being  ten  times  that  of  water  vapor, 
whereas  those  for  HNO^  and  N20,-  were  assumed  to  be  identical  to 
those  of  water  vapor  (Widhopf,  1978). 

Figures  15  and  16  show  the  water  vapor  distributions  as  a 
function  of  latitude  during  the  fall  season  as  given  by  the 
Crutzen  and  Widhopf  models,  respectively.  The  latter  figure 
includes  the  corresponding  Oort  and  Rasmusson  (1971)  data  for 
the  troposphere  and  the  Harries  (1976)  data  for  the  stratosphere. 
These  results  indicate  that  either  parameterization  gives  an 
adequate  description  of  the  water  vapor  distribution  in  the 
troposphere.  In  the  reference  stratosphere,  both  models  give 
similar  results  of  a relatively  dry  and  nearly  constant  verti- 
cal gradient  of  water  vapor  mass  mixing  ratio  at  every  latitude. 
This  parameterization  of  water  vapor  has  so  far  been  tested  in 
an  atmosphere  perturbed  by  H20  effluents  only  in  the  Widhopf 
2-D  model  (Widhopf  and  Glatt,  1978). 


3.4  CRITICAL  CHEMICAL  REACTIONS  FOR  AVIATION  EFFECTS 


The  previous  section  identified  the  five  most  important 
reactions  involving  hydrogen-oxide  species  for  the  assessment 
of  SST  effects  on  the  ozone  column  (Duewer  et  al.,  1977).  In 
chronological  order  of  the  revision  of  their  reaction  rates  they 
are  as  follows: 


(R15)  OH  + N02(+M)  -*•  HNG3(+M) 
(R21)  OH  + H02  - H20  + 02 

(R14)  H02  + NO  -<■  OH  + N02 

(R13)  03  + H02  -►  OH  + 2 02 

(R  7)  03  + OH  H02  + 02 


log1Qk  = -AT/(B+T)-0.5  log1Q(T/280) 
3 x lCf 11 
8 x 10-12 

7.3  x lO-1*1  exp  (-1275/T) 

1.5  x 10-12  exp  (-1000/T) 


where  the  corresponding  reaction  rate  coefficients  are  those 
recommended  most  recently  by  NASA  (1977).  The  parameters  A and 


.... 
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CRUTZEN  2-D  MODEL  H.O  DISTRIBUTION  DURING  FALL  (OCTOBER  30) 


FIGURE  15.  Simulation  of  water  vapor  distribution  in  the  reference  atmosphere  by  the  Crutzen 
model  as  a function  of  latitude  in  the  Northern  Hemisphere  (Source:  Hidalgo  and 
Crutzen,  1977) 


FIGURE  16.  Simulation  of  water  vapor  distribution  in  the  reference  atmosphere  by  the  Widhopf  model 
as  a function  of  latitude  in  the  Northern  Hemisphere.  (Source:  Widhopf,  1978) 


as  a third  body. 

In  addition  to  the  foregoing  set  of  five  critical  reactions 
for  SST  effects,  NASA  (1977)  recommendations  include  two  other 
reactions  that  deserve  top  priority  consideration.  These  reac- 
tions and  rate  coefficients  are  as  follows: 


(R12)  0 + H02  -►  OH  + Op  3.5  x 10-11 

(R17)  H02  + HOp  - H202  + 02  2.5  x 1CT12 

The  NASA  recommendations  are  based  on  the  grounds  that  existing 
data  for  R21,  Rl4,  and  R17  are  mutually  inconsistent;  while  R13, 
R7,  and  R12  contribute  to  the  total  uncertainty  of  stratospheric 
models.  It  is  noted  that  several  of  the  above  rate  coefficients 
have  omitted  any  temperature  dependence  that  may  exist  in  such 
rates,  because  of  the  lack  of  such  measurements  at  the  time  of 
the  NASA  recommendations. 

Preliminary  estimates  have  been  made  of  the  sensitivity  of 
the  rate  coefficients  of  the  reactions  R21,  R13,  and  R17  for 
changes  of  the  stratospheric  ozone  column  caused  by  aircraft 
effluents  (Crutzen  and  Howard,  1978).  The  range  of  uncertainty 
in  these  rate  coefficients  were  considered  to  be  given  by 

2.1  x 10'11  * k21  * 5*1  x 10"11 

10_13exp  (-1525/T)  * k13  * 2 x 10_13exp  (-1250/T) 

5 x 10_12exp  (-500/T)  s k1?  s 1.5  x 10_11exp  (-500/T) 

which,  for  example,  at  room  temperature  (296°K)  yield 

5.8  x 10~16  s k,,  * 2.9  x 10"15  and  9-2  x 10-13  £ k s 2.8  x 

-12  1 < 

10  . The  temperature  dependency  for  the  Howard-Evenson  reac- 

tion rate  was  assumed  to  be  given  by 
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I 


I 


= 3 x 10“  exp  (-390/T) 


These  sensitivity  estimates  considered  combinations  of  the 
above  end  values  for  each  range  using  (a)  a steady-state,  one- 
dimensional model  that  did  not  include  the  troposphere,  (b)  an 


oxygen-hydrogen-nitrogen  stratosphere,  and  (c)  assumed  NO  and 

o c*  ^ 


o cr  ? 

H^O  emissions  of  10  and  2 x 10  molecules/cm  -sec , respectively, 
at  altitudes  of  either  16  or  20  km.  Although  these  preliminary 
investigations  could  not  overcome  the  problem  of  underdestruc- 
tlon  of  ozone  in  the  reference  stratosphere  (Section  3.1),  the 
results  indicated  that  the  sign  of  the  change  of  the  ozone  col- 
umn by  SST  flight  at  16  km  or  20  km  depended  mainly  on  the 
relative  magnitudes  of  the  rate  coefficients  for  R21  and  R13* 
These  results  are  summarized  in  Table  12.  The  fourth  column 
in  this  table  shows  that  the  sign  in  the  change  of  the  ozone 


column  by  the  assumed  NO  and  H„0  effluents  at  either  16  km 

A (— 


or  20  km  changes  from  negative  to  positive  for  a combination 

S'  / n -i  „ i a™  11  \ ...  i mi  . . j i.  i j i m 


of  the  lower  (2.1  x 10  ) value  for  R21  with  the  higher  value 

(2.9  x 10“15)  of  R13.  The  higher  value  of  R13  is  an  upper  limit 
derived  in  the  original  sensitivity  analysis  by  Duewer  et  al . , 
1977.  This  table  includes  also  the  NASA  recommended  values 
for  these  rates,  where  the  room  temperature  value  has  again 
been  used  for  R13  to  provide  a reference  when  comparing  it  with 
the  values  in  these  preliminary  investigations.  In  spite  of  the 
NASA  recommendations  for  the  rate  coefficient  of  R13>  recent 
(1978)  experimental  results  by  Howard  indicate  that  this  rate  is 
actually  closer  to  the  higher  value  used  in  these  calculations.* 


As  indicated  earlier,  recent  measurements  of  the  temperature 
dependence  of  the  reaction  rate  coefficient  of  R13  is 

k13=(1.4  ± 0.4)  10-1\xp  [-(580  ± 100 )/T]  cm^  molecule-1  sec-1 

(Howard,  1978) . At  296°K,  this  rate  coefficient  yields 

k13  = 20  x 10"16. 


2-1(-ll)  lu  3.0(-ll) 

29.0(-16)  9.8(-16) 


It  must  be  emphasized  that  these  preliminary  results  are 
subject  to  the  additional  constraints  of  a lack  of  coupling 
between  the  N0x  and  C10x  cycles  in  the  stratosphere  (M2,  Table  1) 
below  35  km  through  reactions  such  as  NOp  +C10  + M CINO^  + M 
and  NO  + CIO  Cl  + NO^.  Furthermore,  a current  open  question 
in  the  chlorine  chemistry  is  the  difference  by  an  order  of 
magnitude  between  the  higher  observed  and  (1-D)  C10x  values  in 
the  reference  atmosphere  (Crutzen  et  al.,  1978).  It  is  current- 
ly unknown  if  this  difference  is  due  to  experimental  errors  in 
the  observations  or  to  an  unknown  source  of  C10x  in  the  strato- 
sphere . 
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4.  ASSESSMENT  OF  THE  2-D  FORMULATION 
OF  THE  ATMOSPHERIC  DYNAMICS 


I 


The  previous  section  emphasized  considerations  of  the 
chemical  mechanisms  in  Fig.  2 and  Table  1 for  the  reference  and 
perturbed  stratospheres  and  tropospheres.  Because  of  the  empi- 
ricism in  the  2-D  formulation  of  the  dynamics  for  the  reference 
atmosphere,  it  becomes  necessary  to  (a)  review  the  basic  assump- 
tions in  the  available  2-D  parameterizations  of  the  dynamics  for 
the  reference  upper  troposphere  and  stratosphere,  (b)  assess  the 
applicability  of  such  formulation  to  extensions  of  such  para- 
meterizations to  an  atmosphere  assumed  to  be  perturbed  by  large 
continuous  injections  of  aircraft  effluents  in  the  flight  cor- 
ridors, and  (c)  explore  newer  concepts  for  a more  correct  theo- 
retical 2-D  formulation  of  the  dynamics  for  the  reference  and 
perturbed  upper  tropospheres  and  lower  stratospheres. 

4.1  BASIC  ASSUMPTIONS  IN  THE  2-D  PARAMETERIZATIONS  OF  THE 

DYNAMICS 

A basic  aim  in  the  2-D  formulation  of  the  dynamics  of  photo- 
chemical models  is  to  describe  the  mass  transports  of  trace  chem- 
ical species  in  the  reference  lower  atmosphere  (0  < z < 50  km) 
as  well  as  the  enhanced  mass  transports  of  engine  effluents  for 
assumed  large  aircraft  traffic  in  the  flight  corridors  of  the 
perturbed  upper  troposphere  and  lower  stratosphere.  It  should 
be  noted  that,  because  of  the  relatively  long-time  scales  in 
Fig.  3 (i.e.,  months),  the  interest  is  on  the  eddy  mass  transports 
by  the  large-scale  motions.  The  available  2-D  parameterizations 
of  the  atmospheric  dynamics  have  utilized  two  basic  assumptions 
which  are  (a)  the  decoupling  of  the  temperature  and  dynamics 
from  the  chemistry,  and  (b)  the  extension  of  the  mixing  length 


hypothesis  for  momentum  transports  in  microscale  turbulent  flow 
to  the  mass  transports  in  the  macroscale  turbulence  of  the 
stratosphere.  Each  of  these  assumptions  is  described  below. 

4.1.1  Decoupling  of  Temperature  and  Dynamics  from  the  Chemistry 

A first  basic  approximation  in  the  2-D  formulation  of  the 
Crutzen  and  Widhopf  2-D  models  is  to  decouple  the  conservation 
of  the  chemical  species  from  the  so-called  primitive  equations 
for  the  general  circulation  of  air  in  the  reference  troposphere 
and  stratosphere.  Since  the  primitive  equations  for  the  refer- 
ence atmosphere  essentially  describe  the  temperature  and  dynamics, 
a fundamental  assumption  in  these  2-D  models  is  then  that  anthro- 
pogenic changes  in  the  distributions  of  radiative  and  chemically 
active  trace  gases  (e.g..  Fig.  6)  will  not  produce  first-order 
effects  on  the  initial  temperature  field  and  air  motions  of  the 
reference  atmosphere.  Hence,  such  decoupled  formulation  cannot 
consider  feedback  effects  on  the  calculated  ozone  changes  from 
any  change  in  the  background  stratospheric  temperature  and  dy- 
namics which  would  result  from  the  changes  in  both  the  absorp- 
tion of  solar  radiation  and  emission  of  infrared  radiation  that 
would  be  induced  by  the  calculated  changes  of  ozone  (Figs.  7, 

10).  These  feedback  effects  on  the  ozone  change  would  take  place 
through  (a)  the  temperature  dependence  of  the  reaction  rate  co- 
efficients (e.g..  Table  3),  (b)  the  coupling  between  the  tempera- 
ture changes  and  changes  in  the  air  motions,  and  (c)  changes  in 
the  distributions  of  chemical  species  due  to  the  changes  in  the 
dynamics.  Obviously,  the  validity  of  the  decoupling  of  the  tem- 
perature and  dynamics  from  the  chemistry  becomes  increasingly  open 
to  question  as  the  magnitude  of  the  changes  in  ozone  by  all  anth- 
ropogenic effects  becomes  large,  a condition  that  stems  from  the 
fact  that  such  anthropogenic  effects  on  the  temperature  and  dyna- 
mics would  become  part  of  the  atmospheric  background  for  the  air- 
craft effluents.  In  this  context,  it  is  of  interest  to  note  that 
current  assessments,  based  on  1977  chemistry,  of  the  effect  on 


the  global  (1-D)  ozone  column  from  recent  release  rates  of 
chlorof luoromethanes  at  the  earth's  surface  (Pig.  2 as  well  as 
M17  and  M2  in  Table  1)  forecast  an  eventual-  decrease  of  the  1-D 
ozone  column  by  as  much  as  about  1*4  percent  by  the  first  quarter 
of  the  next  century  (e.g.,  NASA,  1977). 

The  coupling  between  the  temperature  field  and  dynamics 
with  the  chemistry  in  the  primitive  equations  for  the  reference 
stratosphere  is  taken  into  account  by  a parameterization  of  the 
net  heating  of  the  stratospheric  air.  Hence,  the  remaining 
coupling  between  the  temperature  and  dynamics  can  be  illustrated 
by  considering  the  change  in  the  magnitude  of  the  vertical  eddy 
transports  in  the  tropopause  region,  which  are  of  crucial  impor- 
tance in  the  assessment  of  the  magnitude  of  aviation  effects  on 
ozone.  The  change  in  the  vertical  eddy  transports  in  the  tropo- 
pause region  becomes  important  because  engine  effluents  deposited 
above  the  tropopause  are  then  trapped  and  allowed  eventually  to 
propagate  upwards  to  the  ozone  production  regions  in  the  middle 
and  upper  stratosphere  (Fig.  5)* 

The  change  in  the  magnitude  of  the  vertical  eddy  transports 
in  the  tropopause  region  depends  on  the  change  in  atmospheric 
stability,  which  depends,  in  turn,  on  the  change  in  the  sign  of 
the  vertical  temperature  gradient  in  this  region.  The  latter 
change  is  due  to  the  vertical  distribution  of  the  ozone  mixing 
ratio  in  the  stratosphere  (Fig.  6)  and  the  consequent  absorp- 
tion of  UV  solar  radiation  by  ozone.  As  indicated  by  the  U.S. 
Standard  Atmosphere  (1976),  for  example,  the  undisturbed  am- 
bient temperature  (a)  decays  with  Increasing  altitude  in  the 
troposphere,  (b)  becomes  constant  at  the  tropopause  (z^)  and 
in  the  lower  stratosphere  (z^  < z < 20  km),  and  (c)  increases 
with  increasing  altitude  in  the  middle  and  upper  stratosphere 
(20  < z < 50  km).  A fundamental  problem  in  the  2-D  description 
of  the  dynamics  is  then  the  formulation  of  the  eddy  transports 
of  chemical  species  between  the  lower  stratosphere  and  upper 
troposphere.  However,  basic  difficulties  in  such  formulation 
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are  (a)  the  lack  of  observations  for  the  statistics  of  the  very 
small,  vertical  transports  in  the  tropopause  and  lower  strato- 
sphere regions  (Oort  and  Rasmusson,  1971),  and  (b)  the  complex 
three -dimens  ioyial  dynamics  of  the  tropopause  gaps  produced  by 
the  jet  streams  at  the  northern  middle  latitudes  which  are  im- 
portant for  aviation  effects  (Fig.  7). 

The  most  general  approach  for  the  simulation  of  the  inter- 
actions of  the  temperature  field,  the  dynamics  and  chemical 
processes  in  the  reference  and  perturbed  atmospheres,  would  be 
based  on  the  coupling  of  the  3-D  primitive  equations  for  atmos- 
pheric motions  with  those  for  the  conservation  of  chemical 
species.  This  coupled  system  of  equations  would  yield  solutions 
for  the  vertical  and  horizontal  transports  of  chemical  species  as 
a function  of  time  and  the  three  space  coordinates  of  the  atmos- 
phere (latitude,  altitude,  longitude).  However,  the  complexity 
of  such  a system,  coupled  with  the  basic  limitations  given  in 
Section  1,  has  so  far  prevented  the  practical  application  of  such 
a system  to  either  the  reference  stratosphere  or  the  study  of  the 
effects  of  high-altitude  subsonic  and  SST  flight  on  atmospheric 
ozone.  Nevertheless,  it  is  of  interest  to  describe  such  a system 
in  order  to  bring  out  the  subsequent  2-D  approximations. 

The  3-D  or  primitive  equations  for  the  air  motions  without 
chemical  interactions,  in  a frame  of  reference  fixed  to  the  ro- 
tating earth,  are  as  follows  (e.g.,  Lorenz,  1967): 


dU/dt  = - fk  x 5 - g?z  + f (1) 
dT/dt  = kToi/p  + Q/Cp  (2) 
$ • U + 3w/3p  = 0 (3) 
3z/3p  = -RT/gp  (4) 


Equation  (1)  represents  the  conservation  of  momentum  for  the 
quasi-horizontal  motions  (U)  and  contains  two  scalar  equations 
for  the  eastward  (u)  and  northward  (v)  air  velocities.  This 
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equation  gives  the  total  or  incJ  vidual  acceleration,  defined  by 
d)(/dt  = + ^ • x U + 3(x^)/3p,  in  terms  of  the  Coriolis 

acceleration  due  to  the  earth's  rotation  (-fk  x U),  that  due  to 
the  geopotential  (gz)  gradients  along  surfaces  of  constant  pres- 
sure (p),  and  the  acceleration  due  to  friction  (?).  Equation  (2) 
represents  the  conservation  of  energy  and  gives  the  total  time 
rate  of  change  of  the  air  temperature  (T)  on  the  left-hand  side, 
whereas  the  variables  cj(=dp/dt)  and  Q are  related  to  the  vertical 
winds  and  the  net  heating  of  air  per  unit  mass,  respectively. 

In  the  relatively  dry  stratosphere  (Fig.  6),  the  net  heating  is 
the  result  of  (a)  absorption  of  solar  radiation  by  ozone,  carbon 
dioxide,  and  water  vapor;  and  (b)  emission  of  IR  radiation  by  0^, 
C02,  and  H20,  which  is  a mechanism  dominant  in  the  lower  strato- 
sphere. Equation  (3)  represents  the  conservation  of  air  mass 
and  Eq.  (4)  the  condition  of  vertical  hydrostatic  equilibrium, 
which  allows  conversion  of  the  altitude  to  the  pressure  coordi- 
nate. The  condition  of  hydrostatic  equilibrium  filters  out 
sound  waves  that  are  unimportant  for  the  macroscale  motions, 
whereas  the  condition  of  geostrophic  equilibrium  [i.e.,  balance 
of  the  Coriolis  acceleration  with  those  from  pressure  gradients 
In  Eq . (1)]  at  middle  and  high  latitudes  filters  out  the  unim- 
portant gravity  waves  (Lorenz,  1967). 


The  interactions  of  the  stratospheric  temperature  and  dy- 
namics with  the  chemistry  are  given  by  the  conservation  of  the 
chemical  species;  i.e.. 


(5) 


where  R^  denotes  the  mixing  ratio  for  any  of  the  given  i species, 
and  as  well  as  the  production  and  loss  chemical  mechanisms 
In  Table  1.  It  should  be  noted  that  Eq.  (5)  is  also  applicable 
to  phase  changes  such  as  those  of  H20  in  the  troposphere,  and 
that  the  production  term  for  aviation  effects  must  include  the 
engine  effluents.  The  dynamic  terms  in  Eq.  (5)  may  be  identified 


through  the  expansion  of  dF^/dt,  which  in  pressure  coordinates 
becomes  (Lorenz,  1967): 


where 


9ojR  . 
+ 9p“ 


(6) 


and  "a”  denotes  the  earth's  radius. 

The  primitive  equations  (1  through  4),  without  the  use  of 
the  conservation  of  chemical  species,  are  the  basis  for  the 
available  general  circulation  models  (GCMs)  of  the  stratosphere 
and  troposphere.  The  GCMs  must  consider  the  conservation  equa- 
tion for  H^O  and  extend  the  parameterization  of  the  net  heating 
(Q)  to  include  the  heating  by  condensation  of  H^O  in  the  lower 
and  middle  troposphere.  The  GCMs  must  incorporate  further  para- 
meterizations  for  the  effect  of  clouds  on  the  incoming  solar 
radiation  as  a function  of  season,  the  sinks  in  the  tropo- 

sphere, and  appropriate  boundary  conditions  at  both  the  earth's 
surface  and  stratopause. 

Either  observations  or  solutions  of  the  primitive  equations 
provide  data  for  the  foregoing  dependent  variables  x(^>  ©>  P,  t) 
as  a function  of  longitude  (X),  latitude  (0),  pressure  level  (p), 
or  height  (z),  and  the  time  (t).  However,  such  variables  may  be 
expressed  as  x = [ x 3 + X*>  where  the  operation  [ ] denotes  longi- 
tudinal or  zonal  averages  and  x*  deviations  from  such  averages. 
Hence,  the  zonal  mean  is  given  by  [ x D = f(0,  p,  t);  whereas  x*  = 
g(A,  0,  p,  t)  but  Cx* 1 = 6*  Similarly,  the  dependent  variables 
may  also  be  expressed  in  terms  of  time  averages  x = X + X1 » where 
X and  x'  are  associated  with  the  so-called  standing  and  transi- 
ent motions,  respectively.  Hence,  any  dependent  variable  such 
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as  the  northward  wind  (v)  may  be  expressed  as  v = [v]  + [v]'  + 
v#  + v*'.  The  poleward  transport  of  x (say  0^)  is  then  given 
by  (Lorenz,  1967): 


Lxv]  = ExHv]  + Ixl'M’  + [x*v*  ] + [x* ' V*  ’ ] (8) 


where  the  first  two  terms  on  the  right  side  are  the  standing 
and  transient  transports  by  the  meridional  circulation,  and  the 
latter  two  the  corresponding  standing  and  transient  eddy  trans- 
ports. This  split  of  the  atmospheric  motions  considers  a zon- 
al symmetric  flow  plus  eddy  or  wave  motions  superposed  on  the 
zonal  symmetric  flow.  The  wave  motions  are  characterized  by  a 
spectrum  of  (for  example)  zonal  wave  numbers,  i.e.,  number  of 
wavelengths  in  a latitudinal  circle  around  the  earth's  axis  of 
rotation. 


A general  2-D  formulation  of  the  couplings  of  the  tempera- 
ture, dynamics  and  chemistry  would  involve  taking  the  time  and 
zonal  averages  of  the  system  of  Eqs . (l)  through  (5).  However, 
from  the  decoupling  of  the  chemistry  from  the  temperature  and 
dynamics,  the  basic  2-D  equation  for  the  conservation  of  species 
may  then  be  obtained  by  (a)  performing  the  zonal  average  of 
only  Eq.  (5),  (b)  using  Eq.  (8)  to  split  the  transports  into 
the  zonal  mean  and  eddy  components,  and  (c)  using  the  zonal 
average  of  the  continuity  Eq . (3).  The  result  is  as  follows:* 


acy 

3t“~ 


30^] 

30“ 


+ [«] 


30^] 

W 


— — pr  [v*R,  * ]cos0 

a cos0  30  L i 


3p 


- [pp  - [4] 


(9) 


The  decoupling  of  the  chemistry  from  both  the  air  tempera- 
ture and  motions,  as  used  in  the  Crutzen  and  Widhopf  2-D  models, 
eliminates  the  need  to  consider  the  time  and  zonal  averages  of 
the  primitive  equations  (1)  through  (3)  for  the  reference  and 

u 

The  time-averaging  notation  in  Eq . (8)  is  subsequently  dropped 
for  simplicity.  Likewise,  the  difference  between  standing  and 
transient  eddies  is  omitted. 
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perturbed  atmospheres.  These  2-D  models  must  instead  specify  as 
a function  of  month  or  seasoxi,  the  meridional  air  temperature 
(so  as  to  determine  the  reaction  rate  coefficients,  e.g.  , 

Table  3)  and  motions  for  the  reference  stratosphere  and  tropo- 
sphere. The  role  of  the  2-D  formulation  of  the  dynamics  is 
therefore  defined  by  the  second  to  fifth  terms  in  the  left-hand 
side  of  Eq.  (9).  The  second  and  third  terms  involve  the  trans- 
ports of  mixing  ratio  by  the  meridional  circulation,  whereas 
the  fourth  and  fifth  are,  respectively,  the  poleward  and  verti- 
cal mass  transports  by  the  eddy  motions.  Hence,  the  first  term 
in  the  left-hand  side  of  Eq.  (9)  is  the  net  local  rate  of  change 
of  the  mixing  ratio  of  any  of  the  i species  as  a result  of  the 
transports  by  the  dynamics  as  well  as  the  chemical  mechanisms 
in  the  right-hand  side  of  the  equation.  It  is  important  to 
note  that  the  determination  of  the  meridional  profiles  of  ozone 
and  related  i species  (e.g.,  Fig.  6)  requires,  in  general,  the 
simultaneous  numerical  solution  of  a set  of  i equations  (9). 

While  the  Widhopf  2-D  model  solves  for  each  individual  Eq . (9), 
the  Crutzen  model  utilizes  equilibrium  approximations  for 
species  with  relatively  short  characteristic  times  so  as  to 
reduce  the  number  of  equations.  These  simultaneous  solutions 
must  utilize  appropriate  boundary  conditions  at  the  earth's 
surface  and  stratopause,  and  include  the  heterogenous  mechanisms 
in  the  tropospheric  sinks  (M9,  Table  1)  in  the  term  for  the 
lower  troposphere. 

4.1.2  Hypothesis  of  Large-Scale  Diffusion  in  the  Stratosphere 

The  decoupling  of  the  chemistry  from  the  air  temperature 
and  motions  allows  the  simultaneous  solution  of  the  system  of 
Eq.  (9)  with  specified  meridional  temperature  profiles  and  dy- 
namics for  the  reference  atmosphere  as  a function  of  month  or 
season.  The  dynamics  is  given  by  (a)  the  meridional  circulation, 
given  in  Fig.  17,  which  was  determined  by  Louis  (1974)  and 
Louis  et  al.  (1974)  from  temperature  observations  and  the  use 
of  the  energy  equation  as  well  as  the  conservation  of  air  mass; 
and  (b)  the  poleward  and  vertical  eddy  mass  transports,  which 


FIGURE  17.  Meridional  Circulation  as  a function  of  season.  (Source:  Louis  in  CIAP  Monograph  3 
1975;  Chapter  4). 


are  essentially  unknown  in  the  stratosphere  due  to  the  lack  of 
observations  for  the  statistics  of  such  motions  (Oort  and  Ras- 
musson,  1971). 

A second  basic  approximation  in  the  formulation  of  the 

Crutzen  and  Widhopf  2-D  models  is  the  extension  of  a modified 

Prandtl  (1925)  mixing  length  hypothesis  for  the  eddy  transports 

of  momentum  in  microscale  turbulence  to  the  eddy  transports  of 

mass  in  the  macroscale  turbulence  in  the  reference  atmosphere. 

In  the  usual  space  coordinates  of  microscale  flow,  the  poleward 

eddy  flux  in  Eq . (9),  for  example,  would  be  given  by  the  Prandtl 

mixing  length  (Jt)  hypothesis  as  [v*x*3  = -K  9[x3/<*y;  where 

y (=  aO)  is  the  horizontal  coordinate,  taken  as  positive  (Ay>0) 

in  the  poleward  direction.  The  eddy  diffusion  coefficient , 

K (y,  z,  t),  would  be  a specified  function  of  latitude  (y  or 
j <y 

0),  altitude  (z),  and  month  or  season  (t).  The  Prandtl  mixing 
length  concept  would  thus  eliminate  the  eddy  fluxes  in  Eq.  (9) 
in  terms  of  latitudinal  and  vertical  gradients  of  the  mixing 
ratio  with  [ x ] = [R^].  However,  the  Prandtl  mixing  length 
hypothesis  breaks  down,  for  example,  in  the  lower  stratosphere 
when  [x]  denotes  either  ozone  or  temperature;  i.e.,  for  the  eddy 
fluxes  of  ozone  ([v*0*^])  or  heat  ([v*T*]).  Since  the  eddy 
diffusion  coefficient  is  always  positive  (Kyv>0),  the  reason 
for  the  breakdown  of  the  Prandtl  mixing  length  hypothesis  in 
the  stratosphere  is  that  9[x3/9y  would  have  to  be  negative  for 
a positive  poleward  transport  of  0^  or  temperature;  i.e.,  con- 
trary to  observations,  the  latitudinal  distribution  of  [ x 3 along 
surfaces  of  constant  pressure  would  have  to  decrease  with  in- 
creasing latitude.  Thus,  the  Prandtl  mixing  length  hypothesis 
requires  that  the  direction  of  the  eddy  flux  must  always  be 
from  high  to  low  values  of  [x3- 

The  application  of  the  Prandtl  mixing  length  concept  had 
then  to  be  modified  If  it  was  to  take  into  account  poleward 
countergradient  eddy  fluxes  in  a direction  from  low  to  high 


values  of  [ x ] - The  countergradient  eddy  fluxes  In  the  lower 
stratosphere  were  explained  by  Newell  (196*0  and  formulated 
analytically  by  Reed  and  German  (1965).  The. compression  of  the 
atmosphere  with  increasing  latitude  towards  the  cold  polar  re- 
gions produces  a downward  inclination  in  the  poleward  slopes  of 
constant  pressure  surfaces.  Because  the  temperature  is  constant 
in  the  relatively  dry  (Fig.  6)  undisturbed  lower  stratosphere, 
the  slope  of  constant  entropy  tends  then  to  be  the  same  as  the 
downward,  poleward  slopes  of  the  isobaric  surfaces.  The  crux 
of  the  argument  to  explain  countergradient  eddy  fluxes  of  heat 
in  the  lower  stratosphere  was  that  the  subsiding  air  parcels 
moving  northward  from  the  cold  (altitude-wise)  tropical  regions 
do  so  along  inclined  paths  with  a downward  slope  larger  than 
those  of  the  undisturbed  isentropic  surfaces.  Similarly,  as- 
cending air  parcels  would  be  transported  southward  along  paths 
with  inclinations  similar  to  those  of  the  subsiding  parcels. 

The  foregoing  concept  for  countergradient  eddy  fluxes  is 
illustrated  in  Fig.  18  for  the  heat  flux,  where  [\]  denotes  the 
potential  temperature  as  defined  by  an  assumed  isentropic  com- 
pression from  the  parcel  temperature  and  pressure  to  a standard 
pressure  (1000  mbar).  The  use  of  the  potential  temperature 
separates  the  effects  of  adiabatic  and  nonadiabatic  heating 
during  the  displacement  of  the  air  parcels  across  a mixing 
length  distance  (approximately  100  km).  The  figure  shows  (1) 
a downward  slope  a for  both  subsiding  and  ascending  parcels 
arriving  at  the  point  defined  by  the  origin  of  the  horizontal 
and  vertical  axes,  and  (2)  the  smaller  slope  P for  the  isen- 
tropes.  The  local  gradient  vector  $[xl  in  the  figure  gives 
9[x]/<*y>0  along  the  horizontal  axis;  i.e.,  if  the  eddy  flux  is 
poleward,  it  will  take  place  from  low  to  high  values  of  f x 3 • 

The  countergradient  flux  of  potential  temperature  at  the  origin 
becomes  evident  by  considering  the  net  eddy  transport  of  [v*x*3* 
If  the  eddy  fluxes  of  heat  are  quasi-i sentropic , the  northward 
subsiding  parcels  arriving  at  the  origin  will  be  anomalously 
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FIGURE  18.  Schematic  representation  of  the  slope  (a)  of  the  mixing 
length  path  necessary  for  countergradient  eddy  fluxes  in 
the  Prandtl  mixing  length  hypothesis.  (Source:  Reed  and 
German,  1965) 


warm;  whereas  the  southward  ascending  particles  arriving  at  the 
same  point  will  be  anomalously  cold.  Hence,  the  net  eddy  trans- 
port of  heat  in  the  subsequent  mixing  of  the  air  parcels  will 
be  northward  in  spite  of  the  fact  that  9[x]/3y>0;  i.e.,  the 
poleward  eddy  transports  of  heat  become  countergradient  when 
a > S.  Similar  consideration  indicates  that  when  a = 3,  the 
eddy  fluxes  tend  to  vanish,  and  when  a = 0 the  poleward  eddy 
flux  becomes  as  in  the  Prandtl  mixing  length  hypothesis.  These 
results  are,  therefore,  a function  of  a for  the  fixed  condition 
of  3[X]/3y>0  along  the  y-axis.  Hence,  the  magnitude  of  the 
slopes  a(y,z,t)  of  the  eddy  motions  as  well  as  the  corresponding 
variability  (i.e.,  variance  [a*  ])  of  such  slopes  determine  the 
nature  (with  or  countergradient)  as  well  as  the  magnitude  of 
the  horizontal  and  vertical  eddy  fluxes  in  Eq . (9)-  The  modi- 
fication of  the  Prandtl  mixing  length  hypothesis  must  then  take 
into  account  the  slopes  a and  3 of  the  mixing  length  and  isen- 
tropic  surfaces. 

The  formulation  of  the  eddy  horizontal  and  vertical  trans- 
ports in  Eq.  (9)  are  usually  given  in  the  space  coordinates 
(y,z)  instead  of  those  in  Eq.  (9).  Hence,  the  concepts  illus- 
trated in  Fig.  18  indicate  that 

X*  - -t  • 'txl  - - (*y  ^ + . O0) 

where  1 is  the  mixing  length  vector  with  horizontal  and  vertical 

components  given  by  SL  and  Jt  . The  horizontal  and  vertical  eddy 

«/  ^ 

fluxes  are  then  given  by 

[v*y*]  = - (k  d Ex J + k j.jr-X-l  ^ 

L X J \ yy  oy  + *yz  32  ) 

and 

- - (x2y  iCxi  ♦ kZ2  sc*])  _ 


(11a) 


(lib) 
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where  Kyy  E CV*]*  Kyz  H [V*3’  Kzy  = and  Kzz  = [*zw*]* 

Using  the  approximation  that  the  mixing  length  components  (£  , £ ) 

z y 

are  proportional  to  the  respective  air  velocities  (w,  v)  and  the 
fact  that  the  magnitude  of  a is  small  (approximately  1 0-i|),  then 

v*  = V cos  a = V 

£ = £ cos  a = £ /no 

y (12 

w*  = V sin  a = Va 

£ = £ sin  a = £a  , 
z 

where  V is  the  wind  along  the  mixing  length  direction.  With 

Eqs.  (12),  the  eddy  coefficients  K . K , K , and  K in  the 

yy  yz’  zy’  zz 

Cartesian  coordinates  of  Fig.  18  can  be  expressed  in  the  prin- 
cipal axes  given  by  t and  as  follows: 


yy 

^zy 


K 


zz 


[£yV*]  = [V  £]  = K 

[£zv*]  = [V£<x]  * [V£ ( [a]  + a*)]  * [a]  K 
[£yw»]  = [V£a]  = Kyz 

[£zw*]  = [V£x2]  = [V£([ct]  + a*)2]  * ([a]2  + [a*2])  K, 


(13) 


where  it  is  assumed  that  [a]  and  a*  are  independent  of  V and  ?. 
The  results  in  Eqs.  (13)  indicate  a coupling  among  K , K , and 
K2Z.  Substitution  of  Eqs.  (13)  into  (11a)  and  (lib),  and  using 
[8]  - [tan  8]  = - ( 3 [x ]/9y )/( 3 Cx ]/3z ) , alternative  expressions 
for  the  eddy  fluxes  are  given  by 


[v*X*  ] 


9y 


[w#x*]  = _Kzz  (l 


[a][B]  \ 3 [ x ] 

[a]2  + [a*2]/  9z 


(14) 


Equation  (14)  indicates  that  when  [a]  = [ 8] , the  horizontal 
eddy  fluxes  vanish,  and  when  [a]  = 0,  the  horizontal  and  verti- 
cal eddy  fluxes  are  as  in  the  Prandtl  formulation.  The  Crutzen 
2-D  model  utilized  Eqs.  (13)  in  the  form  (Crutzen,  1975): 
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(15) 
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K11 
[a]  K 

K22  + 
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where  K22  = [a*  1 K^.  The  Crutzen  model  used  the  seasonal  dis- 
tributions of  ozone  and  water  vapor  in  the  lower  stratosphere, 
where  the  effects  of  chemistry  are  not  dominant  (Section  2.1  and 
Pig*  5),  so  as  to  determine  the  eddy  parameters  K,,,  [a],  and 
K22*  a test  of  this  parameterization  of  eddy  coefficients  was 
made  by  simulating  the  late-time  history  of  the  dispersion  of 
gaseous  tracers  deposited  in  the  lower  stratosphere  by  nuclear 
detonations.  The  results  indicated  depletions  of  the  strato- 
spheric burdens  of  such  tracers  that  were  faster  than  observa- 
tions by  about  30  percent  (Crutzen,  1978).  This  discrepancy  was 
not  considered  to  be  significant  in  view  of  the  uncertainty 
associated  with  the  accuracy  of  this  type  of  tracer  data. 


The  Widhopf  2-D  model  utilizes  the  height  coordinate  and 
Eqs.  (11a)  and  (lib).  As  indicated  earlier,  the  Widhopf  para- 
meterization of  the  eddy  coefficients  K , K , and  K was  based 

yy  yz5  zz 

on  (a)  extrapolations  to  the  stratosphere  by  Luther  (1973),  of 
the  Oort  and  Rasmusson  (1971)  data  for  the  eddy  transports  of 
heat  in  the  upper  troposphere  up  to  latitudes  of  50°N.  Such 
extrapolations  have  to  take  into  account  directly  the  complex 
3-D  dynamics  of  the  tropopause  gaps  (Section  4.1.1)  at  the  impor- 
tant middle  latitudes  for  aviation  effects  (e.g..  Fig.  7),  and 
(b)  the  matching  of  segments  of  the  time  history  of  excess  C-14 
data  from  the  1961-1962  nuclear  detonations  in  the  atmosphere 
at  approximately  70°N  latitude,  a procedure  that  is  constrained 
by  the  accuracy  and  scope  of  the  C-l4  data.  More  specifically, 
this  procedure  consisted  in  the  reproduction  of  the  C-14  profiles 
at  30°N  between  April  1963  and  January  1966  with  adequate  accuracy 
(Widhopf  et  al.,  1977).  Subsequent  independent  tests  of  this 
parameterization  included  numerical  simulations  of  the  distri- 
butions of  other  tracers  in  the  stratosphere  such  as  W-I85  from 
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a nuclear  detonation  at  approximately  11°N  latitude  In  the  summer 
of  1958  (Friend  et  al.,  1961)  and  radioactive  Zr-95  from  a nu- 
clear detonation  at  approximately  40°N  in  the  summer  of  1967 
(Telegadas,  197*0.  These  tests  of  the  parameterization  of  eddy 
transports  also  gave  adequate  results  when  compared  with  these 
observations  of  limited  accuracy.  Since  the  Widhopf  parameteri- 
zation of  the  eddy  coefficients  is  thus  independent  of  the  chem- 
istry, the  subsequent  successful  simulation  of  the  meridional 
ozone  distributions  could  be  assumed  to  be  a consistency  test  of 
the  formulation  of  the  1976  chemistry. 

The  magnitude  of  the  eddy  coefficients  is  illustrated  in 
Table  13,  which  shows  the  eddy  coefficients  for  the  Crutzen 
model  for  the  lower  stratosphere  and  upper  troposphere  during 
summer  and  winter  in  the  Northern  Hemisphere.  The  results  in 

Table  13  indicate  that  the  horizontal  eddy  coefficient  (K  ) is 

10  2 7^ 

of  the  order  of  10  cm  /sec,  and  its  magnitude  at  every  latitude 

is  significantly  larger  during  winter  than  in  summer.  The  slope 

-4 

of  the  mean  diffusion  axis  ([a])  is  of  the  order  of  10  , and  it 

is,  in  general,  negative,  as  indicated  in  Fig.  18.  The  eddy 

4 2 

coefficient  K is  of  the  order  of  10  cm  /sec  in  the  lower  strat- 
zz 

osphere,  and  its  increase  with  decreasing  altitude  near  the  tropo- 
pause  is  as  expected  from  the  previous  considerations  of  the  re- 
lative stabilities  of  the  lower  stratosphere  and  upper  troposphere. 

Table  14  illustrates  the  eddy  coefficients  for  the  Widhopf 
model,  which  are  given  as  a function  of  month  instead  of  season, 
as  in  the  Crutzen  model.  Because  of  the  poorer  data  base  for 
the  Southern  Hemisphere,  the  eddy  coefficients  are  taken  to  be 
there  the  same  as  in  the  Northern  Hemisphere;  but  shifted  by 
6 months  to  take  into  account  the  corresponding  shift  in  the 
solar  inclination  as  a function  of  latitude  (Widhopf,  et  al., 

197?).  Although  the  sample  data  in  Table  14  are  given  in  the 
model  for  the  month  of  August  in  both  hemispheres,  the  winter 
southern  latitudes  for  August  have  been  changed  in  this  table 


98 


Y PARAMETERS  FOR  CRUTZEN  MODEL  DURING  SUMMER  AND  WINTER  IN  NORTHERN 
HEMISPHERE*  (Source:  Hidalgo  and  Crutzen,  1977) 


Ed'S / mrireters  are  M /er  as  a function  of  season  (instead  of  nort^l  in  tfe  Crutzen  rwviel . 


back  to  the  corresponding  winter  northern  values  during  February. 
The  reason  for  this  latitudinal  shift  is  to  facilitate  a com- 
parison of  the  Widhopf  and  Crutzen  eddy  coefficient  data  for 
winter  in  the  Northern  Hemisphere.  This  shift  from  southern 
(August)  back  to  northern  (February)  latitudes  has  required  a 
repetition  for  winter  of  the  zero  latitude  column  during  sum- 
mer. Although  a more  consistent  comparison  of  the  data  in 
Tables  13  and  14  would  require  the  seasonal  averages  of  the 
Widhopf  data,  these  tables  suggest  rather  significant  differences 
between  these  two  sets  of  eddy  coefficients  for  the  lower  strat- 
osphere in  either  season. 

A third  parameterization  of  the  2-D  dynamics,  based  on  the 

Louis  meridional  circulation,  was  given  by  Danielsen  and  Louis 

(1977).  This  parameterization  utilized  the  observed  transports 

of  ozone  in  the  lower  stratosphere,  based  on  arguments  that 

(a)  ozone  may  be  considered  to  be  a passive  tracer  below  about 

21  km  (i.e.y  Section  2.1  and  Fig.  5),  and  (b)  ozone  is  the  best- 

documented  tracer  in  the  lower  stratosphere . The  calculations 

of  the  K , K , and  K coefficients  utilized  the  following 
y y y z z z 

assumptions:  (1)  the  horizontal  eddy  flux  of  ozone  is  propor- 

2 

tional  to  its  mean  flux,  and  (2)  the  product  [a]  K in  the 

«/  J 

vertical  flux  (Eq.  15)  must  be  a minimum  so  as  to  avoid  large 

values  of  [a],  which  would  give  vertical  fluxes  of  ozone  that 

are  greater  than  the  ozone  destruction  at  the  ground.  These 

two  assumptions  gave  two  relationships  between  K and  Kyz. 

The  third  coefficient,  K , was  obtained  by  solving  Eq.  (9) 

z z 

for  ozone  in  the  lower  stratosphere.  Table  15  illustrates 
results  from  these  calculations,  where  the  K values  are  given 

J “ 

in  terms  of  [a]  through  the  ratio  K /K  . These  results  are 

J " J J 

for  December-February ; i.e.,  for  winter  in  the  Northern  Hemi- 
sphere and  summer  in  the  Southern  Hemisphere.  However,  if  the 
eddy  coefficients  in  either  hemisphere  are  assumed  to  depend 
only  on  season  (as  in  the  Widhopf  model),  the  December-February 
summer  data  for  the  Southern  Hemisphere  also  becomes  applicable 
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for  the  summer  (June-August)  in  the  Northern  Hemisphere.  The 
validity  of  this  assumption  for  the  Danielsen-Louis  parameteri- 
zation may  be  somewhat  constrained  by  the  difference  in  the  ozone 
distributions  in  the  two  hemispheres  as  reflected  by  the  data 
in  Fig.  3.  Based  on  this  assumption,  a comparison  of  the  results 
in  Tables  13  and  15  indicates  rather  drastic  differences  for 
either  season  between  the  Crutzen  and  Danielsen-Louis  param- 
eterizations  of  the  eddy  coefficients,  both  of  which  utilized 
ozone  distributions. 

The  Danielsen-Louis  parameterization  of  the  eddy  transports 
was  also  tested  using  numerical  simulations  of  segments  of  time 
history  of  the  meridional  distributions  of  tracers  deposited 
in  the  lower  stratosphere  by  nuclear  detonations.  These  tests 
simulated  the  transports  of  Zr-95  from  a series  of  Chinese  nu- 
clear tests  (1968-1972),  and  a high-yield  1 9 6 8 French  test  as 
well  as  the  transport  of  W-I85  from  the  U.S.  Hardtack  test  at 
11°N  latitude  during  1958.  The  results  of  these  simulations 
indicated  that  use  of  the  Louis  meridional  circulation  (Fig.  17) 
and  the  eddy  coefficients  in  Table  15  gave  depletion  rates  of 
the  stratospheric  burdens  of  such  tracers  that  were  too  large 
compared  with  observations.  Significant  improvements  in  these 
simulations  were  obtained  by  decreasing  by  a factor  of  2 both 
the  air  velocities  in  Fig.  17  and  the  eddy  coefficients  in 
Table  15.  This  procedure  was  justified  by  citing  the  uncer- 
tainty of  a factor  of. about  2 in  the  derivation  of  the  Louis 
meridional  circulation. 

Since  Fig.  18  indicates  that  for  countergradient  eddy  trans- 
ports (i.e.,  [a]  >>  [6]),  [a]  = [w*R* ]/[ v*R* ] , which  from  Eqs. 

(13)  is  also  given  by  [a]  = K /K  , then  K /K  = [w#R* ]/[v*R»] 

y ^ y y y^yy 

Hence,  Eqs.  (11a, b)  yield  (Mahlman,  1975): 


[v»R»r 

[v*R*  ] ^51  + [w*R*  ] ILS  ’ 


(16) 
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TABLE  15.  SAMPLE  EDDY  PARAMETERS  DERIVED  FROM  POLEWARD  TRANSPORT  OF  OZONE* 
(Source:  Danielsen  and  Louis,  1977) 
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♦Data  is  applicable  to  December-February;  i.e.,  for  summer  in  the  Southern  Hemisphere  and 
winter  in  the  Northern  Hemisphere.  See  text  for  the  switch  in  summer  latitudes  from 
south  (December-February)  to  north  (June-August). 


and 


K = _Lw*H«]fc 

ZZ  [v»R»]  ijS_+  t w*R*  ] 

o y d 

Equations  (16)  and  (17)  also  identify  the  parameters  that  would 
have  to  be  deduced  from  observations  for  a reliable  2-P  param- 
eterization of  the  eddy  coefficients.  These  parameters  involve 
the  eddy  horizontal  and  vertical  mass  fluxes  as  well  as  the 
corresponding  gradients  of  the  mixing  ratio  as  a function  of 
altitude,  latitude,  and  month  or  season. 


Equations  ( 16 ) and  (17),  together  with  [a]  = [w*R* ]/[ v*R* ] , 

p 

yield  K = [a]cK  , a result  that  is  consistent  with  Eqs . (13) 

J ZZ  yy’ 

because  the  condition  [at]  = [w*R#  ]/[  v*R*  ] also  implies  from  Eqs. 
(1*0  that  [a*‘]/[a]  <<  1.  The  foregoing  result  can  be  put  in  the 
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[a]2K 


yy 


i. 
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which  is  a condition  that  can  be  used  to  test  the  internal  con- 
sistency of  the  parameterizations  in  Tables  13  through  15.  In- 
spection of  these  tables  at  either  constant  altitude  (rows)  or 
constant  latitude  (columns)  Indicates  that  this  consistency 

p 

condition  is  violated  (i.e.,  K^/Ca]  Kyy  ¥■  1)  whenever  there  is 
mainly  a rapid  change  of  a alone.  This  consistency  test  reveals 
the  following: 

(a)  The  Crutzen  2-D  parameterization  (Table  13)  is  con- 
sistent within  itself  only  at  middle  and  high  latl- 

p 

tudes  within  the  stratosphere , since  K /[a]‘  K >>  1 

zz  yy 

at  low  latitudes  (2Q°N  £ O s;  20°S)  in  both  the  lower 

stratosphere  and  upper  troposphere.  Furthermore, 

2 

K„„/[arKyy  increases  rapidly  above  unity  in  the  up- 
per troposphere  at  any  latitude. 

(b)  The  Widhopf  2-D  parameterization  (Table  1*0  is  con- 
sistent within  itself  only  at  middle  latitudes  within 

10*4 


■ U 


p 

the  stratosphere,  because  Kzz/[a]"K  > 1 at  high  (80°) 
and  low  (6  < 20°)  latitudes  in  the  lower  stratosphere. 
Again,  K /[a]  K increases  rapidly  above  unity  in  the 
upper  troposphere  in  the  20°  < 6 < 60°  latitudinal 
range . 

(c)  The  Danielsen-Louis  2-D  parameterization  is  consistent 
within  itself  only  at  middle  latitudes  (40°  £ 0 s 60°) 

p 

within  the  stratosphere  below  20  km,  because  K /[a]  K 

zz  yy 

is  otherwise  either  larger  (>  1)  or  much  larger  (>>1) 
than  unity . 

The  criterion  given  by  Eq . ( 1 8 ) indicates  that  the  Widhopf 

parameterization  is  the  most  consistent  at  low  latitudes  in  the 

2 2 

stratosphere  because  K /[a]  K > 1 instead  of  K /[a]  K >>  1 

z z yy  z z y y 

as  in  the  other  two  parameterizations . These  consistency  tests 
reveal,  then,  an  extra  degree  of  freedom  used  in  these  param- 
eterizations as  the  result  of  apparent  violations  of  the  inter- 
relationships among  the  eddy  parameters  as  given  by  Eqs . (13). 
Furthermore,  the  inconsistency  at  low  latitudes  becomes  impor- 
tant for  the  proper  interpretation  of  the  transport  of  engine 
effluents  from  the  Northern  to  the  Southern  Hemisphere;  i.e., 
for  the  Impact  of  high-altitude  flight  in  the  flight  corridors 
of  the  Northern  Hemisphere  on  the  ozone  column  in  the  Southern 
Hemisphere . 

It  is  important  to  indicate  that  the  eddy  coefficients 
(e.g..  Tables  13  through  15)  are  assumed  to  be  properties  of 
the  atmosphere  which  are  taken  to  be  the  same  in  the  calcula- 
tion of  the  eddy  transports  in  Eq.  (9)  for  each  of  the  chemical 
species . 

Equations  (9)  and  (11a,  b)  also  identify  model-dependent 
effects  in  the  Crutzen  and  Widhopf  2-D  photochemical  models, 
which  involve  uncertainties  in  the  magnitudes  of  the  photolysis 
rates  and  chemical  rate  coefficients  (e.g.,  Table  3)  as  well  as 
eddy  coefficients  (Tables  13  and  1^4).  These  model-dependent 
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effects  are  a consequence  of  the  basic  constraints  indicated 
previously  (Section  1)  and  they  have  been  responsible  for  the 
poor  forecast  record  of  aviation  effects  on  ozone  (Section  3). 
For  example,  since  Eqs.  (11a,  b)  with  x = show  that  the  eddy 
fluxes  in  Eq.  (9)  are  proportional  to  the  gradients  of  the  mix- 
ing ratio  of  the  i species,  it  becomes  clear  that  the  Howard- 
Evenson  rate  coefficient  in  the  revised  1977  chemistry  has  dis- 
torted the  previous  calculated  transports  of  mixing  ratio  by 
the  meriodional  circulation  and  eddy  motions  through  distortions 
of  the  meridional  distributions  of  the  mixing  ratios  . Hence, 
any  compensation  of  the  imprecise  knowledge  of  photolysis  and 
chemical  kinetic  rates  through  artibrary  adjustments  in  the 
parameterizations  of  the  eddy  transports  or  vice  versa  would 
not  necessarily  be  justified. 


4.2  EVALUATION  OF  THE  EXTENSION  OF  THE  DIFFUSION  HYPOTHESIS 

TO  THE  PERTURBED  ATMOSPHERE 

The  applicability  of  arbitrary  extensions  of  the  2-D  param- 
eterizations of  the  dynamics  for  the  reference  upper  troposphere 
and  lower  stratosphere  to  describe  the  mass  transports  of  engine 
effluents  from  an  assumed  large,  high-altitude  aircraft  traffic 
requires  considerations  of  (a)  the  degree  of  empiricism  in  the 
available  2-D  parameterizations  of  the  dynamics  as  revealed  by 
the  dynamic  meteorology  of  the  region,  and  (b)  the  mass  trans- 
ports of  passive  tracers  in  the  upper  troposphere  and  lower 
stratosphere  as  determined  from  3-D  solutions  of  the  equations 
of  atmospheric  motions  for  instantaneous  (e.g.,  nuclear  detona- 
tions) and  continuous  (i.e.,  aircraft-like)  sources  of  such 
material  in  the  upper  troposphere  and/or  lower  stratosphere. 

These  considerations  are  given  below. 

4.2.1  Dynamic  Meteorology  of  Eddy  Transports  in  the  Stratosphere 

The  foregoing  2-D  empirical  parameterizations  of  large- 
scale  eddy  transports  of  mass  in  the  reference  stratosphere 
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(Tables  13  through  15)  utilize  two  basic  assumptions:  (1)  the 
meridional  circulation  is  independent  of  the  eddy  motions,  and 
(2)  the  eddy  components  of  the  mixing  ratios  (R*)  are  propor- 
tional to  the  meridional  gradients  of  the  zonally  averaged  mix- 
ing ratio  (Eq.  10).  The  validity  of  these  assumptions  is  not 
supported  by  the  fundamental  primitive  equations. 

The  differential  equations  for  the  meridional  circulation 
and  eddy  motions  can  be  obtained  from  the  primitive  Equations 
(1)  through  (3)  by  the  following  steps:  (1)  splitting  the 
dependent  variables  (x  = u,  v,  w,  T)  into  their  zonal  mean  and 
eddy  components  (Eq.  8);  (2)  taking  the  zonal  average  of  the 
primitive  equations,  and  (3)  subtracting  the  zonally  averaged 
equations  from  each  respective  primitive  equation.  Using  the 
altitude  instead  of  the  pressure  as  a vertical  coordinate,  the 
results  are  given  in  Tables  16  and  17  for  the  zonal  mean  cir- 
culations and  eddy  motions,  respectively.  The  equations  for  the 
zonal  mean  circulations  consist  of  four  equations  for  the  four 
dependent  variables  [v],  [w],  [u],  and  [T],  while  those  for 
the  eddies  give  the  corresponding  components  v*,  w# , u*,  and  T*. 

The  results  in  Table  16  indicate  that  the  meridional  cir- 
culation is  driven  by  the  eddy  motions  through  the  terms  [Fx], 

[F  ],  and  [G],  whereas  those  in  Table  17  show  complex  couplings 
among  the  eddies  and  the  meridional  circulation.  Therefore, 
the  system  of  equations  in  Tables  16  and  17  implies  that  the 
mass  transports  by  the  eddies  in  Eq . (9)  are  not  as  assumed  by 
Eqs.  (11a,  b),  since  the  poleward  and  vertical  eddy  velocities 
are  coupled  in  a very  complex  way  not  only  with  the  meridional 
circulation  ([v],  [w]),  but  also  with  the  zonal  circulation 
([u])  and  other  eddy  variables  such  as  the  eddy  component  of 
the  net  heating  (Q*).  Furthermore,  since  tracer  material  in 
the  reference  stratosphere  is  imbedded  in  the  eddy  field,  the 
results  in  Table  17  imply  that  the  eddy  mixing  ratio  (R£)  can- 
not be  as  is  assumed  by  Eq.  10. 
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TABLE  16.  PRIMITIVE  EQUATIONS  FOR  MERIDIONAL  CIRCULATION  (Source:  e.g.,  Holton,  1975) 


The  symbol  h in  Eqs.  (a)  and  (b)  denotes  the  earth's  rotation,  and  [$]  in  Eq.  (b)  is  given 
aUi/Sz  = R[T]/H;  where  R is  the  gas  constant  for  dry  air  and  H the  scale  height  (5  RTs/g, 
where  Ts  is  a constant  stratospheric  mean  temperature).  T0  in  Eq.  (c)  is  a basic  state 
temperature  and  p = ps  exp  (-z/H)  In  Eq.  (d),  where  p5  is  the  density  at  sea  level. 


TABLE  17.  PRIMITIVE  EQUATIONS  FOR  EDDY  MOTIONS  (Source:  e.g.,  Holton,  1975) 


i-  { (v*T*  - [v*T*])cosGl  + -57  (P0(w*T*  - [w*T*])} 


It  should  be  noted  that  results  in  Table  16  have  been  used 
to  determine  the  meridional  circulation  shown  in  Fig.  17.  Con- 
siderations of  the  scales  of  the  meridional  velocities  indicate 
that  [v]  ~ 1.5  m/sec  and  [w]  ~ 1 cm/sec.  Because  of  the  small- 
ness of  the  vertical  velocity,  its  magnitude  cannot  be  deter- 
mined from  direct  observations.  Therefore,  the  procedure  used 
by  Louis  consisted  of  using  Eqs . (c)  and  (d)  in  Table  16,  to- 
gether with  observed  temperature  fields,  and  estimates  of  the 
eddy  fluxes  of  heat  in  the  term  [G],  and  the  net  heating  rate 
per  unit  mass  [Q].  The  uncertainty  in  these  results  depends, 
then,  on  the  parameterization  of  the  eddy  transports  of  sensible 
heat  as  well  as  the  radiation  calculations. 

The  meridional  circulation  as  well  as  the  eddy  variables 
in  Table  17  can  be  determined  from  solutions  of  the  (unaveraged) 
primitive  equations  themselves,  which  can  be  solved  either  by 
finite  difference  or  spectral  numerical  techniques.  Figure  19 
shows  solutions  of  such  system  of  equations  based  on  the  latter 
technique  and  the  assumption  of  geostrophic  equilibrium  (Cunnold 
et  al.,  1975).  Although  the  analytical  results  in  Fig.  19 
resemble  those  in  Fig.  17,  which  were  obtained  from  observations, 
this  agreement  is  not  surprising;  because  the  model  specified 
the  radiative  equilibrium  temperature  distribution  and  the  static 
stability  instead  of  using  only  the  flux  of  solar  radiation 
across  the  outer  boundary  of  the  atmosphere. 

The  splitting  of  the  atmospheric  motions  into  average  cir- 
culations and  eddy  components  (Eq.  8),  as  used  in  the  results 
of  Tables  16  and  17,  turns  out  to  be  more  important  in  the 
stratosphere  than  in  the  troposphere.  A main  reason  for  this 
condition,  of  course,  is  the  absence  in  the  stratosphere  of  the 
topographical  and  thermal  forcings  of  the  troposphere.  For 
example,  it  has  been  shown  that  the  annual  energy  budget  of  the 
large-scale  eddy  transports  in  the  lower  stratosphere  is  driven 
primarily  by  planetary  waves  propagating  upwards  from  the 
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troposphere  during  winter  (Dopplick,  1971).  This  mechanism  can, 
surprisingly,  be  described  in  terms  of  linear  wave  mechanics 
and  requires  consideration  of  the  effects  of  the  meridional  cir- 
culation on  the  zonal  flow  [u].  The  Coriolis  torque  exerted  by 
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FIGURE  19.  Simulation  of  meridional  circulation  by  a three-dimensional 
model  of  the  troposphere-stratosphere-mesosphere  (0  < z < 

70  km)  utilizing  the  condition  of  geostrophic  equilibrium  and 
the  spectral  representation  of  the  dependent  variables.  Note 
the  simulation  of  the  three  cells  in  the  troposphere;  i.e., 
the  direct  Hadley  cell  at  low  latitudes,  the  indirect  Ferrel 
cell  at  middle  latitudes  and  the  polar  cell.  Note  also  the 
complex  downward  stratospheric  motions  at  the  tropopause  gaps 
between  the  Hadley  and  Ferrel  cells  at  about  30°  latitude 
during  winter  and  summer.  Corresponding  eddy  motions  are 
not  shown.  (Source:  Cunnold  et  al.,  1975) 
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the  meridional  flow  (Fig.  19)  generates  stratospheric  mean  zonal 
winds  ([u])  in  the  westward  direction  during  summer  and  in  the 
eastward  direction  during  winter.  This  change  in  the  sign  of 
the  stratospheric  zonal  wind  allows,  then,  vertical  propagation 
of  tropospheric  planetary  waves  to  the  stratosphere  during 
winter.  This  mechanism  was  demonstrated  by  Charney  and  Drazin 
(1961)  who  showed  that  quasl-geostrophic  waves  generated  in  the 
troposphere  cati  propagate  vertically  into  the  stratosphere  when- 
ever their  phase  speeds  are  (a)  westward  relative  to  the  [u] 
wind  and  (b)  smaller  than  a critical  velocity  which  decreases 
with  increasing  wavelength. 

4.2.2  MassTransports  of  Inert  Material  in  the  Reference 
and  Perturbed  Atmospheres 

The  only  available  procedure  for  assessing  the  applicability 
of  the  mixing  length  hypothesis  to  the  perturbed  upper  tropo- 
sphere and  lower  stratosphere  is  to  solve  numerically  the  primi- 
tive equations  for  injections  of  passive  tracers  in  the  upper 
troposphere  and/or  lower  stratosphere.  Although  chemically 
inert  and  without  absorption  or  emission  radiative  character- 
istics, such  tracers  must  be  assumed  to  be  water  soluble  so  as 
to  stimulate  their  heterogenous  removal  in  the  lower  tropospheric 
sinks.  To  determine  the  validity  of  the  mixing  length  hypo- 
thesis in  the  perturbed  upper  troposphere  and  lower  stratosphere, 
the  injections  of  tracer  materials  must  consider  two  cases,  as 
discussed  in  the  following  paragraphs. 

The  first  case  would  involve  instantaneous  injections  at  a 
point  in  the  lower  stratosphere  to  simulate  the  dispersion  of 
debris  from  nuclear  detonations  in  the  reference  upper  tropo- 
sphere and  lower  stratosphere;  i.e.,  by  considering  solutions 
of  Eq . (5)  for  the  conservation  of  the  passive  tracer  mixing 
ratio,  dR/dt  = 0,  subject  only  to  the  removal  of  the  tracer 
material  in  the  tropospheric  sinks.  These  solutions  must  con- 
sider an  initial,  local  distribution  of  tracer  and  provide  the 
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subsequent  time  evolutions  of  the  mixing  ratio  in  the  strato- 
sphere and  troposphere.  More  important,  these  numerical  solu- 
tions can  provide  insight  into  the  relative  roles  of  the  merid- 
ional circulation  and  eddy  motions  in  the  overall  transport  of 
tracer  mixing  ratio  in  the  stratosphere.  Another  result  would 
include  an  evaluation  of  the  validity  of  the  mixing  length 
hypothesis  itself  anywhere  in  the  reference  upper  troposphere 
and  stratosphere.  This  latter  result  would  be  derived  from 
(a)  comparisons  of  the  simulated  and  observed  dispersions,  as 
a function  of  time,  of  tracer  material  from  nuclear  detonations;  f 

and  (b)  comparisons  of  the  eddy  parameters  in  Tables  13  through 
15,  with  those  derived  from  the  eddy  transports  and  gradients 
of  mixing  ratio  obtained  from  the  primitive  equations  and  the 
assumption  that  Eqs . (11a,  b)  are  indeed  valid. 

The  second  case  would  require  continuous  injections  at  a 
typical  location  of  the  flight  corridors  to  simulate  the  dis- 
persion of  aircraft  engine  effluents  (NO  and  H^O)  in  the  per- 
turbed upper  troposphere  and  stratosphere;  i.e.,  by  considering 
solutions  of  Eq.  (5)  for  the  tracer  mixing  ratio,  dR/dt  = S(A, 

0,  z),  subject  to  (a)  a continuous  source  (S)  in  the  typical 
flight  corridor  as  a function  of  longitude,  latitude,  and  alti- 
tude; and  (b)  the  removal  of  tracer  material  in  the  tropospheric 
sinks.  Because  of  the  consideration  of  a passive  tracer,  the 
numerical  results  are  applicable  to  the  conservation  of  the  air- 
craft NC>x  and  species  during  their  dispersion  in  the  atmos- 

phere, i.e.,  to  NOy  (N0x  + HNO^)  as  well  as  to  exclusion  of  H^O 
from  mechanism  Mllab  in  Fig.  2 for  H^O.  These  solutions  must 
consider  numerical  time  integrations  for  a time  long  enough  (sev- 
eral years)  to  approach  equilibrium  conditions  (i.e.,  3R^/3t  -*•  0 
for  a given  month  or  season)  between  the  continuous  source  and 
the  tropospheric  sinks  in  the  Northern  and  Southern  Hemispheres. 

These  solutions  can  again  provide  the  eddy  parameters  based  on 
results  from  the  primitive  equations  together  with  the  mixing 
length  hypothesis  [Eqs.  (11a,  b)]  for  the  perturbed  upper  tropo- 
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sphere  and  stratosphere.  These  results  would  thus  establish  the 
applicability  of  the  eddy  parameters  for  the  reference  atmosphere 
when  extended  to  the  perturbed  conditions. 

Numerical  solutions  of  the  above  nature  have  been  carried 
out  at  the  Geophysical  Fluid  Dynamics  Laboratory  (GFDL)  by 
Mahlman  ( 197 3 > 1975)  and  Mahlman  and  Moxim  (1978).  The  basic 
approach  consisted  of  using  available  GCM  solutions  for  the 
primitive  equations  (1  through  4)  as  inputs  to  the  dynamical 
terms  in  the  expansion  of  dR/dt  (Eqs.  6 and  7).  The  numerical 
solutions  of  the  passive  tracer  mixing  ratio  are  then  based  on 
the  use  of  the  3-D  dynamics  for  one  annual  cycle  in  Eqs.  (6)  and 
(7);  i.e.,  the  tracer  integrations  proceeded  beyond  a year  by 
recycling  the  same  input  data  each  year.  Although  this  con- 
straint stemmed  from  otherwise  prohibitive  investments  in  com- 
puter time,  it  should  be  noted  that  this  procedure  did  keep 
all  the  terms  in  Tables  16  and  17  unchanged  from  one  annual 
cycle  to  the  next.  The  result  is,  then,  a tracer  behavior 
during  several  annual  cycles  that  is  solely  due  to  the  tracer 
Itself . 

The  tracer  solutions  provide  the  3-D  distributions  of  the 
tracer  mixing  ratio  as  a function  of  time,  but  they  can  be 
reduced  to  the  2-D  form  of  Eq.  (9)  by  performing  their  zonal 
averages.  These  operations  yield,  then,  each  of  the  four 
dynamic  terms  in  the  left  side  of  Eq.  (9),  which  may  be  written 
for  the  instantaneous  injection  case  as: 

• <19> 

where  M„  and  denote  the  terms  related  to  the  horizontal  and 
vertical  transports  of  mixing  ratio  by  the  meridional  circula- 
tion, E^  and  Ey  the  corresponding  eddy  terms,  and  ET  those  due 
to  the  sinks  as  well  as  small-scale  (subgrid)  diffusion  and 
filling  (i.e.,  a correction  to  the  advective  terms  whenever 


truncation  errors  lead  to  a local  negative  mixing  ratio).  Hence, 
it  becomes  possible  to  (a)  identify  the  relative  roles  of  each  of 
the  above  four  dynamic  terms  on  the  local  time  evolution  of  the 
mixing  ratio  during  any  season  as  a function  of  time  (years)  or 
configurations  of  the  tracer,  and  (b)  calculate  the  correspond- 
ing eddy  parameters  [a],  K , and  K from  the  ratio  of  the  eddy 

yy  z z 

fluxes  as  well  as  the  use  of  Eqs.  (16)  and  (17). 

Since  the  tracer  equation  may  be  made  dimensionless  by 
dividing  each  term  of  the  equation  by  a constant  reference 
value  of  the  mixing  ratio,  the  numerical  results  may  be  scalable 
to  any  other  amplitude  or  peak  value  of  the  initial  tracer  dis- 
tribution. Hence,  the  mixing  ratio  can  be  expressed  in  arbi- 
trary units,  denoted  as  R units.  Likewise,  the  rate  of  change 
of  mixing  ratio  can  be  based  on  a time  characteristic  of  the 
zonal  averaging  (e.g.,  a day);  hence  the  rate  of  change  of  the 
mixing  ratio  can  be  expressed,  in  general,  in  terms  of  R-units/ 
day. 

The  numerical  experiments  for  the  instantaneous  injection 
simulated  a nuclear  detonation  in  the  lower  stratosphere  at 
about  mid-latitudes  in  the  Northern  Hemisphere.  The  initial 
tracer  distribution  started  in  January  1 and  was  centered  at 
65  mbar  (approximately  19  km)  and  36°N.  The  initial  distribu- 
tion had  a peak  value  of  about  0.31  R units  at  65  mbar  and  36°N, 
a value  that  decayed  exponentially  with  either  increasing  or  de- 
creasing altitude  as  well  as  latitude.  The  numerical  integra- 
tions used  a time  step  of  24  minutes  for  integrations  of  five 
annual  cycles.  The  3-D  results  were  then  averaged  zonally  to 
obtain  the  equivalent  2-D  values,  and  monthly  to  determine  the 
time  evolutions  of  the  tracer  mixing  ratio. 

Figure  20  shows  the  meridional  distribution  of  tracer 
mixing  ratio  for  January  or  winter  (top)  and  July  or  summer 
(bottom)  during  the  first  (left)  and  third  (right)  years  of 
integrations.  The  vertical  shape  of  the  mixing  ratio  for  Janu- 
ary of  the  first  year  near  30°N  reflects  that  of  the  initial 
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FIGURE  20.  Simulation  of  8-lines  (Fig.  18)  by  a three-dimensional  passive  tracer  model  utilizing 
finite  differences  of  the  primitive  equations  for  an  instantaneous  injection  of  mate- 
rial at  36°N  and  65  mbar  during  January.  The  dashed  closed  contours  denote  the  loca- 
tion of  the  peak  initial  distribution.  (Source:  Mahlman  and  Moxim,  1978) 
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profile,  which  had  mixing  ratio  values  at  36°N  equal  to  one- 
third  the  6 5 mbar  peak  value  at  both  10  mbar  and  190  mbar.  The 
peak  value  of  100  (10-^  R units)  for  January'  of  the  first  year 
shows  that  the  mixing  ratio  there  has  decayed  by  more  than  a 
factor  of  about  30  from  its  initial  peak  value  (0.31  R)  during 
the  first  month  after  the  assumed  detonation.  The  July  contours 
for  the  first  year  indicate  that  this  decay  of  mixing  ratio  has 
proceeded  by  another  factor  of  five  after  the  first  half  year, 
when  the  formation  of  the  [g]  lines  (Fig.  18)  becomes  evident. 

The  January  and  July  meridional  distributions  of  mixing  ratio 
during  the  third  year  are  typical  of  those  for  the  preceding 
and  subsequent  years.  Figure  21  shows  a cross  plot  of  results 
for  the  end  of  the  fourth  year  similar  to  those  in  Fig.  20,  for 
constant  pressure  levels  during  winter  (December) . Figure  21 
indicates  that  the  mixing  ratio  at  the  location  of  the  peak 
initial  value  of  mixing  ratio  has  dropped  from  the  value  at  the 
beginning  of  the  third  year  by  a factor  of  about  three  during 
the  next  2 years. 

Figure  22  shows  the  dynamic  interactions  (Eq.  19)  at  the 
65  mbar  injection  level  corresponding  to  the  meridional  distri- 
butions in  Fig.  20.  The  3-D  results  shown  in  Fig.  22  indicate 
considerable  compensation  between  the  mean  meridional  circula- 
tion and  the  eddy  fluxes.  As  pointed  out  by  Holton  (1978), 
these  compensating  effects  are  not  accidental,  because  Andrews 
and  McIntyre  (1976)  have  already  shown  that  the  mean  meridional 
circulation  can  be  divided  into  two  parts:  a portion  due  to 
the  diabatic  heat  sources  and  a portion  driven  by  the  eddies. 

For  stationary  eddies,  the  latter  exactly  cancels  the  eddy 
fluxes.  Most  of  the  mean  meridional  circulation  is  thus  driven 
by  the  eddies.  All  of  this  would  be  much  less  confusing  if  a 
Lagrangian  instead  of  an  Eulerian  point  of  view  were  to  be  used. 
Thus,  it  is  to  be  expected  that  eventually  McIntyre’s  "Lagrangian 
mean  formulation"  will  replace  the  present  awkward  parameteri- 
zations  in  2-D  models.  As  pointed  out  again  by  Holton  (1978), 
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ZONAL  MEAN  TRACER  MIXING  RATIO,  (10  R units) 
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FIGURE  21.  Stratospheric  tracer  mixing  ratio  as  a function  of  latitude, 
altitude  and  season  for  winter  (December)  of  the  fourth  year 
after  an  assumed  instantaneous  release  of  tracer  material  at 
36°N  and  65  mbar.  (Source:  Mahlman  and  Moxim,  1978) 
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FIGURE  22.  Interactions  at  the  Injection  level  (65  mbar)  among  the  eddy 
and  meridional  motions  in  the  dispersion  of  tracer  material 
(Eq.  19)  in  the  stratosphere  during  winter  (top)  and  summer 
(bottom)  (Source:  Mahlman  and  Moxlm,  1978) 
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these  considerations  are  limited  to  the  stratosphere,  where  2-D 
models  are  expected  to  be  more  useful  than  in  the  troposphere; 
i.e.,  the  2-D  model  approach  is  suspect  in  the  troposphere  due 
to  the  presence  there  of  very  large  zonal  asymmetries  in  terms 
of  surface  sources  and  sinks. 

It  is  of  interest  to  note  that  the  dynamic  interactions 
shown  in  Pig.  22  are  different  for  the  two  annual  cycles  Indi- 
cated there,  even  for  the  July  interactions  when  the  correspond- 
ing isopleths  of  the  mixing  ratio  are  more  similar  to  each  other 
(Pig.  20).  Just  as  the  difference  in  these  dynamic  interactions 
between  January  and  July  of  the  third  year  with  somewhat  similar 
[f$]-profiles  suggest  that  the  eddy  parameters  (Tables  13  through 
15)  would  be  a function  of  season,  the  difference  in  these  dyna- 
mic interactions  for  July  of  the  first  and  third  years  suggests 
that  the  eddy  parameters  would  also  be  a function  of  the  configu- 
ration of  the  tracer  during  a given  season.  Since  these  results 
are  based  on  fixed  dynamics  (I.e.,  no  variability)  during  each 
annual  cycle,  they  are  caused  only  by  the  transports  of  the  tracer 
(Tables  16  and  17).  Because  of  the  differences  between  the  tracer 
configurations  for  nuclear  debris  and  ozone  distributions  in  the 
stratosphere  during  a given  month  of  successive  annual  cycles,  the 
results  in  Pigs.  20  and  22  already  suggest  that  (a)  the  test  of 
the  Danielsen-Louis  parameterization  of  the  eddy  transports  for 
the  reference  atmosphere  baspd  on  the  dispersion  of  nuclear  debris 
may  be  inconsistent  with  their  derivation  from  the  ozone  distribu- 
tions, (b)  parameterization  of  the  eddy  transports  for  the  refer- 
ence atmosphere  based  on  the  dispersion  of  nuclear  debris  may  have 
lower  accuracy  than  those  utilizing  observed  ozone  as  an  inert 
tracer  in  the  lower  stratosphere,  as  proposed  by  Danielsen  and 
Louis  (1977),  and  (c)  the  assumption  of  the  same  eddy  coefficient 
for  every  one  of  the  chemical  species  with  different  configura- 
tions may  be  invalid.  Furthermore,  because  of  the  expected  dif- 
ferences between  the  tracer  configurations  of  aircraft  effluents 
and  those  of  either  ozone  or  nuclear  debris  in  the  reference 


atmosphere,  the  available  parameterizations  of  the  eddy  trans- 
ports for  the  reference  atmosphere  may  not  be  accurate  for  their 
extensions  to  the  dispersion  of  engine  effluents  in  the  perturbed 
stratosphere . 

The  numerical  experiments  for  the  continuous  injections 
simulate  SST  traffic  at  middle  latitudes  in  the  lower  strato- 
sphere of  the  Northern  Hemisphere.  The  simulation  of  the  air- 
craft traffic  was  assumed  to  be  given  by  a source,  S,  of  mixing 
ratio  given  by 


/ 0 -48°N  ^ 2 

t z-16.5  Y 

L \ 12  / J 

exp 

\~~rrs-  ) J 

(20) 


i.e.,  by  assuming  that  the  peak  source  (SQ)  from  aircraft  traf- 
fic occurs  at  0Q  = 48°N  and  z = 16.5  km.  Thus,  most  of  the  SST 
tracer  material  was  deposited  at  the  110  mbar  level  and  in  the 
altitude  range  14  < z < 18.5  km.  At  either  end  of  this  altitude 
range,  the  tracer  was  assumed  to  drop  to  zero.  The  longitudinal 
(A)  dependence  of  the  SST  source  was  assumed  to  be  unity  over 
the  North  Atlantic  (0°W  < A < 70°W)  and  East  Pacific  (130°W  < 

A < 140°E),  and  zero  otherwise. 

The  initial  conditions  for  this  continuous  source  experi- 
ment utilized  available  results  for  the  buildup  of  tracer  mix- 
ing ratio  as  obtained  from  the  following  assumptions  (Mahlman, 
1973):  (a)  the  SST  problem  could  be  represented,  approximately, 

by  a superposition  of  individual  tracer  time  histories  of  dif- 
ferent age;  which  are  determined  from  a long  series  of  single 
sources  pulsing  with  a frequency  of  once  every  few  hours,  and 
(b)  the  time  history  of  mixing  ratio  from  each  pulse  is  inde- 
pendent of  that  of  the  others,  an  assumption  that  was  based  on 
the  fact  that  the  ultimate  global-scale  distribution  of  a given 
tracer  release  is  nearly  independent  of  the  initial  (January) 
meteorological  conditions  of  the  single  source  experiment  (e.g.. 
Pig.  20).  It  is  important  to  indicate  that  the  main  motivation 
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for  using  these  initial  conditions  is  avoidance  of  prohibitive 
investments  in  computer  time  when  reaching  the  equilibrium 
tracer  distributions  given  by  the  eventual  balance  between  the 
continuous  injections  and  the  tropospheric  sinks.  Since  this 
superposition  procedure  had  been  carried  out  for  the  first  4 
years  of  assumed  SST  operation,  the  subsequent  numerical  experi- 
ments using  a continuous  source  (e.g.,  Eq.  20)  extends  currently 
to  3 additional  years  (fifth  to  seventh)  of  assumed  SST  operations, 

Figures  23  and  2b  indicate  results  for  the  enhancement  in 

the  tracer  mixing  ratio  during  winter  (December)  and  summer 

(August),  respectively,  from  SST-like  sources  as  given  by  Eq . 

20.  Since  the  tracer  equation  may  again  be  made  dimensionless 

by  dividing  it  by  the  strength  of  the  source  term,  these  results 

are  expressed  in  an  arbitrary  but  convenient  unit  that  allows 

their  scaling  to  any  specified  source  intensity.  The  units  in 

Figs.  23  and  2b  represent  the  enhancement  of  mixing  ratio  from 

one  day’s  emissions  over  a reference  value  given  by  a uniform 

dilution  of  the  daily  tracer  injection  into  the  whole  atmosphere. 

a 

For  example,  for  an  annual  emission  rate  of  10  kg  of  tracer 

material  per  year,  one  day’s  emission  would  be  2.7  x 10^  kg  per 

day.  The  uniform  dilution  of  this  daily  emission  into  the  mass 

of  the  atmosphere  (5.1  x 10  g of  air)  gives  a reference  value 
-13 

of  5.4  x 10  per  day.  Since  Fig.  2b  gives  an  enhancement  of 
about  3800  over  the  reference  value  at  40°N  and  65  mbar  (19*0  km) 
during  summer,  the  predicted  enhancement  in  mixing  ratio  at  this 
location  would  be  2.0  ppbm.  Because  the  enhancements  in  Fig.  2b 
are  approaching  equilibrium  conditions  between  the  continuous 
source  and  the  tropospheric  sinks,  this  enhancement  in  mixing 
ratio  at  40°N  and  65  mbar  is  also  approaching  the  maximum  per- 
turbation in  the  tracer  mixing  ratio  for  a given  intensity  of 
the  continuous  source.  The  results  in  Figs.  21,  23,  and  2b  in- 
dicate the  following: 

• There  is  a significant  difference  in  the  tracer 

configurations  of  the  reference  and  perturbed  lower 


ENHANCEMENT  OF  ZONAL-MEAN  TRACER  MIXING  RATIO,  [R] 


NORTH  LATITUDE  SOUTH 


in  in 


FI6URE  23.  Enhancement  of  tracer  mixing  ratio  as  a function  of  latitude 
and  altitude  during  winter  (December)  of  the  fifth  year  for 
a simulated  SST -like  continuous  release  of  tracer  material  as 
per  Eq.  (20).  Note  the  tendency  towards  a uniform  distribu- 
tion of  tracer  mixing  ratio  between  38  mbar  (standard  altitude 
of  22.4  km)  and  10  mbar  (31.2  km)  in  the  Northern  Hemisphere. 
(Source:  Mahlman  and  Moxim,  1978) 
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ENHANCEMENT  OF  ZONAL-MEAN  TRACER  MIXING  RATIO,  [R] 
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FIGURE  24.  Enhancement  of  tracer  mixing  ratio  as  a function  of  latitude 
and  altitude  during  summer  (August)  of  the  sixth  year  for  a 
simulated  SST-like  continuous  release  of  tracer  material  as 
per  Eq.  (20).  (Source:  Mahlman  and  Moxim,  1978) 
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stratospheres,  for  example,  during  winter  in  the  North- 
ern Hemisphere.  Figure  21  for  the  reference  stratosphere 
shows  significant  vertical  gradients  in  the  tracer  mixing 
ratio  at  40°N  above  15.6  km,  whereas  Fig.  23  for  the  per- 
turbed stratosphere  indicates  corresponding  negligible 
vertical  gradients  and  a nearly  uniform  mixing  ratio  at 
every  northern  latitude  above  22.4  km.  The  significance 
of  these  results  is  that  the  mixing  length  hypothesis 
would  then  predict  negligible  vertical  eddy  motions  in 
a uniformly  mixed  perturbed  stratosphere  during  winter 
[i.e.,  a = B = = 0 in  Fig.  18  and  w*  = 0 in  Eq.  (12)'], 

if  it  were  used  with  the  perturbed  transient  tracer  con- 
figuration derived  from  the  primitive  equations . Con- 
versely, the  eddy  coefficients  derived  from  the  mixing 
length  hypothesis  from  the  tracer  configuration  in  the 
reference  stratosphere  ( i.e .,  [a]  ^ 0 , [B]  ^ 0 , [a]  >> 
[S])  could  not  predict  the  transient , uniformly  mixed 
conditions  of  the  perturbed  stratosphere  obtained  from 
the  primitive  equations . 

• Figures  23  and  24  for  the  perturbed  winter  and  summer 
seasons  indicate  that  the  tracer  emitted  at  the  contin- 
uous source  propagates  to  the  ozone  production  region 
(Ml,  Table  1)  at  the  higher  altitudes  and  lower  latitudes 
in  the  Northern  Hemisphere.  Hence,  the  upward  vertical 
transports  of  tracer  mixing  ratio  from  the  continuous 
source  pressure  level  to  the  middle  stratosphere  (20  < 
z < 30  km)  are  also  important  in  the  assessment  of  avi- 
ation effects  on  the  ozone  column. 

The  significant  difference  in  the  tracer  configurations  of 
Figs.  21  and  23  can  further  be  illustrated  by  examining  the  rela- 
tive roles  of  the  mass  transports  (Eq.  19)  by  the  meridional  cir- 
culation and  eddy  motions  in  the  reference  and  perturbed  strato- 
spheres. These  interactions  are  shown  in  Figs.  25  and  26,  which 
imply  a rather  drastic  difference  in  the  relative  roles  of  the 
meridional  circulation  and  eddy  motions  in  the  transports  of 
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FIGURE  25.  Eddy  transports  and  the  local  rate  of  change  of  mixing  ratio 
(Eq.  19)  at  65  mbar  during  winter  (December)  of  the  fourth 
year  after  an  instantaneous  release  of  passive  tracer  material 
(Source:  Mahlman  and  Moxim,  1978) 


mixing  ratio  at  65  mbar  in  the  reference  (Ftg.  25)  and  perturbed 
(Fig.  26)  stratospheres  even  though  the  continuous  source  itself 
is  at  a lower  pressure  level  110  mbar).  The  significance  of 
these  results,  utilizing  the  same  dynamics  (Tables  16,  17)  for 
the  reference  and  perturbed  atmospheres,  is,  again,  that  exten- 
sions of  empirical  parameterizations  of  eddy  mass  transports  in 
the  reference  stratosphere  to  describe  the  eddy  transports  of  en- 
gine effluents  in  the  perturbed  stratosphere  cannot  be  justified 
from  the  conservation  equations  for  the  stratospheric  motions. 

Returning  to  the  basic  considerations  concerning  the  hypo- 
thesis of  large-scale  diffusion  in  the  stratosphere  (Section 
4.1.2),  it  now  becomes  important  to  emphasise  that  the  Newell 
(1964)  conceptual  explanation  of  countergradient  eddy  fluxes 
of  ozone  or  heat  could  not  be  based  on  actual  observations . 

The  reasons  for  this  situation  are  (a)  the  physical  impossibility 
of  observing  the  slope  (a)  of  the  mixing  length  path  (Fig.  18) 
due  to  the  fact  that  the  eddy  transports  are  not  actually  de- 
coupled from  the  meridional  circulation  (Tables  16  and  17),  and 
(b)  the  smallness  of  the  slope  [a]  as  indicated,  for  example,  by 
the  Crutzen  and  Widhopf  models  themselves  (Tables  13  and  14). 
Furthermore,  since  the  numerical  solutions  of  the  primitive 
equations  for  the  instantaneous  and  continuous  injections  yield 
the  vertical  and  horizontal  eddy  transports  of  tracer  mixing 
ratio,  it  becomes  possible  to  calculate  a as  stipulated  in  Fig.  18; 
i.e.,  [a]  = [w*R*]/[v#R*] . These  results  are  illustrated  in  Table 
18  for  winter,  corresponding  to  the  conditions  in  Figs.  21  and  23, 
and  it  includes  those  for  summer  (June)  for  the  SST  case.  A com- 
parison of  the  winter  [a]  values  for  the  instantaneous  injection 
in  Table  18  with  the  corresponding  Crutzen  and  Widhopf  values  in 
Tables  13  and  14  indicates  that  the  absolute  magnitudes  of  [a] 
in  the  lower  stratosphere  at  the  middle  northern  latitudes,  for 
example,  would  be  larger  for  the  solutions  of  the  primitive  equa- 
tions. This  discrepancy  in  the  ratio  of  vertical  to  horizontal 
eddy  transports  leads,  then,  to  drastic  corresponding  differences 

in  the  K__  and  K„„  coefficients,  which  can  be  determined  from 
zz  yy 
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FIGURE  26.  Eddy  transports  and  local  rate  of  change  of  mixing  ratio 
(Eq.  19)  at  65  mbar  during  winter  (December)  after  5 years 
of  an  SST-like  continuous  release  of  passive  tracer  material. 
(Source:  Mahlman  and  Moxim,  1978) 
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TABLE  i8.  SLOPE  OF  MIXING  LENGTH  PATH  AS  GIVEN  BY  PRIMITIVE 
EQUATIONS,  [a],  (10-4)* 

(Source:  Mahlman,  1978) 


For  instantaneous  injection  at  38  mbar  and  72°N,  for  example,  [a] 


-*7 
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Eqs . (16)  and  either  (17)  or  (18),  using  results  from  the  numeri- 
cal solutions  of  the  primitive  equations  for  the  eddy  fluxes  and 
gradients  of  the  mixing  ratio.  The  results  are  as  given  in  Tables 

19  and  20,  which  indicate  that  K , and  K for  the  reference 

yy  zz 

stratosphere  would  have  to  be  either  positive  or  negative.  As 

indicated  by  Eq.  (18),  the  change  in  the  signs  of  K and  K 

zz  yy 

must  take  place  simultaneously  at  a given  location.  Furthermore, 
Eqs.  (16)  and  (17)  show  that  this  change  of  sign  must  come  from 
the  denominator.  Hence,  every  time  there  is  a change  of  sign  in 
Kyy  and  Kzz,  the  denominator  must  go  through  zero;  i.e.,  a condi- 
tion that  implies  extremely  large  .local  values  for  K and  K as 

zz  yy 

both  become  infinity. 

The  results  in  Tables  18-20  include  the  eddy  coefficients 
derived  from  the  primitive  equations  for  the  perturbed  winter 
and  summer  stratospheres.  The  results  in  Table  18  indicate  that 
the  slopes  of  the  mixing  length  path  for  winter  conditions  in 
the  reference  and  perturbed  stratospheres  are  not  identical, 
for  example,  in  the  lower  stratosphere,  at  middle  northern  lati- 
tudes, which  is  an  assumption  in  the  Crutzen  and  Widhopf  2-D 
photochemical  models.  Similar  considerations  apply  to  the  K 

Z c» 

and  K , eddy  coefficients  as  exhibited  in  Tables  19  and  20. 

yy 

Finally,  it  must  be  recognized  that  the  Mahlman-Moxim 
numerical  solutions  of  the  primitive  equations  may  not  neces- 
sarily be  an  exact  simulation  of  the  dispersion  of  engine 
effluents  in  the  actual  perturbed  stratosphere.  However,  it 
is  also  important  to  emphasize  the  following: 

• The  evaluation  of  the  extension  of  the  mixing  length 
hypothesis  to  the  perturbed  stratosphere  is  based  on 
a relative  instead  of  absolute  basis  (e.g.,  Fig.  25 
vs.  Fig.  26). 

• The  primitive  solutions  do  duplicate,  first,  the  ob- 
served behavior  of  the  slope  [6]  (Figs.  18  and  20)  for 
the  tracer  mixing  ratio  at  middle  northern  latitudes 
in  the  reference  stratosphere;  and  second,  the  charac.- 
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For  instantaneous  injection  at  38  mbar  and  72°N,  for  exairple,  K 


teristics  of  the  other  parameters  (e.g.,  potential  vor- 
ticlty)  which  are  strongly  correlated  with  the  ozone 
distributions  in  the  stratosphere  (Mahlman  and  Moxim, 
1978) . 

• The  configuration  of  tracer  mixing  ratio  in  the 
reference  stratosphere  is  in  qualitative  agreement 
with  the  rather  rough  observations  of  the  dispersion 
of  nuclear  debris  in  the  lower  stratosphere  (Tele- 
gadas,  197^;  Mahlman  and  Moxim,  1978). 

• The  empirical  parameterization  of  eddy  coefficients 
based  on  the  Reed  and  German  (1969)  use  of  the  Frandtl 
mixing  length  concept  for  the  eddy  mass  transports  in 
the  reference  atmosphere  is  not,  by  itself,  a critical 
issue.*  It  should  be  noted  that  such  parameter i za- 
t-ions  could  be  based  on  other  hypotheses  such  as  the 
G.  I.  Taylor  (193?)  formulation  for  microscale  tur- 
bulent transports  of  vorticity  and  heat  (instead  of 
momentum,  as  in  the  mixing  length  hypothesis)  or  even 
the  arbitrary  modifications  of  Eq . (10)  for  R* . 

Since  the  empirical  parameterizat ions  of  the  (unknown) 
eddy  coefficients  in  the  lower  stratosphere  do  not 
require  a unique  mathematical  model,  their  validity 
must  rest  on  the  accuracy  of  t lie  appropriate  observa- 
tions as  well  as  on  the  internal  consistency  of  their 

V 

In  fairness  to  the  memory  of  a late  scientist,  and  with 
permission  of  G.  D.  Robinson,  it  becomes  appropriate  to 
quote  his  review  comment  in  regard  to  the  application 
of  the  Prandtl  mixing-length  concept  to  the  mass  trans- 
ports in  the  atmosphere:  "Since  I once  met  Frandtl,  my 
sympathies  go  out  to  him,  revolving,  as  he  must  be,  in 
his  grave.  What  he  said  was — if  an  entity  is  conserved 
in  a fluctuating  transport  pnoce ss,  then  the  averaged  flux 
in  a given  direction  is  proportional  to  the  component  in 
that  direction  of  the  gradient  of  the  averaged  concent  ra- 
tion. Time  averages  are  implied,  over  times  long  compared 
with  the  time  scale  of  fluctuations.  There  are  not  many 
entitles  or  locations  In  the  atmosphere  which  meet  these 
conditions . " 
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calculations.  The  basic  issue,  then,  is  the  extension 
of  the  same  empirical  parameterizations  of  eddy  co- 
efficients for  the  reference  atmosphere  to  the  dis- 
persion of  engine  effluents  in  the  perturbed  atmo- 
phere,  an  issue  that  comes  about  from  the  fact  that 
the  assessment  of  aviation  effects  on  the  ozone  column 
involves  an  accurate  calculation  of  the  difference  of 
two  large  numbers,  as  described  in  Section  1. 

• The  2-D  formulation  of  the  dynamics  for  the  perturbed 
upper  troposphere  and  stratosphere  must  not  be  based 
on  the  Reed  and  German  (1965)  formulation  that  invokes 
the  modified  Prandtl  mixing  length  hypothesis,  which 
must  always  assume  a combination  of  countergradient 
and  gradient  components  of  the  eddy  transports  whenever 
a >>  (3  ( Fig . 18).  It  must  be  noted  that  more  modern 
concepts  for  the  formulations  of  the  2-D  dynamics 
involve  statistical  considerations  of  the  eddy  param- 
eters in  Table  17  instead  of  using  parameterizations 
based  on  K coefficients. 

4.3  STATISTICAL-DYNAMICAL  FORMULATION  OF  2-D  DYNAMICS 

The  splitting  of  the  dependent  variables  (\  = u,  v,  w,  T, 

R)  into  zonal  averages  ( [ x D ) and  eddy  components  (x*)  intro- 
duces the  classical,  unsolved  problem  of  closure  of  the  equa- 
tions for  turbulent  motion.  This  closure  problem  comes  about 
because  the  splitting  procedure  cannot,  of  course,  provide 
additional  equations  to  calculate  the  eddy  transports.  For 
example,  the  derivation  of  Eq.  (9)  for  9[R^]/3t  brought  about 
the  double  correlations  [x*R*3  for  the  poleward  and  vertical 
eddy  transports  of  mixing  ratio.  This  first  averaging  proce- 
dure is  also  known,  in  this  case,  as  the  first  moment  of  the 
continuity  equation.  To  obtain  additional  equations  for  the 
double  correlations,  it  is  possible  to  take  the  second  moment 
of  the  continuity  equation  by  multiplying  it  by  R^.  However, 
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this  second  moment  now  generates  triple  correlations,  which 
would,  in  turn,  require  the  third  moment  of  the  continuity 
equation.  Since  this  procedure  would  continue  in  the  same  man- 
ner through  the  successive  higher  moments,  it  is  seen  that  it 
is  not  possible  to  close  the  overall  system  of  equations.  Since 
the  conservation  equations  are  closed,  it  is  a fundamental 
question  to  inquire  how  the  averaging  process  has  managed  to 
convert  them  into  an  open  set.  Batchelor  (1971)  has  provided 
two  complementary  concepts  to  answer  this  question:  (a)  the 
conservation  equations  may,  indeed,  be  considered  to  be  open  in 
time,  because  any  instantaneous  configuration  or  "realization" 
must  be  integrated  over  the  whole  (future)  time  domain  before 
they  can  give  exact  values  for  the  various  correlations;  and 
(b)  the  additional  correlations  from  the  higher  moments  provide 
less  and  less  additional  information , so  that  only  a finite 
number  of  correlations  is  required  to  determine  the  properties 
of  the  turbulence  field  with  some  specified  accuracy . It  Is 
this  second  concept  that  is  of  interest  in  the  statistical- 
dynamical  formulations  of  turbulence;  which  attempts  to  close 
the  equations  of  motion  by  parameterizations  of  the  triple, 
instead  of  the  double,  correlations  as  in  the  Prandtl  mixing 
length  approach  (i.e.,  Eqs . 11a,  b).  The  main  intent  of  such 
statistical  approach  is  to  improve  the  physics  of  turbulent 
models  by  taking  into  account  additional  energy  transfer  pro- 
cesses among  the  eddy  and  mean  flow  components. 

The  history  of  the  application  of  the  statistical  approach 
to  the  troposphere  includes  efforts  by  Smagorinsky  (1964),  who 
calculated  the  eddy  transports  of  momentum  and  heat  by  using 
phenomenological  properties  of  the  Rossby  waves  involved  in 
such  transports.  This  approach  was  able  to  calculate  the  maxi- 
mum values  of  the  transports  of  momentum  and  heat  at  the  correct 
latitudes  for  the  separation  of  the  tropospheric  Hadley  and 
Perrel  cells  near  30°,  as  well  as  the  Ferrel  and  polar  cells 
near  60°  (Fig.  19).  It  is  important  to  emphasize  that  these 
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results  did  not  use  the  modified  Prandtl  mixing  length  hypo- 
thesis (Eqs.  11a,  b)  for  the  poleward  eddy  transports  of  momen- 
tum and  heat,  and  that  the  triple  correlations  in  this  approach 
can  be  shown  to  be  identically  zero  (Hidalgo,  1972).  Other 
examples  of  the  application  of  statistical-dynamical  formula- 
tion to  the  tropospheric  motions  are  those  by  Saltzman  (1968), 
Kurlhara  (1970),  and  more  recently  by  Vernekar  and  Chang  (1978). 

The  2-D  statistical-dynamical  formulation  must  be  based  on 
the  respective  equations  for  the  meridional  circulation  (Table 
16)  and  eddy  motions  (Table  17)-  However,  the  results  in  these 
tables  must  consider  also  the  continuity  equation  for  the  tra- 
cers so  as  to  derive  the  corresponding  terms  for  [R]  and  R*. 

Such  a formulation  is  being  considered  at  GFDL  by  Yen-Huei  Lee 
and  Mahlman  and  it  consists  of  five  equations.  The  first  equa- 
tion is  the  conservation  of  the  mean  tracer  mixing  ratio  (Eq.  9). 
The  next  three  are  second  moments  of  the  momentum  equation  for 
the  poleward  direction,  the  energy  equation,  and  the  continuity 
equation  for  R^,  respectively.  These  equations  yield,  respec- 
tively, the  local  rate  of  change  of  the  northward  eddy  tracer 

flux  [v*R*],  the  vertical  eddy  tracer  flux  [u)*R*],  and  the  local 

2 

rate  of  change  of  the  eddy  mixing  ratio  [R*‘  ].  The  fifth  one 
is  a relationship  for  the  term  [T*R*]  that  appears  in  the  energy 
equation.  Important  points  in  regard  to  this  2-D  statistical- 
dynamics  formulation  are  the  following: 

• It  does  not  require  the  use  of  the  modified  Prandtl 
mixing  length  hypothesis  (Eqs.  11a,  b),  but  must 
utilize  meteorological  know-how  and  experience  based 
on  both  observations  and  3-  model  results. 

• It  will  have  to  establish  parameterizatlons  for  triple 
correlations  involving  eddy  tracer  terms  such  as 
[v»2R»],  [v»w»R»],  etc. 

• The  solution  of  this  system  of  equations  would  require 
inputs  for  the  mean  (i.e.,  [v],  [ut]»  [T])  as  well  as 


double  correlations  (i.e.,  [v#2],  [v#w*],  [v*T»], 

[w*T*],  [ T* 2 ] in  the  stratosphere. 

The  potential  payoff  of  such  statistical-dynamical 
formulations  for  the  reference  atmosphere  would  be  a 
greater  confidence  in  the  validity  of  its  results 
when  extended  to  the  mass  eddy  transports  [x*R*]  in 
the  perturbed  stratosphere.  However,  because  of  the 
difficulty  of  the  endeavor,  the  advent  of  a valid 
statistical-dynamical  formulation  of  the  stratospheric 
dynamics  may  be  a few  years  away,  depending  upon  the 
priority  placed  on  such  effort  by  dynamic  meteorologists 
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Note  added  July  1979  to  Report  No.  FAA-AEE-78-23 , 
Status  of  Representative  Two-Dimensional  Models  of  the 
Stratosphere  and  Troposphere  as  of  Mid- 1 97 8 , 
by  Henry  Hidalgo  (See  pp.  102-105) 


Eqs.  16-18  may  be  valid  only  at  middle  latitudes.  Efforts 

are  currently  being  made  to  calculate  the  relevant  Reed  and 
2 

German  [a  ]/[a]  parameter  from  wind  data  of  a general 
circulation  model.  Results  are  expected  to  be  published 
in  the  near  future. 

Available  results  indicate  that  the  three  2-D  parameteri- 
zations  of  stratospheric  dynamics  (pages  104,  105)  do  not 

* p 2 

always  yield  valid  values  for  [a  ]/[a]  . 

In  view  of  (1)  above,  and  contrary  to  the  statements  on 
pages  104,  105,  the  three  2-D  parameterizations  of  strato- 
spheric dynamics  may,  in  fact,  be  consistent  within  them- 
selves. However,  in  view  of  (2),  the  three  parameteriza- 
tions remain  inconsistent  with  each  other  and  may  also  be 
inconsistent  with  the  basic  Reed  and  German  formulation. 


